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Beiträge zurnumerischen Lösung linearer Eigenwertprobleme. II’) 
Von N. Joachim Lehmann in Dresden 


Die in einer gleichnamigen Arbeit?) begonnenen Untersuchungen über die Eingrenzung vorgeschriebener 
Bigenwerte werden fortgesetzt und abgeschlossen. Von den Ergebnissen der ersten Veröffentlichung aus- 
gehend, wird zuerst ein Verfahren angegeben, das von unten und oben gegen den interessierenden Bigenwert 
konvergierende Näherungen liefert. Eine weitere Konstruktion führt nochmals auf untere Eigenwertschranken. 

The inwestigations of a previous paper of Ihe same, heading?) concerning the limitation of eigenvalues 
of a perassigned order are continued and terminated. Using the results of the first paper, a method is given 

‚ yielding approximations that are converging from below and from above to the fixed eigenvalue. A further 
construction leads to lower boundaries of the eigenvalues likewise. 


Les recherches, commencees dans un traite du m&me en-töte?) et qui concernent les limites d’une valeur 
propre, d’un ordre arbitraire, sont continuees et lerminees. ‚Partant des resultats de la premiere publication on 
arrive d’abord & une methode d’approximations successives, convergentes 4 la valeur propre en question de 
bas en haut et vice versa. Une construction nouvelle nous fournit une autre fois des limites inferieures des 
valeurs propres. 

HeenezoBauna IO BbIYHCJIeHNI ABYCTOPOHHHX OIEHOK COÖCTBEHHBIX yncel, HayaTble B 
OAHOUMEHHOH Pa6oTe ?),, 31eCb IPOAOMKAMTCH U 3araHmyuBamTch. Ha OCHOBaHHuH pesysIbTaTOB 
nepboiü yacru paboTkI cHayama paspabaTsıBaercaH MeTon BEIUUCHeHNA HPHÖALmKeHnK, 
CTPeMAINUXCAH CHUSY U CBEPXy K HCKOMOMy COÖCTBEHHOMy yuciay. JlarıpHeüniee ImocTpoeune 
APHBOAUT ONATb K HISKHeÜ TPAHHIE OMEHKH COÖCTBEHHOTO YHCIA. 


Den Überlegungen und Beweisen in dieser Arbeit liegt wieder die in Gl. (1)* der ersten 

_ Veröffentlichung?) vorgelegte lineare homogene Integralgleichung mit reellem symmetrischen 

Kern mittlerer Stetigkeit zugrunde; den Ansatzpunkt für die beabsichtigte Eigenwertberechnung 

geben die vorangegangenen Untersuchungen über die neue Eigenwertdarstellung, insbesondere 

die Gl. (20)*. Während bisher in den Abschn. 7*, 8* lediglich die mit Gl. (20)* aufgezeigten 

-  Minimum- und Maximum-Eigenschaften verwertet wurden, sollen die Min-Max- und Max-Min- 

Bi» Sätze jetzt ganz allgemein ausgenutzt werden. Damit lassen sich bisher noch bestehende 
E.; Schwierigkeiten für die Eigenwertberechnung beseitigen. 


I. Näherungsfolgen, die von unten ‚und oben gegen einen bestimmten Eigenwert konvergieren 
“ 1. Die Lösung der Men rox. Aufgabe nach Gl. (20)* für eine 


spezielle Funktionenklasse. 
Analog zum Vorgehen von Ritz?) soll jetzt in der Extremalaufgabe nach Gl. (20)* nur eine 
Funktionenklasse {p}; zur Konkurrenz zugelassen sein, die als lineare Mannigfaltigkeit in 
der Form 


k 
WoR= 2m 9%) a nur dh) 


definiert ist. Dabei sind die 9,(x) k voneinander linear unabhängige s-,,‚Koordinatenfunktionen“, 
die 9, beliebige Konstanten. Die @,(x) sollen derart allgemein gewählt werden, daß lim {p}. die 
_ Gesamtheit der s-Funktionen umfaßt. k—>» 
Die bei diesem Funktionenbereich im 7T-Quotienten '7,(y, t;) erreichbaren Extremwerte 
seien mit A;4;(k, ti) bezeichnet: 
en: | Max Min (>#;) Bela, 
x b Min Max (<t;) ZW ir Ass (k; t;) =, 1 , res  -- (2). 
v Wrıla, t 
1) Diese Arbeit gibt die Teile II und III der von der TH-Dresden 1948 genehmigten Dissertation 
Verfassers wieder. j 
2) Abgedruckt in Z. angew. Math. Mech. Bd, 29, S. 341—356. Hinweise auf Gleichungen, Abschnitte 
und Fußnoten aus der ersten Veröffentlichung werden hier durch Anfügen eines Sterachens gekennzeichnet. 
3) Vgl. z.B. Collatz, Eigenwertberechnung (Fußnote #) *). 


a N 
Lehmann, Beiträgs zur ni linearer Eigenwertprobleme., ! 


Die Argumente deuten überäll auf die Abhängigkeit von der Zahl k der Koordinatenfunktionen Ri} Ti 
E und der Größe des Parameters t;. Man berechnet die A:i+;(k, t;), wie anschließend bewiesen wird, 
B als Wurzeln der k-zeiligen Determinante 

hr 0=D,;(ti, A) Er e —t,a” Aare ara ae—(A+#) are +At; R ®) 
7 Ole. 


<e mit % 

Be | = [oy(® 2) [m (0) —% [Kt (2,s) @(s) ds] dx,. | 
er Rn J[K(&, s) 9(s) Pe —u[K (8, 5) 9. (8) ds] dads Ei. (da) 
Er bzw. 

©, = [mia)geia)dz, =/[Kas)miagelsdeds, we STE 
B:; are = [[[K(a, s)9,(s) ds) [| K(&, 5) Pe(s) ds] de = [[ K% (a, 5) Pu(z) pe (s) dads w). 


j 
? Die a”, stellen eine Verallgemeinerung der Schwarzschen Konstanten aus Gl]. (26)* | 
wo dar. Ist mit o,(z 2) = po? x) 0” (2) %) RR X, 8) ) )ds=|K® ( 2, 3)o,(s) ds (n— 1, Ze zB 
: die nach Gl. (25)* n-mal iterierte Koordinatenfunktion, so wird unabhängig von m a E 


A )ya)de=ar? O<m<n=01h2%....... (de) i 


' Wegen der Kernsymmetrie ist Ayg= Ag» Bye = Bo, usw., so daß auch die Determinante b 
D,(t;,, A) symmetrisch wird. Die Lösungen A;+;(k, t;) von Gl. (3) sind alle reell und von &; | 
verschieden, sie werden unter Beachtung ihrer Vielfachheit in der Reihe 


.<Aulbt) <Alkk)< 4< Ani) <Annt)z... 


geordnet und beziffert. Uneigentliche Lösungen 1/A = 0 bleiben immer unberücksichtigt. 
Bisher wurde entsprechend Gl. (12)* ,< t;< A;+ı vorausgesetzt; im Falle 4; —%4, oder 
— Airı kann dabei einmal A,(k,t)=4(=4) bzw. Arrı(k, t;) = t; (= Ai+ı) eintreten, sobald 
nämlich A;(k, t;) bzw. A;tı(k, t;) unabhängig von t; ein Eigenwert ist. fr 
Der Nachweis für die mit Gl. (2) behaupteten Eigenschaften der nach Gl. (3)/(5) bestimmten 
Airj(k, t;) läßt sich auf verschiedenen Wegen führen; um den Zusammenhang mit der be- 
kannten Hauptachsentheorie*) quadratischer Formen möglichst eng zu gestalten, a man 
vielleicht am besten an die Gl. (19)* an. Im dort auftretenden Quotienten 1(T,(y,t) —Lt) 
soll jetzt nur die Funktionsklasse {po} ‚ beachtet werden. Eintragen der Darstellung N GL (1) ee 
ergibt mit Gl. (17)* und den Abkürzungen aus Gl. (4a) 


k 
| 1 > Pr Pe Bro 
ET BE EENEÄNGE > 7a (6) 
A, % ER ER AL RR Rear an. ur LE 
Kr ; a Pr Pe (Ave — ti Be) ; Dr 
‚e= 4 


Die zu der Extremalaufgabe Gl. (19)* gehörigen u We Gl. lm nehmen 
für die Funktionenklasse {p}, die Gestalt an: BR oe 


I rl { 3,0. h. (mit = [% (8) %) 9,(%) ) dx) folrnan- Ip n.=0 } al PRROR Ma) 
' oder bei Vektorschreibweise mit D®= (pı, Pa - - > Pk) und 3 = (Mr Yam..o DE ERRT A: 


v=1 i 

Da die v; die Gesamtheit der s-Funktionen vertreten, können die yYyr bei nd 
Yr+ beliebige Vektoren sein. RS HN a 
-Nun stellt die Gl. (6) den Quotienten zweier quadratischer Formen : 

ist nach Gl. (17)* [s. auch unter Gl. (20)*] mit 4; #4, positiv- -definit 
— A] zumindest immer positiv-semidefinit. Die Max-Min- bzw. Min-Max 
men bei linearen Nebenbedingungen entsprechend Gl. (7) sind aber wohlbekanı 
sie als Nullstellen der symmetrischen Determinante 


Bre—r(An—uB)|=0 neo=1, Du Ba 


Die k-Wurzeln x, (mehrfache werden ihrer Vielfachheit entsprechend oft; gez 
sind reell und können daher der Größe nach geordnet werden, wegen der. Defi 
von Gl. (6) ist, von Grenzfällen bei 1; —=4, abgesehen, immer 1/x, # 
einmal mit der kleinsten Wurzel als x, beginnend in der Reihe 


“4 Courant-Hilbert, a.a.0O., insbes. $. 19 u. nn A S. 32 


my = 0 = Ir r=1,2,..., % oder DW) 2; KIN . a. Bi . 
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N Er U BE 
N 3 RT 
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BR $ x) <H., <a zen. <x rr, [Reihe x, mit 5j= 0, —1,-2,..,—k+1] » (9a) 
und einmal mit der größten Wurzel als x, beginnend nach 
Dr, Sa (Reihe, mit SE 2... KR]. oe. . . „000 (0b) 
erfolgen. 
u Dabei ist selbstverständlich x, = %_4+1, X = *%_x+2 usf., jedoch wird darauf im folgenden 
b kein Bezug genommen. Die Zweckmäßigkeit der Doppelbezifferung wird bald ersichtlich. 
Die nach wachsender Größe geordneten x; (also fürj < 0) können als Folge der Max-Min- 
{ Werte des Quotienten von Gl. (6) angesehen werden; die gleichen Werte nur in umgekehrter 
; Reihe (d.h. die Folge x; bei j > 0) ergeben sich nochmals beim Aufsuchen der Min-Max-Werte 
. des Quotienten [?) insbesondere S. 29]. 
j Es ist bzw. 
Max Min = Max Min 
5 7 Buche, BR P» Pe Bye v {pl 10a; 1 2; 5=0,—1,... 10) 
E 2 nn 5 or] nee ae el re A A ( 
4 Min Max ab Min Max a el, 2N... 0b 
r 5 ae, mp Pe (Ave — bi By.) E (mel), ; J 
| le) je 12. RR 
"Ei - Entsprechend Gl. (11)* gibt a;= als für ) die zum jeweiligen x; ge- 
2 hörige Zahl von Nebenbedingungen (0,1,2,...,(k—1)) nach GI]. (7) an. 
* Wie schon im Mittelteil der Gl. (10a, b) angedeutet ist, sind die x, wegen der Beziehungen 
Gl. (6)/(7) auch die Max-Min- bzw. Min-Max-Werte von 1/(T,(y, t;) —t;), wenn lediglich der 
Funktionenbereich {p}; berücksichtigt wird. Die Darstellung der gesamten x; einmal nach 
h Gl. (10a),einmal nach Gl. (10b) ermöglicht es, alle x, >0 als Min-Max-, die x;< 0 aber als 
N Max-Min-Werte anzusehen. Mit dieser Auffassung führt Gl. (10) auf 
| Min Max(>0)\)__ 1 Rt ie 1 IR an 
Max Min <O), Tv) —u Ay GV=0-—1, tirike FIAT 5 


“ {ee 1a, 
y 
Hierbei werden offensichtlich alle von Null verschiedenen x; erhalten; diese bestimmen auch die 
in Gl. (11) eingeführten Größen Ares ist 


ev > 2 
HR sofern 5; #0, und zu 120 gehört ;=|0, ER el ra (12). 


Beachtet man hierzu noch die Ordnung der z; für 5 >0 (j<0) nach GI. (9b) (Gl. (9a)), ergibt 
sich für A;+; die Größenbeziehung 


für 5 >0 De a un Na ran ı6 (13a) 
für ;<oO TE A EN ER ER OREMPERG N 


Diese Reihen können solange fortgesetzt werden, als noch Werte x; > 0 oder < 0 vorhanden sind. 
Bi Die Extremwertdarstellung Gl]. (11) für die Zahlen A;+;, die aber mit Gl. (12)/(13) auch aus 
den Lösungen x; der Determinante Gl. (8) bestimmt werden können, entspricht nun völlig der 
Eigenwertdarstellung in Gl. (19)*. Zunächst ist hier die Gesamtheit der /p}, zur Bestimmung 
er Min-Max-(Max-Min-)Werte von 1/(7,( t,) — t;) herangezogen worden, benutzt werden 
ber im Min-Max-(Max-Min-)Problem wiederum nur die 2; = 1/(A;+ —t;) >0 ‘< 0); genau 
wie bei Gl. (16)* genügt es also, im Min-Max-(Max-Min-)Problem von vornherein nur solche 
Funktionen aus {p}; zur Konkurrenz zuzulassen, für die bzw. 1/(T,({p}a, t) — ti) >0(< 0) 
ausfällt. (Damit wird zwar die zugelassene Funktionenklasse eingeengt, sie enthält aber immer 
noch alle verwendeten Extremalfunktionen des alten Bereiches und ergibt also wieder die- 
ben Extremwerte.) In diesem Sinne wird nun auch Gl. (11) weiterhin verstanden. 
‚Nach den soeben zuendegeführten Betrachtungen unterscheidet sich das Extremalproblem 
(1) von dem fruheren Gl. (19)* nur durch den Bereich der zugelassenen Funktionen und 

eventuell den Werten und Anzahlen der Extremwerte x; bzw. A:+;. Daher folgt aus 
auch sofort entsprechend dem Übergang von Gl. (19)* zu \20)*: 
Max Min (>t; REN 
nn a Ti) = Al)  ... - Sel 
F 23:59 1a, 

ıst über die x; aus Gl. (8), (9) in Gl. (12), (13) definierten A;+; sind nach Gl. (14) 
rte des 7-Quotienten 7,(y,t) im Funktionenbereich {p}; gefunden; die Ay45 
‚ösungen der Gl. (2), d.h. A45 = Air; (k, t;). Es bleibt nur noch zu zeigen, daß 
Base eiHhingen (3)/(5) damit harmonieren. Vergleich der Ordnungsbedingungen 
ENTER } h F $ h % £ 


er ar 


1? 


„ Ta ; a us 
i 4 TE 
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N Gl. (5) und GI. (13) ergibt Übereinstimmung. Außerdem genügen die Ay, nach GI. (12) und (8) 
4 alle der Gleichung 
Bi 
Br 1 3 
De] 47 


und hiermit ist die Beziehung A; 1; = Ai+j(k, t;) nochmals bestätigt. Man bemerkt, daß zu den 
bei Gl. (12) unberücksichtigt gebliebenen x; = 0, bei Gl. (5) ebenfalls weggelassene uneigentliche 
Lösungen 1/4 = 0 gehören. Die unter Gl. (5) angeschriebenen Grenzbeziehungen für 4 —4; 
oder — A;;ı folgen aus den im Abschn. 4* und unter Gl. (20)* für T, (%, 1) 8 ganz allgemein nach- 
gewiesenen entsprechenden Eigenschaften. i 


2. Die Aıs(k,t) als Schranken für die Eigenwerte und ihre Ab- 
hängigkeit von Ä. 

Nachdem die Eigenwerte A;+; und die Größen A;4+;(k, t;) mit Gl. (20)* und Gl. (2) (resp. 
(14)) als Max-Min- bzw. Min-Max-Werte des gleichen 7-Quotienten in verschiedenen Funktions- 
bereichen nachgewiesen sind, können zu ihrem Vergleich die bekannten Prinzipien der Varia- 
tionsrechnung benutzt werden®). Hier wird zuerst benötigt, daß beim Erweitern (Einengen) 
der zugelassenen Funktionenklasse Max-Min-Werte höchstens fallen (steigen), Min-Max-Werte 
allenfalls steigen (fallen). 

Nach Gl. (2) sind die A;ı;(k, t;) für j > 1 Max-Min-Werte, für j <0 Min-Max-Werte von 
T,(Y, t)) im Bereich der Funktionen {g};; erweitert man die konkurrenzfähige ‚Funktionen- 
klasse durch Hinzunahme einer weiteren Koordinatenfunktion +1 (x) (macht man also einen 
(k +1)-gliedrigen Ansatz entsprechend Gl. (1)), so findet man im Extremalproblem Gl. (2) (also 
als Wurzeln von Dy+1(t, 4) = 0 nach GI. (3)/(5)) als Extremwerte A;4;(k + 1,;) und es ist 
nach dem vorher angeführten Vergleichssatz 

für j > ARE, (k, 7) > A; (k — ık t;) >Arz ehe dee ee (16a) 
für <O Air; (k, &,) <A; (k+1, t,) N FE EN 3 NS (16b). 
Für die letzte Beziehung der Ungleichungsfolge wurde beachtet, daß die Eigenwerte /;;; nach 
Gl. (20)* jeweils im gleichen Extremalproblem erhalten werden, wenn schließlich die Gesamtheit 
der s-Funktionen berücksichtigt wird; wegen der Festsetzung über die Allgemeinheit der 
Koordinatenfunktionen bei Gl.(1) enthält lim fp}, alle s-Funktionen, so daß die Grenzbeziehung _ 
e k>o 


za Ayılk,. Ey on 873 B BE) 4177 


gültig ist. 

Wie üblich wurde zunächst A,< t;< A;ı; angenommen, t; — A; oder — A;+ı sind Co 
als Grenzlälle in allen Betrachtungen mit eingeschlossen. 

Mit den Gl. (3)/(5) und (16)/(17) ist die beabsichtigte Konstruktion von Näherungsfolgen, 
die von unten und oben gegen vorgeschriebene Eigenwerte konvergieren, eigentlich beendet. 
In den nächsten Abschnitten wird noch der Einfluß ‚des Parameters it, betrachtet, anschließend 
werden Fragen über die Anwendung der Methode besprochen. 


3. Die Abhängigkeit der Ay;(k,t) vom Parameter t; 


Bei der Erledigung dieser Aufgabe müssen die Extremwerte verschiedener Extremal- 
ne. zueinander in Beziehung gesetzt werden. Der dazu KloeeseHun Satz sei voran- es 
geste u: 

Werden zwei Ausdrücke 4(y) und B(y) für die Funktionen y? eines elaaran Bereiches _ 
miteinander verglichen und ist dabei immer A(y’) > B(y”), so gilt diese Beziehung auch für 
alle einander entsprechenden zum gleichen Funktiohenbersdh gehörenden Mas und 
Min-Max-Werten von A(y) und B(y): 


aus A(y®) > B(ye) folgt Max Min Max Min Bi 
für alle y* Min Max 4(y) 2Min an) B(W) Se Sr See ei 
er iorz vo yellöa RE = x. 
= 0,1; 


(Die Extremwerte entsprechen sich, wenn die Zahl der N .. Orthogonalitätsbedingı i 
übereinstimmt.) 
Mit diesem Satz läßt sich die Abhängigkeit der A;+;(k, t;) vom Parameter Er leic Ve 
blicken. Zuerst wird einmal der Fall 5 <0 vorgenommen. Die Bestiumngg 3 
Breiht dann 
Min Max (<t) T,(y, t)= A;4;(k, t; =0,—1,— 
eh Yy, 6) + (A, &;) 2 “E 
5) Courant-Hilbert, a.a.O., insbes. 8.353. 


4 1 A 1 le x s x { 
Fr r » Y 
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Diese Min-Max-Werte werden-mit entsprechenden (zu denen also die gleiche Zahl von Ortho- 
gonalitätsbedingungen gehört und die damit die gleiche Kennziffer j tragen) der Aufgabe 


Min Max (759.0 = Ay; (k, i) Bell, —2, en. VEREEL2I) 
v {PL lo}jj| 


verglichen. (Die Verwendung des Index : bei der Bezeichnung der letzten Extremwerte ist durch 
die Beziehung A;ı;(k, > Aiız(k, t;) für — t, gerechtfertigt.) Dabei sei 


De Em . S2IR 
Zur Vereinfachung wird ein Zwischenschritt eingebaut und Gl. (19) zunächst mit 
Min Mx T,y)=A, 53=09,—1—2,...; Nebenbedingung: T,(yt)<t;. . . (22) 


v {orıtji 
in Beziehung gesetzt. 
Wegen der bei Gl. (22) besonders angeschriebenen Nebenbedingung sind in den Problemen 
Gl. (19) und (22) offensichtlich die gleichen Funktionen konkurrenzfähig. Nun gilt nach Gl. (28)* 
und (30)*, sofern nicht unabhängig von 2 1/7, (y, {)—=0 ist, für jede zugelassene Funktion 


»(Y) Tal) —us ()] 
2 UN li ea a) 
1(y 5) = maly) Kaly) —& 


Wie unter Gl. (30)* (man beachte hierbei auch das Bild 1*) folgt damit aus der zu Gl. (19) und 
(22) gehörigen Bedingung 7, (y, t;) <t,; sofort u, (y) < t;, beit; < t ergibt daher Gl. (23) immer 
rt < T,(yrt): Nach ae am Anfang dieses Abschnitts genanaten Satz ist somit 
Milk, TE Geht man jetzt von der” Gl. (22) zu (20), so wird die Nebenbedingung 
T,(yt;) <t; durch T,(y, tt) < t abgelöst. Das bedeutet aber höchstens eine Erweiterung des 
zugelassenen Finktionenbereich: (man vergleiche hier wieder Gl. (23) oder einfacher die Folge- 
rungen unter G1.(30)* bzw. Bild 1*), denn für jede Funktion y, die 7, (y, 1)< t; ergibt, istu,(y)<;, 
nach Gl. (21) also auch u3(y)< t und damit wiederum 7,(9, )< i. Beim Erweitern einer 
Funktionsklasse können die Min-Max-Werte aber allenfalls wachsen (vgl. im vorigen Abschnitt), 
d.h. A =< Alk, D), womit 


Artzj(k, t;) < Airj(k, ) für'..3 =Z0 bei ESS TR A a A ee (24) 


miteas(y Lu, (a B3ily)) >07 (23). 


gezeigt ist. 

Wie der Beweis zu Gl. (24) vorführt, existiert mit A;+;(k, t;) auch A;+;(k, t). Geht man 
umgekehrt von A;+;(k, t) aus und verkleinert den Parameter £ bis t;, so hat der Satz Gl. (24) 
nur dann Sinn, wenn die Existenz von A;+;(%, t;) behauptet werden kann. Man braucht dazu 
nur zu bedenken, daß beim Schritt von Gl. (20) zu (22) der konkurrenzfähige Funktionenbereich 
eventuell so weit eingeengt werden muß, daß {ot für Gl. (22) überhaupt keine Funktion mehr 

stellen kann. (Das läßt sich beim eingliedrigen Ansatz an Hand des Bildes 1* sehr bequem 
verfolgen.) 

Zum Vergleich von Re t;) und A;yj(k, £) beij > 0 und t; < t bleibt die Gl. (24) gültig, 
wenn man nur die Existenz von A;+;(k, t) voraussetzen kann. Der Beweis wird durch schritt- 
weisen Vergleich der drei Max-Min-Probleme 

Max Min (>) 7,9, ) = Aitj(k, t) a REN + BULORUFEE AN 5 Ser ar la Me (25) 

®» {la 
Max = Min ST, ty, Az j=1,;2,...; Nebenbedingung: T,(y, t) >t, (26) 

v {Pr 1 WW) 
Maxs Min St), Tut) Ar; lk, t;) Ge N rl Pen 9ck(27) 

a v {Ar Wu-ı) 
geführt. In den ersten beiden Aufgaben sind sichtlich die gleichen Funktionen konkurrenzfähig. 
- Wie bei Gl. (23) und (30)* usw. folgt dabei für jede zugelassene Funktion aus T,(y, t) >t zu- 
nächst u,(y) > tund mit £ > t; über u, (Y) >t; schließlich 7T,(y,t) 2 T,(y, &;). Der Vergleichs- 
satz vom Anfang dieses Abschnitts ergibt daher Ajyj(k, t) > Aiıj. Der Übergang von Gl. (26) 
zu (27) ändert nur die Nebenbedingungen, die zugelassene Funktionenklasse wird höchstens 
vergrößert: Jede in Gl. (26) konkurrenzfähige Funktion erfüllt 7, (9, £) > t, woraus entsprechend 
vo) hal y) >t >t; und schließlich 7, (%, t;) > t; folgt — damit ist sie aber auch in Gl. (27) 
Ben, BuSnEnell ‚neben anderen zusätzlichen Suablionen, Da Be Eryeitern des Funk- 


an need ud t>n..t.n.. 020 (28). 


Ei 
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Für die Untersuchungen dieses Abschnitts wurde eigentlich weder für noch £ die Zwischen- 
lage A;< t,t< Aiyı benötigt. Soll aber Arzz (k, t) eine Ai4z(k, &;) entsprechende Schranke für 
den Eigenwert A;ı; sein, muß selbstverständlich auch A;< t< A;+ı angenommen werden. i 

Für die praktische Anwendung der A:ız(k, t;), die nach Gl. (16) für) >0 obere, fürggs 0 
untere Schranken für die Eigenwerte ergeben, entnimmt man Gl. (24) und (28), daß zur Brech- 
nung der oberen Schranke möglichst kleine, für die untere möglichst große Parameterwerte aus U 2 
m dem Intervall ;...... sr zweckmäßig sind. Anderseits liegt ein sroßer Vorteil der Methode end j 
7 gerade darin, daß die Näherungswerte A:+;(k, £;) zufolge ihrer Extremwerleigenschaften beson- .- 
: ders bei mehrgliedrigen Ansätzen in weiten Bereichen nur schwach vom Parameter i; beein- | 

flußt werden. Man vergleicht hierzu die Beispiele aus dem Abschn. 6, insbesondere B mit Bild 6. 


>. 4.Die VerallgemeinerungdesEinschließungssatzes vonKryloff. 
Beim Vergleich von Gl. (16) mit der Ordnungsbeziehung Gl. (5) stellt man fest, daß zwischen 
.. jedem beliebigen Parameterwert £ (der untere Index « kann wegbleiben, solange man sich auf n. 
Bi; keinen bestimmten i-ten Eigenwert festlegen will) und den zugehörigen Lösungen A,(k, t) der y 
u Determinante D,(t, A)=0 nach GI. (3) zumindest ein Eigenwert liegen muß — gelegentlich 
kann dieser mit einem A,(k, t) übereinstimmen. Damit liegt der Gedanke nahe, genau wie m 
BE ®. Abschn. 8* den Parameter t—= ty zu bestimmen, für den die durch ty und An(k, tar) für einen 

beliebigen Eigenwert gelieferten Schranken möglichst eng werden. Eine notwendige Bedingungs- 

gleichung hierfür ist entsprechend Gl. (42)* 


Be. dA(k, t)jdt=1 Rn 9 . (29): 

' Diese Beziehung tritt neben die Bestimmungsgleichung (3) für A,(k, t): D.(t, A) = 0. Differen- 
tiiert man D,(t, A) = O nach t und beachtet Gl. (29), dann erhält man zur Bestimmung von ty 
das handlichere Gleichungssystem: 


0D4(t, A) , 2 Dil, A) | Saat 
+ 0 er 


D.(t,A)= 0 und 


Zufolge der Symmetrie dieses Systems in tund A [vgl. Gl. (3)] muß die Reihe der Lösungen fg; 
gleichzeitig die Werte Ay (k, iy,) enthalten, für die [Ay(k,ty,)—!n;|” =Extr. wird. Benutzt Ki 
können alle reellen Wertepaare tyr,, An(k, iy;) werden, deren Schrankenweite [An (k, u) — tm] a 
brauchbar scheint. Für die Praxis ist wichtig, daß man mit den iy aus D;(ty, A)=0 ein 
System von Lösungen A,(k, ty) erhält, welches mit ty an Stelle &;, nach Gl. (5) geordnet 
wird, und für das damit wiederum die Beziehungen nach Gl. (16) verwendet werden können. > 
Meist genügt die Ermittelung der Lösungen A,(%, ty) für den größten und kleinsten der iy, 2) 
um die günstigsten Schranken festzustellen. Be; 
Bei der Lösung der Gl. (30) wird die Symmetrie ausgenutzt und zunächst £+A=2,td=y_ 
gesetzt. Elimination von z. B. zergibt dann als Resultante allgemein ein Polynom k2-ten Grades 
in y, d.h. es gibt %® Lösungspaare &,, y; für das Gleichungssystem. Aus =t+4A, „=td 


erhält man alle zusammengehörenden tyr,, Ayı (k, tar,). | 24 
R ‚.  Sehwierigkeiten bei der Zuordnung der Schranken zu bestimmten Eigenwerten lassen. Rh 
Ei. sich wie bei der Anwendung des einfachen Kr yloffschen Satzes im Abschn. 11*, Beispiel ER 


vermeiden, wenn man von einem geeigneten Nachbarproblem ausgeht. A 


Alle diese Fragen werden am praktisch wichtigen Fall eines zweigliedrigen Ansatzes im 
Abschn. 6, unter C und an Hand von Bild 5 demonstriert, dort wird auch ein einfaches Näherungs- 
verfahren zur Bestimmung der tr angegeben. Hinweise auf den Anwendungsbereich, Genauig- 
keit usw. finden sich im nächsten Abschnitt. RR 7 00 £ 


D. Bemerkungen zurÄAnwendungder vorgeführten Sätze. 2 
Den einfachsten Überblick über die Anwendungsmöglichkeiten der im 2. Abschn. erhaltene 
Beziehungen gestattet die schematische Darstellung dieser Ergebnisse in Bild 1. 0 eit 


Bo km Bi Ar 0 Ara 


| Bene | 
TE RES | 
Aılktl,6) Alc+l,6) °  Analetl,u) - Aualet 
A-ılk,bi) Alk) . Arılk, bi) 5 


Bild 1. Die A; ,5(k, t;) als Bigenwertschranken 


32 


ein Parameter 4; zwischen A, und A;+1 bekannt, so lassen AR 
nach Gl. (3)/(5) und (16) aus D;(t;, A) = 0 sofort obere 
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untere für A,< £; berechnen. Je mehr Koordinatenfunktionen dabei verwendet werden, desto 
genauer muß nach Gl. (16) das Ergebnis ausfallen. Im Abschn. 3 wurde überdies gezeigt, daß 
innerhalb der durch A; und A;+ı gegebenen Grenzen für £; zur Berechnung oberer Schranken ein 
möglichst kleines, zur Ermittelung unterer Schranken umgekehrt ein möglichst großes t; am 
günstigsten ist. Falls einmal A; der größte oder A;ıı der kleinste Eigenwert der vorgelegten 
Aufgabe sein sollte, wird man £;— + oo bzw. — — oo wählen; dabei soll unter D,;(, A) das 
Ergebnis des Grenzüberganges 

u 1 See. a. 2 
verstanden sein. Die Determinante Gl. (3) zeigt, daß D,(®, A) = 0 auf das gewöhnliche 
Ritz-Galerkinsche Verfahren zurückführt. 


Soll irgendein Eigenwert 4; eingegrenzt werden, so sind zumindest zwei A; einschließende 
Parameterwerte t, &<i und 15 ß >i erforderlich. Beieingliedrigen Ansätzen wird man damit auf 
den in der ersten Veröffentlichung ausführlich diskutierten Einschließungssatz kommen. Bei 
mehrgliedrigen Ansätzen brauchen nun /, und ts nicht mehr unbedingt den an 4; angrenzenden 
Eigenwertintervallen entnommen zu werden. Bild 2 und 3 zeigen eine schematische Darstellung 
der Verhältnisse beim ein- und zweigliedrigen Ansatz. 


De ki ti; Airı 


Ai-ıl2,ti-ı) Ai-ı(2, ti) Ai2,t) Ail2, ti-ı) 


ha Rt n ti je DL FERN je Ai Aisı 
Al,t) All, ti-ı) Ai-ı(2, ti-a) Ail2, ti-2) 
Bild 2. Die Eigenwerteingrenzung mit Bild 3. Die Eigenwerteingrenzung bei zweigliedrigem Ansatz; 
eingliedrigem Ansatz, Benutzung verschiedener Parameter. 


Spezielle und ins einzelne gehende Anweisungen zur Handhabung der entwickelten Methoden 
können nicht gegeben werden, es läßt sich nicht allgemein vorhersagen, wie im Einzelfall möglichst 
zweckmäßig vorgegangen werden muß. Das hängt davon ab, für welche Eigenwerte in der 
Nachbarschaft von A; man möglichst bequem wenigstens grobe Schranken erhalten kann — 
damit sind die erforderlichen Parameterwerte i, und ?g bestimmt. Zur Aufstellung von 
| D,(t, A) =0 sind dann noch geeignete Koordinatenfunktionen 9, aufzusuchen, ihre Wahl 
j entscheidet wesentlich über den bei der gesamten Rechnung erforderlichen Aufwand. Über prak- 

- tische Fragen zur Bestimmung der Parameterwerte 2; und der Koordinatenfunktionen vergleiche 
man Abschn. 10* und die dort angegebenen Literaturstellen, dazu die Beispiele in Abschn. 11* 
und 6. 

Während die bisherigen Betrachtungen in dieser Arbeit immer an den 7-Quotienten 
A - T,(y,t) anknüpfen, darf nach den Beweisen in Abschn. 5*, 6* auch von 73n+1(9%, t) und bei 
4 positiv-definiten Kernen sogar allgemein von T,(y,t) n=1,2,... ausgegangen werden. Die 
zu T„(y, t) gehörende, D,(t, A) = 0 aus Gl. (3)/(4) entsprechende Gleichung lautet (analog zu 
Gl. (32)*/(33)*) 
Deu, A)= ar, —tar Ale —tar.,)=0 r=1,2.... ... (22). 
Es ist Do (t, A)=D,;(t, A); die a)? sind wieder die verallgemeinerten Schwarzschen Kon- 
 _ stanten von Gl. (4c). 
R- Die aus Gl. (32) erhaltenen A(® (k,?) approximieren bei n— = im allgemeinen unabhängig 


von t die entsprechende Anzahl der absolut genommenen kleinsten Eigenwerte der Aufgabe. 
j Man beachte hierzu die Gl. (39)*/(41)* und die Min-Max- bzw. Max-Min-Eigenschaften der AP; 


letztere verhindern, daß eiw? immer alle A nach dem absolut genompienen. kleinsten Eigen- 


| En; den a ehibereich der entwickelten Methoden unterrichten die Überlegungen 
‚Abschn. 9*. Damit kann das Verfahren z. B. sofort bei Randwertaufgaben vom Typ der 
47 


rel ae = für 41 2,... se. (83 


‘ # ’ u 
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benutzt werden. Geht man von Vergleichsfunktionen g@®(2) »=1,2,...,k aus (die, falls 
7 = 0 kein Eigenwert ist, lediglich U; (EP) —=0(0 und L(g®) — (s-Funktion) erfüllen müssen — 
vgl. bei Gl. (48)*), so können diese wie angegeben alseinmal iterierte Koordinatenfunktionen auf- 


gefaßt werden, es ist nach Gl. (52)* L(p") BL op. Die Determinante Dy(t, A) bekommt 
die Gestalt | 
L'pd) Lio® 
Mi Er ar ara [EI det Zu ornee = DO (34). 
N 0%) | 


Damit ist noch einmal dargelegt, wie die entwickelte Methode für weite Klassen linearer 

selbstadjungierter Eigenwertaufgaben die Eingrenzung beliebiger Eigenwerte mit gewünschter 
h Genauigkeit bei erträglichem Rechenaufwand ermöglicht. Um z. B. Schranken für den 10. Eigen- 
wert eines positiv-definiten Problems zu bestimmen, benötigt man jetzt nur einige grobe Schran- 
ken für die Eigenwerte bei A,, und kann die Aufgabe mit einem etwa 1—4gliedrigen Ansatz nach 
k. Gl. (1)—(5) erledigen; beim üblichen Ritzschen Verfahren wäre ein mehr als zehngliedriger 
Ansatz nötig, der zudem erst eine obere Schranke ergibt. Weiter gestatten die Sätze die Unter- 
suchung der Eigenwertverteilung in der Umgebung eines beliebigen Wertes Z (dessen genaue 
h Numerierung eventuell nicht bekannt ist). Ein entsprechender mehrgliedriger Ansatz gibt 
Ca Schranken für die benachbarten Eigenwerte. Weitere Anwendungsmöglichkeiten brauchen 
nicht diskutiert zu werden, sie ergeben sich von Fall zu Fall. Man beachte nur die in den Bildern 1 
bis 3 dargestellten prinzipiellen Leistungsmöglichkeiten der Methoden und die im nächsten 
ee Abschnitt durchgerechneten einfachen Beispiele. 
| Das Anwendungsgebiet für den erweiterten Kr ylo ff schen Einschließungssatz im Abschn. 4 
ist begrenzter. Er hat den Vorteil, daß mit seiner Hilfe Eigenwertschranken ohne voraus- 
gehende Bestimmung eines Parameterwertes Zi, erhalten werden können, weist aber auch Nach- 
teile auf: Die Rechnung ist recht umständlich (durch ein Näherungsverfahren in Abschn. 6, 
Beispiel C wird sie erträglicher), die Zuordnung der Schranken zu bestimmten Eigenwerten kann 
Schwierigkeiten bereiten und für die erreichbare Genauigkeit gilt sinngemäß die Schlußbemerkung 
von Abschn. 8*, obwohl natürlich die Hinzunahme jeder Koordinatenfunktion die Schranken 
höchstens einengt und bei k— o die Konvergenz gesichert ist. Aus diesen Gründen wird das 
Verfahren wohl allenfalls mit zwei- oder dreigliedrigen Ansätzen und entsprechend dem ursprüng- 
lichen Satz von Kryloff lediglich zur Bestimmung erster Eigenwertschranken herangezogen 
werden, die dann weiter in den Gl. (3)/(5) als Parameterwerte Z, auszunutzen sind. (Wie beim 
Beispiel A in Abschn. 11* wird man immer von einem geeigneten Nachbarproblem ausgehen.) 


6. Beirspiieke: 


A. Als erstes wird das bereits im Abschn. 11* mit einem eingliedrigen Ansatz behandelte 
Sturm-Liouvillesche Randwertproblem Kar 


yY'+(eRa+N)y=0; YO SEIN ee 
in Angriff genommen. Mit den Funktionspaaren 9 = sinz x, gV) =sin2rx bzw. yD = 
sin 322, y) =sin4r x als Vergleichsfunktionen sind die ersten vier Eigenwerte in Abhängig- 


keit vom Parameter ein Schranken einzuschließen; es werden dabei nur zwei zweireihige Deter- 
minanten vom Typ der G]. (34) auftreten. Mit {/ı® = r und A/n? = ö lauten diese 


1 1 us. i 
I-e+ (4-5) -1(7+5)(1-5) #70 rs. | " 
Bear 2 Im: = 
Im (d-8- (r+)) 16-4e+8(4-54)-(r+0) (1-2) +75 ; BR 
ROSN EL Bier ERS | M R 
er 81-9e+e (-45)-179 (9-4) +70 np E-(T+ ö)) i { N Er 
Bl) 6 - en "a ‚0 486h), 
1972 (29 € (T+ ö)) 256-16e+[4-,45)- (r+0)[16-5)+ 70 A ner 


(Der Strich über D,(t, A) deutet auf die benutzten Funktionen yD, yiD.) 


Zur Berechnung der oberen Schranken für die beiden ersten Eigenwerte wird benutz 
daß die Aufgabe höchstens endlich viele negative Eigenwerte hat und daher tu, — oo werde 
kann; dabei sei hier wieder der kleinste Eigenwert ohne Rücksicht auf das Vorzeichen mit A 
bezeichnet. Die damit aus Gl. (36a) erhaltene obere Schranke A,(2, — ©) für A,(e) wird 
Gl. (36b) als t, eingetragen und ergibt A,(2, t,) und A,(2, t,) als obere Schranken für Az(e 
A,(e). Zur Ermittelung unterer Schranken wird in (36b) zunächst die mit Hilfe des Kryloff 

+ F is 


x Y D, (1, A) RR 


Z.angew. Math. Mech. 


'Bd.30 Nr. 1/2 Jan /Febr.1950 Lehmann, Beiträge zur Lösung linearer Eigenwertprobleme. II ) 


Einschließungssatzes in G1.(57)* (die nn für das spezielle Problem stehen unter Gl. (59) *) 


gefundene untere Schranke (5- 5 syı- 
4:2, 4»4,(2 
nach kann t, = A,(2, t,) in Gl. (36a) benutzt 
werden und liefert mit A,(2, t,) und A,(2, t,) 
die fehlende untere Eingrenzung für A,(e) und 
A,(e). Die Ergebnisse dieser Rechnung sind 
im Bild 4 dargestellt (ausgezogene Linien) und 
mit den Resultaten beim eingliedrigen Ansatz 
in Abschn. 11* (vgl. Bild 3*) verglichen (ge- 
strichelte Linien). Obwohl die Durchführung 
des zweigliedrigen Ansatzes im Vergleich zur 
Rechnung in Abschn. 11*, A nur wenig Mehr- 
arbeit erfordert, ist die Güte der Schranken 
doch wesentlich verbessert; die Ergebnisse 
sind in erstaunlich weitem e-Bereich brauchbar. 

Bei :zweigliedrigen Ansätzen sind nun 
aber die Eigenwertnäherungen A;+;(2, t;), die 
nach Gl. (2) als einfache Minimal- (oder Maxi- 
mal)-Werte des T7-Quotienten- 7,(y, t;) (also 
ohne Nebenbedingungen) erhalten werden (es 
sind immer die t,; benachbarten) gegenüber 
den anderen oft bevorzugt. Die übrigen 
A;+;(2, 1;) sollen nämlich nach Gl. (16) A;rs 
oder A;_ı als Max-Min- oder Min-Max-Wert 
annähern, stellen aber im Bereich von {p}, 
analog zu der Doppeldeutigkeit der x, in 
Gl. (10) auch das Maximum bzw. Minimum 
von T,(%, t;) und somit im allgemeinen eine 


ungünstige Näherung für einen Max-Min- 


bzw. Min-Max-Wert dar. 


Nach dieser Überlegung darf erwartet A 


werden, daß die mit A,(2, — x) gegebene 
obere Schranke für A, (e) und die durch A; (2, t,) 
gestellte untere Schranke für A,(e) sogar 
mit einem zweigliedrigen Ansatz noch ver- 
bessert werden kann. Unter Spy Ei 
der Abhängigkeit der A:;+;(2, 1) vom Para- 
meter £; [vgl. Gl. (24) u 


x 


48 € 


55 (le (dt) 


ER E| 
Bl 2 
16-4e+e (3 5% jr (r+öjlı Fr 1) 


1 \ 
81-9e+@ |} A 2le (T+0) (9 5) +TÖ 


2 
nz 


.) <A,(e) als t, eingesetzt, die Lösungen 


‚t,) sind nach Gl. (16) die gesuchten unteren Schranken für A,(e) und A,(e). Da- 


Yr 


ac 


Bild 4. Schranken für die ersten vier Eigenwerte der Aufgabe 
Gl. (35) 


(28)] scheint die Bildung von D (6, A) = 0 mit den Vergleichsfunk- 
tionen VO =sin2r x, iD — = sin 37 x zweckmäßig. 


Es wird 


48 € 


De ‚(138 -(T+0)) 


=0 (36). 


(Das Zeichen — deutet hier auf die verwendeten Funktionen Z, und £,.) Mit {, = A, (2, — %®) 


; h und {, = 4,(2,t,) erhält man daraus A,(2, t,) und Az(2, 
Bi Bänke für A,(e) und Ay(e). 


t,) als günstigste obere bzw untere 


(In Bild 4 strichpünktiert.) - 
ir Nach dieser Verbesserung werden besonders die ersten drei Eigenwerte gleichmäßig gut 
Bu: 


Die Abhängigkeit der A(k, t) von t bei einem mehrgliedrigen Ansatz zeigt Bild5 (S. 11), 


ie & ist einmal der Verlauf der von Gl. Kun für € = 8 definierten nen Kurve as 


= Tür on 172 


yo)=yl)=0 mit el tin, HD = a < ar 


. Dabei wurden mit {,— — o zuerst die oberen Schranken ausgerechnet, so daß dann wegen 


Bild 6 finden sich im nächsten Abschnitt. j 


‚sten Schranken ty ....Asr(2, tar) führen, alle in der Reihe der a 
En? r a2 5 3 
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Ku 45, N ra 
Mit den Nebenproblemen y’+4° y=0 und y’+24y=0, y(0) = y(1) = 0 liefert a} 
der bekannte Vergleichssatz®) zunächst die erforderlichen groben Eigenwertschranken 
m >h>1Rm . IM ER >45n2 
m >, >2m: 624 28". 


Diese Werte genügen zur Festlegung der notwendigen Parameter t;. Sr ; Br, * 
Einfache Vergleichsfunktionen sind die Eigenfunktionen des „Nachbarproblems“ 

y"+Ay=0, y()=y(l)=0; mit W’=sinza, y”=sin2rz, p) =sin3rz wird 

D,(t, 4) = 0, wo wieder t{/®= r, A/n? = Ö gesetzt ist: | 


N=Rm>ı ZeEmR2—=A also 


3—(T+d)4+6rÖ (2— ro) 16/30 | a0 [" ARE 
D,A)=| 2—rö)16B8n 8—(T+N)16+6rÖ (18 — rö)16/5r |=0 (38). | 
f) (18_70)165m 243 —(T+6)36 +66 4 


Für t,> — » und t, = 4,572 liefert die in Gl. (38) hervorgehobene zweireihige Deter- 
minante folgende Schranken und Mittelwerte “ 
0,646? <A, <0,6502 A, = 0,6487? +0,30], B. 

2,77 m <h,<2N7 mM = 287 n® + 3,50], a ; 

Die Mitnahme der dritten Koordinatenfunktion, also Benutzung von D,(t, A) = 0 verbessert die er E 
Ergebnisse wesentlich, es folgt » 


0,6470 7? <A, < 0,6499 }ı = 0,6485 2? + 0,20, 
2,837 22 < 1, < 2,890 n2 771, — 2,86417220,977: 2 
4,36, 02 <A, <ı6,37 8 


}, < 6,3702 und 87° <A, (s. unter Gl. (37)) {5 = 8m? genommen werden durfte. Trotz des ein- 
fachen Ansatzes — das herangezogene ‚„‚Nachbarproblem‘‘ weicht sehr stark vom betrachteten 
ab — sind die Ergebnisse ausgezeichnet; soll auch der dritte Eigenwert eng eingegrenzt werden, 
wäre zumindest noch eine weitere Koordinatenfunktion hinzuzunehmen — evtl. kann dabei 
die erste oder zweite weggelassen werden. Wählt man die Koordinatenfunktionen im vorliegen- 
den Fall einmal nach physikalischen Gesichtspunkten, kann man mit noch einfacheren Ansätzen 
bessere Resultate erzielen. — Weitere Bemerkungen zu diesem Beispiel und ein zugehöriges 


C. Als Anwendungsbeispiele für den erweiterten Kryloffschen Satz nach Ab- 
schnitt 4 sollen die unter A. und B. behandelten Aufgaben dienen. In der Gl. (35) wird speziell 
€ = 8 gesetzt. Zur Aufstellung von Schranken nach Gl. (30) werden die Vergleichsfunktionen 
pP) = sinz x, 9) = sin 27 x herangezogen, das Gleichungssystem (30) wird dann nach Gl. (35a), 
wo e=8 sein soll, mit t{® = r, Al = ö y j 


33 32m +38) tr BB TH sg 
—128(3+ + ö)/9n? 16/3 —8/®+rö |” 
und (schon zusammengefaßt) 


6+ +5 —18(3 +7 +6)/Im| 143/73 —32/e + 3(t Hd) Hr6 —256/972 
— 236/90 1688476 |T| —1288+7+ 89m +8 


Von den 8 Lösungen dieses Systems ergeben natürlich nicht alle ein Minimum der Schranken- : 
weile £....A(k, t), es sind darin auch Maximalschranken enthalten. Zunächst liefert die Auf- 
lösung von Gl. (39) vier brauchbare Wertepaare ty, Ay(2, tu); in Gl. (39a) findet man aus _ 
D (tur, A) =0 alle zugehörigen A,(k, tr). Die Zuordnung der erhaltenen Schranken zu be- 
stimmten Eigenwerten bereitet keine Schwierigkeiten, wenn man z.B. die Ergebnisse von Ab- 
schnitt 11*, Beispiel A, insbesondere bei Gl. (57)*, berücksichtigt — wie dort werden auch 
die Eigenwerte ohne Rücksicht auf ihr Vorzeichen einfach der Größe nach geordnet und nı 


D,(t, 4) -|' 


0... 
a 


=0 ob 


riert. Man erhält so | 


u = = 4,14 <ı, <A, tr = t)=—-3,0090°.. .. 
A,2, tu, =t)= 4,134 <i, < tu) tip — 3,0 
Al, te, =t)=-3679<ı, .....mmtb =-0,8680 << 

tn, =t)=-36WR< .... Art t)=— 0,8680? <i 


In Übereinstimmung mit den Betrachtungen an Gl. (30) sind die We 


Fr 


enthalten. | 
.. %) Collatz, Eigenwertberechnung, a. a. O., S. 134. 
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Als günstigste Schranken findet man im Gesamtlösungssystem: 
— 3,676 <A, <— 3,0097; —0,868n? <A, < 1,094 n8. 


| h Obwohl diese Schranken schlecht sind, stellen sie gegenüber den aus dem einfachen Kr y- 
_ loffschen Satz mit yD=sinaz x und yÜ) = sin 2x x erhaltenen Werten —4, 4m?<A, < 0,67? 
und —2,1n? <}, <2,1r, eine erhebliche Verbesserung dar. 
| Einen anschaulichen Einblick in den Mechanismus dieser Eigenwerteingrenzung gibt Bild’5. 
Dort ist der Verlauf der durch Gl. (39a) definierten algebraischen Kurve A(2, t) aufgezeichnet. 
me Als charakteristisch fallen die Symmetriegerade A= t und die zwei Asymptotenkreuze auf, 
. Die hier interessierenden Parameterwerte /y sind dabei die der Symmetriegeraden £—= A nächst- 
gelegenen Kurvenpunkte, wie man aus den Bedingungsgleichungen (29)/(30) entnimmt. 
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Bild 5. Die durch Gl. (39a) definierte algebraische Funktion A (2,t) mit nach Abschn. 4 bzw. 6 konstruierten Schranken. 


Die bisher vorgeführte Berechnung der /y war recht mühsam. Es scheint daher wichtig, 
daß sich für diese Parameterwerte sehr leicht Näherungen finden lassen. Dazu betrachte man 
nochmals die oben aus Bild 5 entnommenen Charakteristika des Verlaufs von A(2, t); diese sind 
nämlich für die Gesamtheit aller nach Gl. (3) durch D;(t, A) = 0 bestimmten algebraischen 
Kurven typisch. Die Symmetrie in 2 und A zeigt schon Gl. (3), während D;(, A) = 0 [vgl. 
Gl. (3) u. (31)] im allgemeinen %k Asymptotenkreuze liefert. Außerdem kann die algebraische 
Br Kurve Ak 0) die Symmetriegerade = A nur dann einmal schneiden, wenn im Ansatz Gl. (1), 

also in {p} eine Eigenfunktion enthalten ist. Aus A,(k, t) = t folgt nämlich nach der Bemerkung 
unter Gl. (20)*, die anzuwenden ist, weil nach Gl. (2) A,(k, £) immer als spezieller Wert eines 
T-Quotienten 7,(y, it) betrachtet werden kann, daß {= 4, ein Eigenwert sein muß, der von 
 T,(y,t) nur für eine Eigenfunktion angenommen wird. 
Je besser die im Ansatz Gl. (1) enthaltenen Eigenfunktionsnäherungen sind, desto enger 
_ werden die nach Gl. (16) bzw. Bild 1 möglichen Eigenwertschranken — d.h., um so besser 
 schmiegen sich die Kurven A(k, £) ihren Asymptoten an. Falls einmal die Eigenfunktionen als 
oordinatenfunktionen benutzt werden: %, = p,, erhält die Determinante D,(t, A) = 0 wegen 
der Orthogonalität der ‚Eigenfunktionen Diagonalgestalt 


0 faye—(6+ A)ave+ Atare | = IT rlı-3)(1-4) = I 1 Re 


in N Y 
zung des Kronekersymbols ” ist doch 
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A,(k, t) zerfällt also unter diesen Umständen in ein Geradensystem. Das gleiche Verhalten zeigt 
nun nach Gl. (17) bei k>® A,(k, t) ganz allgemein, die Kurve neigt immer zum Zerfallen. 

Und gerade diese Tatsache läßt sich zur Ermittlung von Näherungswerten für iyr aus- 
nutzen. Die Stellen minimaler Schrankenweite liegen annähernd auf der Geraden t-A= 
2 A,(k, ©). Dabei sind die A,(k, ®©) die Koordinaten der Asymptoten. In D;(t, A) = ein- 
getragen: 

D,(t, 2 A, 00, — Ü) = Oar - e.  e: 
gibt dieser Ansatz mit den A,(k, ©) nächst gelegenen Wurzeln erste Näherungswerte ty für ty. 
Setzt man in Gl. (41) zur Abkürzung (2Av(k, oo) —t) t = u, so so bekommt man an Stelle des 
exakten Systems Gl. (30) hier eine Gleichung von nur k-tem Grade, die entsprechend einfacher 
zu lösen ist. 

Die Werte t4, können sofort an die Stelle der ty treten; in D;(tir, A) = 0 eingesetzt er- 
\ geben sie alle zugehörigen A,(k, tjr). Wieder muß zwischen tjr..... A,(k, ty) mindestens ein 
Eigenwert liegen, die Gl. (16) bleiben in Kraft. Gegenüber der Verwendung von iy wird lediglich 
die Schrankengüte verschlechtert, jedoch oft nur so wenig, daß die Berechnung der genauen iy 
überhaupt nicht lohnt. a 
Im Beispiel Gl. (39) erhält man nach Gl. (41) 


u, = 41437? im, = — 3,018 7° iu, = — 0,856 7? im, = 3,016 
(—4,134 72?) (— 3,009 7?) (—0,868 2?) " (2,004. 7°) . 
Zum Vergleich sind in den Klammern die oben festgestellten genauen Werte iyr angegeben. 


Die Abweichungen sind so klein, daß sie in Bild 5 gar nicht berücksichtigt werden konnten. 
Mit der Näherungs- 

VAR formel (41) wird der für 
mehrgliedrige Ansätze 
erweiterte Kryloff- 
sche Satz für nume- 
rische Rechnungen dis- 
kutabel, über sein An- 
wendungsgebiet vgl. am 
Schluß von Abschn. 5. 
Tee RN Das Bild 6stelltden 
Verlauf der aus Gl. (38) 
(bzw. der hervorgeho- 
benen _Unterdetermi- 
nante) bestimmten, zur 
Aufgabe Gl. (37) gehö- 
rigen A (2, t) dar; die ty 
und die zugehörigen 
A,(2, tu) sind mit ein- 
gezeichnet. Zufällig zer- 
fällt diese algebraische 
Kurve in zwei von zwei- En. 
Be. | A $ , ter Ordnung, nämlich a 

BP en a a 

- und A,„(2,t) (man er- . 

kennt das an der Zahl der Doppelpunkte oder rechnet es auch leicht nach). Damit ist aber 
tr = tar. Weiter sind die mit einzelnen Vergleichsfunktionen g»=sinzzundg® =sin2z2 
erhaltenen A®(1, 1) bzw. AP(1,t) mit eingetragen; sie verdeutlichen die durch mehrgliedrige 4 
Ansätze erreichbaren Verbesserungen der Eigenwertschranken. A 
D. Als letztes Beispiel wird die Randwertaufgabe TE 
y"+Al+D)y=0; y()= yll), y’(0) = y'(l) (kurz: yperiod) ..... a) 2 
untersucht. Infolge der Periodizitätsforderung als Randbedingung treten einige Besonderheiten 
auf: Einmal ist A= 0 Eigenwert, außerdem ist zu erwarten, daß die Eigenwerte paarweise eng 
beieinander liegen, da beim ungestörten Problem Y’+AY=0 (Y period.) alle Eigenwerte 
außer A= 0 doppelt sind. Bei der Ermittelung der Eigenwerte nach den bewiesenen Sätzen x 
stört das jedoch überhaupt nicht. &f 2. 
Es sollen wieder die beiden kleinsten aber von Null verschiedenen Eigenwerte bestim 
werden. Da also A,=0 ein nichtinteressierender Eigenwert ist, muß man nach den Überlegun, n 
Gl. (47)*/(48)* von den zur Eigenwerteingrenzung benötigten Vergleichsfunktionen Orthog 
tät zur Eigenfunktion y,= const. verlangen. Dazu werden die Eigenfunktionen Y, der „un- 
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gestörten‘ Aufgabe Y’'’’+AY=0(Y period.) herangezogen; im Ansatz 9 = Y,+ (, wird 
die Konstante C, durch die Orthogonalitätsforderung bestimmt. 

Verwendet werden nur: g{P = cos2rzz und g® = sin2rxr2-+ 1/3rx. Diese Funktionen 
seien gleich nach Gl. (52a)* nochmals iteriert, aus 


ga"—=—(14 eg) go  (P@ period. und | y, im Grundgebiet 0...1) . .. (43) 
folgt ‘ | 
al - sin2n& 10 2 
it) cos27% > © 75 a 
2 1 :  cos2rne 2% % In 168% 
Blut msin2nr+ = elıt2)r gr rn .4051° 


Damit können neben der an 7, anschließenden Determinate D, (t, A) auch Dr—V(t, A) n= 2,3 
(nach Gl]. (32)) und zwar letztere hier besonders bequem erhalten werden. Die Zahlenwerte für 
die Schwarzschen Quotienten von Gl. (4c) sind 


ar 2m EEE, a? =2n2 
20:75 as — — 0,039789 a 01733113 
03, — 0,283751 m a3 = — 0,029364/° a; = 0,269909/r? 
a4 = 0,107770/* a —= — 0,016348 /mt 2% — 0.099697 4: 


Nach Gl. (32) sind die beiden Determinanten DW(t, A) und D@(t, A) aufzustellen, beide können 
benutzt werden, da die Aufgabe ohne den Eigenwert A, = 0 positiv-definit ist. Mit 4 >— 
ergeben sich daraus als Wurzeln obere, mit t, = 8? untere Schranken für },, Ag (ft, = 8x? wurde 
mit dem bekannten Vergleichssatz®) aus den Nebenproblemen y’+ y=0, y’+24M.y=0 
(y period.) erschlossen: 2n? <A,/Ay, <An# Sn? <A,4A, <167?). Durch Lösung einer ein- 
fachen zweireihigen Determinante erhält man die folgenden Resultate, deren Brauchbarkeit . 
offensichtlich ist: 
aus DD(t,A)=0: 2,5479? > A, 2 2,5436? - A, = 2,5467? + 0,089), 
2,8416? > 4, > 2,8351? _A,= 2,8387? + 0,1 /,. 
aus DO (t,A4)=0: 2,5461? > A, > 2,5453 n2 A, = 2,5457? + 0,029), 
2,8387 2 > 1, > 2,8371? A 2,8379? + 0,03 %9. 
I. Die Minimum-Maximum-Eigenschaften der Eigenwerte und eine weitere 
untere Eigenwertschranke 
7. Bemerkungen zur Darstellung des Eigenwertspektrums eines 
Kernes mit höchstens endlich vielen positiven Eigenwerten 
durchein Min-Max-Problem. 


Besitzt ein Kern Eigenwerte beiderlei Vorzeichens, so werden im allgemeinen zur Charakteri- 
sierung der positiven Eigenwerte Max-Min-Eigenschaften, für die negativen Eigenwerte Min-Max- 
Eigenschaften gewisser quadratischer Formen ausgenutzt. Von diesem Wege wurde zwecks 
numerischer Eingrenzung der Eigenwerte bereits in den vorangegangenen Teilen dieser Arbeit 

2 abgewichen. Bei nur endlich vielen positiven Eigenwerten (ein Kern mit endlich vielen negativen 
ordnet sich diesem Fall nach Vorzeichenumkehr unter) lassen sich nun sogar alle Eigenwerte 
fortlaufend aus einer einzigen Min-Max-Aufgabe gewinnen. 

Es sei einmal A; der größte Eigenwert des Kernes X (z, s) (mit nur endlich vielen positiven 
Eigenwerten), dann darf in der Eigenwertdarstellung nach Gl. (20)* der Parameterwert 1; — & 
gewählt werden, und man findet mit Gl. (28)* die bekannte Beziehung”) 


x ? K(z, s) y(x) y(s) dx ds 
Min Max —= Min Max JJKta 9) ya) yi)deds —= Au; 
v a v vita SUR 5) yis) ds? dan U 


- (konkurrenzfähig ist wieder die Gesamtheit der s-Funktionen) und hiermit sind offensichtlich 
schon alle Eigenwerte gewonnen. 
Anschließend soll gezeigt werden, wie auch dieses Min-Max-Problem z. B. zur Berechnung 
unterer Schranken wenigstens bei positiv-definiten Kernen P(x, s) herangezogen werden kann. 
Bi: Das wird durch einen Hilfskern ermöglicht, der zunächst diskutiert wird. 


Teen, (44), 


Fir > * ” . . 
Br. : 8 Die Konstruktion eines Hilfskernes. 


Wie in-den vorangegangenen Teilen der Arbeit, insbesondere in Gl. (12)*, wird wieder ein 
arameter f£, eingeführt: 


a 7. 
» Man) 


E rad). 


_ 2) Sofern der Kern endlich viele negative Eigenwerte hat, setzt man an Stelle von Gl. (44) meist eine 
Max-Min-Darstellung. 


Bar 


Be 
Se 14 
Sr Um Verwechselungen mit dem oben bei Gl. (44) durch den größten Ei igen bestimmten und 
“r schließlich nach > ® vergrößerten t; zu vermeiden, ist hier ein neuer Parameterindex r gewählt. 
= Die Eigenwerte A, gehören jetzt zum positiv-definiten Kern Pi%,'s). 2 | 
Be: Damit wird der Hilfskern H(x, s) gebildet: Ei 
tt, — A 2) | "ae 
Be: Ha, )= 4 Pia )-Pa 32: van) =- IT we a 
be PR (7,8 Pk (HE s) dE ist ein iterier Kern.) ; 
” die angeschriebene Äquivalenz !5)* zeigt, hat dieser Kern H(x, s) ie Eigenwerle ; 
Er; \,= — %,/(t, — A,), die aber in dieser Reihe noch nicht alle der Größe nach geordnet sein müssen. 
Bi: ge t außerdem, daß nur endlich viele Eigenwerte positiv sind — solange nämlich 
u BE have rask 
B Um einen Überblick über den Zusammenhang der Eigenwerte A, von H(z, s) mit Aeri 
Re“ ki ursprünglichen A, von P(&,s) zu bekommen, wird die 
N: | A) Funktion | x ua 
Bi | Pr A FO) = Rs NN: 28.» > 
KR | insbesondere im Gebiet A>0 betrachtet. Sie wächst vn 
| (0) = 0 bis zur ungeraden Polstelle A= {, monotonmit dm 
: | Argument. Für größere A-Werte verläuft sie ganz im Nega- 
SR, tiven mit einem Maximum bei A= 2t,; asymptotisch gt 
1 fa>x)——(A+1t,). In Bild7 ist dieser Verlauf dar- | 
IR i gestellt. 
| Wegen des monotonen Wachstums der Funktion f(2) drh 
i mit Aim Gebiet O<A< t,sind mit den Eigenwerten A, we- ar 
nigstens auch die positiven A, in der nach Gl. (2)* vorge- 
Bild 7. Verlauf der Hilfsfunktion nach 61. (47) schriebenen, der Größe De Folge geordnet use Br 
numeriert, wenn bei TE RT ar 
Ay = Arlle — Ay) (> 0) 2 Rn Te . (48) ee 
gesetzt wird. Die der Größe nach geordneten und numerierten negativen Eigenwerte de oe<0 Sk 
lassen sich nicht mehr so einfach durch die A, ausdrücken, da dann f(A) ER EN keine ein- 0% 
deutige Lösung hat. Jedenfalls ist aber sicher er 
Ace AR (vgl. Bild HER LAT? . (49) H 
Eh 


Die Eigenwerte A, sind nach Gl. (44) und (46) durch eine Min Mas Adtgahe bestimmt, NY 
aus Gl. (48) sieht man, daß hier A, der größte Eigenwert des Hilfskernes H (z, s) sein muß: : 


Min Max le 9 yo) u@deds_ 7, RN Ser 
evt [IH (es) y(s) ds) da a ER: ER Et, BR 


3. della 
9. Die Bestimmung unterer Schradken für die Eigenwerte Fe 


Zufolge der Herleitung bei Gl. (44) ordnet sich das Min- -Max-Problem nach Gl. (50) voll- 
ständig der Extremaldarstellung von Gl. (20)* unter, wenn dorb ;— gesetzt wird. Wie bei 
Gl. (1) und folgende kann die Aufgabe Gl. (50) für die Funktionenklasse Apr gelöst werden — 


damit sind dann schon die gesuchten Schranken für 7, gefunden. Die Lösungen der Min-Max- : | 
' Aufgabe Gl. (50) im Funktionenbereich {p}; seien A,(k, t,). Ihre Bestimmungsgleichung ergibt 
die Determinante D;(t, A)=0 nach Gl. (3), wenn darin entsprechend der Herleitung. von a 
Gl. (44)/(50) zunächst der Grenzübergang t;— » ausgeführt: Dy;(t, A) D;(® N und. K 
dann der Hilfskern aus Gl. (46) eingetragen wird. Das Ergebnis ist DUEN # 
m—Aai=0 . vo=12. ao 


Dabei sind die überstrichenen .a?®, 'a’® noch mit dem Hilfskern H(x, s) zu] bile 
auf Schwarzsche Konstanten, die mit dem Kern Pr (2, s) der Norge 
Integralgleichung bestimmt werden, führt auf 


@e= [[H (x, s) 9,(a) 9. (s sel‘ P%) (m, O—Pins s)] 
1.02 — are 

’ | we /[[H( (2, 5) 9 (s)ds- DR (2, 8 DL 2 

FE Wieder ist nach Gl. (de) / a Kg 

= mm an [Plz s) Qy( on deds 
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Die Ordnung der Wurzeln A,(k, t,) folgt nach Gl. (5) wenn dort 1;— oo strebt und dann an 
Stelle A;4z(k, 0) entsprechend den gewählten Bezeichnungen A,ı;(k, t,) gesetzt wird: 


| N .. n. a). 
Aus den Gl. (16b) und (17) folgt ebenso 
A,y(k,t) <A (ek 1,6) <A; und im Arssi(k,t,) = Any (54a, b). 
"up 


37-0, 1,2...) 


Die hiermit berechneten Größen liefern aber sofort auch die gesuchten unteren Schranken 
für die Eigenwerte des positiv-definiten Kernes P (x, s). 


10.7 Untere Schranken für die Eigenwerte von P(z,s). 


Sobald > A; und solche Koordinatenfunktionen 9, bzw. in solcher Anzahl benutzt 
werden, daß einige der aus Gl. (51)/(53) erhaltenen Hilfsgrößen A,+;(k, t,) positiv sind 


Beier tkm <lr— 1). 22.8... (55) 
(da es nur r positive Eigenwerte A; gibt, muß nach Gl. (54a) m < (r—1) sein), können damit 
(positive)’ untere Schranken für die Eigenwerte ),1; 3=0, —l,...., —m angegeben werden. 
Zu diesem Zwecke sei gesetzt: | 
ar +; (k, t (k Ban TEAR HE FE 
Arzi(k, t,) = en ) > V:also 0 <l, 1, ) = Ara; Arzilkitr) {| 1 EN ‚| (56a, b) 
br —br+jlk, tr) 2 73 I la 
Arrj(k, t,) 
3=0, —1,..,—m 


GI. (54a) erhält die Gestalt | 
brilk,t,) - buk+lhm) __ An 
, — bı+j(k, 4) — b»+;(k +1,45)" u — Ari 
Da nun nach GI. (56) sicher 0 < 1,ı;(k, t,) < t, ist und in diesem Gebiet die Funktion f(A) 
nach Gl. (47) monoton mit 4 wächst, folgt daraus 
bl, t,) <urj(k +1 ,1,)< Arry ar, IF 0 2.2, Mm. u (08). 


Damit ist schon die gesuchte von unten gegen Eigenwerte des Kernes P(z, s) konvergierende 
Näherungsfolge gewonnen. 


0< Bl EM OT 


11. Bemerkungen zur HandhabungundAnwendungder 
bewiesenen Schranken. 


+ Um die unteren Schranken nach Gl. (58) aufzustellen, muß zuerst mit geeigneten Ko- 
 ordinatenfunktionen 9, die Determiinante nach Gl. (51)/(52) zusammengestellt werden; dabei 

sind die Koordinatenfunktionen zweimal zu iterieren, damit die erforderlichen Konstanten 
A — [gm ve dz O<m<sn=0,1,... durch einfache Quadratur gebildet werden kön- 


nen. Über Aie Durchführung dieser ‚Rechnung und die Bestimmung des Parameters t, 
informiert man sich etwa im Abschn. 5 und bei den Beispielen in Abschn. 6. Die Wurzeln 


Er, Ay(k, {,) der charakteristischen Determinante Gl. (51) werden nach Gl. (53) numeriert und 
FR dienen ME — soweit sie positiv sind — in der zweiten Formel von Gl. (56) zur Berechnung 
= der gesuchten unteren Eigenwertschranken I,+;(k, t,). 


‚Die notwendige zweimalige Iteration der Koordinatenfunktionen erschwert die Anwendung 
der Schranke wesentlich; ein Vorteil des Verfahrens ergibt sich bei der Berechnung von Eigen- 
en ınktionsnäherungen. [In diesen Arbeiten wurde allerdings auf diese Aufgabe nicht weiter ein- 
gegangen; es sei nur bemerkt, daß die zu den Werten A; 1;(k, t;) bzw. I,+;(k, t,) gehörigen Extre- 
malfunktionen auch als Näherungen für die zugehörigen Eigenfunktionen angesehen werden; 
; Bestimmung erfordert die Berechnung der Konstanten p, von Gl. (1) mit einem linearen 
nogenen Gleichungssystem, dessen Koeffizientendeterminante die Gl. (3) bzw. (51) darstellt.] 
I Se besteht darin, daß infolge der gleichzeitigen Brneepe verschieden ner Iterations- 
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‚bare Schranken. 
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darstellen. Ist y aus {yt; die zugehörige Extremalfunktion, und benutzt man hier die mit Y 
und dem Kern P(z, s) entsprechend Gl. (26)*, (27)* gebildeten Schwarzschen Konstanten 
und Quotienten, so findet man den Zusammenhang 


m % 


A, (kt) {bg P@ (x, s) — P(x, s)| y(x) y(s) dx ds ER L a,(4) — a, (y) Bi 
Be ff PO iz, s) — Pi S)] Y (3) ds}?d« L a,(Y) a t,Ag(Yy)-+az(Y) (59), 
us(y) way) Ir —Mely)] 7, %) T3(y 6) 


—ayE +) la) m Taly, %) 
Wird dem Gl. (56) gegenübergestellt, dann erhält man die Abschätzung 7,(9, %,) .... urslk, tr)... 
T (9 t) 5 Ir (k, t,) wird durch eine Mittelbildung aus zwei Temple schen Quotienten erhalten. 


1 BerispierR 
Es wird das bereits im Teil I, Abschn. 6 unter D. betrachtete Randwertproblem 
y’+til+an)y=0;  4y0)=yd), y()=yll)..... . (60) 
benutzt. Da die Aufgabe positiv-definit ist (vom Eigenwert A,= 0 wird abgesehen), so kann 
die Eigenwertschranke nach Gl. (58) versucht werden. Sie wird mit dem gleichen zweigliedrigen 
Ansatze wie in Abschn. 6, D. durchgeführt, die Koordinatenfunktionen sind zu %, orthogonal, 
so daß der Eigenwert A, = O in den Formeln überhaupt nicht stört (vgl. in 6. D.). Die zur Auf- 
stellung der Determinante Gl. (51)/(52) erforderlichen Konstanten «?® sind schon im Abschn. 6D. 
ausgerechnet, so daß diese sofort hingeschrieben werden kann: 
Ay = 2° —1t, 0,75 — A[0,75—2t,.0,283751/r2? + 170,10 7770/r2] 
—=0 mit A,,= t,0,039789 — A.[—0,039789 + 2 1,0,0293644 /? — 17 0,016348/r°] (61). 
Ay,= 27° -1t,0,733113- 4[0,733113— 2 t, 0,269909/r2 + 120, 099697 [1 
Mit t%=8n2 (vgl. Abschn. 6, D.) gibt das zunächst A, (2, 1,) = 1,18693 7; 
A, (2, t,) = 1,5579 2 und dann nach Gl. (56) b 
1, >21, = 2,54463 7? 
is 2b, = 2,83631 7", 
Ein Vergleich mit den Ergebnissen in Abschn. 6, D. zeigt die Güte dieser unteren Schranken 
und außerdem die bei Gl. (59) behauptete Zwischenlage T,...4,+j... Tz. 


Aıı An 
Ars Ay 


Zusammenfassung 

Das Ziel dieser und der vorhergegangenen Arbeit ist die Eingrenzung der Eigenwerte be- 
liebiger (nicht notwendig positiv-definiter) linearer und selbstadjungierter Eigenwertprobleme, 
insbesondere von Integralgleichungen; es wird mit Hilfe einer neuen Extremwertdarstellung 
der Eigenwerte erreicht. Unter Ausnutzung leicht erhältlicher grober Eigenwertschranken werden 
Näherungsfolgen konstruiert, die von unten und von oben gegen einen vorgegebenen Eigenwert 
konvergieren. Die erforderliche Rechenarbeit hängt dabei wesentlich von der vorgeschriebenen 
Eingrenzungsgenauigkeit ab, bleibt aber bei passend geführter Rechnung wohl immer in er- 
träglichen Grenzen. 

In vielen Fällen, sofern z. B. „‚gestörte‘‘ Probleme vorliegen oder in denen die Lösungen 
„ähnlicher“ Aufgaben bekannt sind, liefert bereits ein einfachster Einschließungssatz brauch- 


Die dabei benutzten Extremwerteigenschaften der Eigenwerte zeigen außerdem einen 
engen Zusammenhang verschiedener bisher nur nebeneinander benutzter Extremal-Prinzipe und. 
der darauf gegründeten Verfahren zur Berechnung oberer Schranken für die Eigenwerte. Benutzt 
man den T-Quotienten 7, (y, t), so ergibt sich für das Grundprinzip des (Ritz)- Galerkin- 
schen F Verfahrens Max Min #2(y) = MaxMin T,(y, ©) und für das von Grammel_ 
Max Min 4,(y) = Max Min T,(y, 0), sie unterscheiden sich also nur durch die Wahl des Para- 
meters £;. Gleichzeitig erhält man volle Klarheit über die Tragweite dieser älteren Methoden. 

Nebenher bietet sich hier ein interessanter und verallgemeinerungsfähiger Zugang zu ö 
Kr yl off schen (=Weinsteinschen) Einschließungssatz, der damit auch im Rahmen vo: 
Maximum-Minimum-Prinzipien verstanden werden kann — außerdem ist gezeigt, daß er 
gleichen Gültigkeitsbereich hat, wie z.B. das Verfahren nach Galerkin. Baker: 

. In einem abschließenden Teil wird vorgeführt, wie sich noch auf anderen Wegen [hier über 
einen einfachsten Hilfskern vom „Polynomtyp“ a,K (2, s)+ a, KM (x, s)+...+@,K® (x, s)] 
brauchbare untere Eigenwertschranken gewinnen lassen. i ES 
Eingegangen am 20. März 1949. Bi ? + = BR 
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| Die Geschwindigkeitsverteilung 
bei lokalen Überschallgebieten an flachen Profilen 


Von K. Oswatitsch in Emmendingen. 


Die Arbeit gibt eine erste kurze Darstellung einer Theorie zur Berechnung lokaler Überschallgebiete an 
flachen symmetrischen Profilen. Nach Vereinfachung der Ausgangsgleichungen wird eine Integralgleichung 
für die Geschwindigkeütsverleilung in der Strömungsebene aufgestellt. Diese Integralgleichung wird durch 
einen Nüherungsansatz auf eine Integralgleichung für die Gesshwindigkeitsverteilung am Profil reduziert. Die 
Lösung erfolgt durch Parameleransälze, wobei die Gleichung in einer Anzahl von Punkten exakt erfüllt 
wird. Es ergeben sich die experimentell bekannten Geschwindigkeitsverteilungen mit Verdichtungsstoß. 

The paper gives a first condensed description of a theory of local supersonic aeras on thin symmetrical 
airfoils. Having simplified the basic equations, the author formulates an integral equation for the distri- 
bulion of velocities in Ihe whole field of flow. This integral equation is reduced by Ihe aid of a formula of 
approzimation to un integral equation for Ihe distribution of velocities on the airfoil, The solution is given 
by supposed parasmeler-formulas which will exactly meet the equalion for a certain number of points. The 
results show the distribution of velocities with shock-wave well-known by experimental researches. 

L’auleur communique un premier expose suceinct d’une theorie pour le calcul de domaines locauz super- 
soniques sur profils plais et symetriques. Apres la simplification des formules de depart, il &lablit une 
Equation integıale pour la repartition de vitesses dans le plan de l’Ecoulement. Cette &quation integrale est alors 
reduite & une Equation intögrale pour la röpartition de vitesses sur le profil a l’aide d’une formule d’approwi- 
malion. La solution s’oblient par formules de parameitres ou l’Equation est exactement definie pour un cer- 
tain nombre de points. Le resultat met en Evidence les r&parlitions de vilesses wvec onde de choc qui sont 
connues par les recherches experimentales. 

B 3aToü pabore naerca TNepBoe KpaTKkoe H3NOKEHHEe TEOPHH PAcyeTa MECTHEIX CBepX- 
3BYKOBEIX OOAACTeü y TOHKUX CHMMETPH4YHBIX Ipobmneii. Jlocıe ympomenna UCXOAHEIX 
ypaBbHeHunfä aBTOP BbIBONuHT HHTErPalbHoe ypaBHeHne IA PpacımpezeneHnn CKOPoCcTeh B 
INIOCKOCTH HOTOKA. OTO HHTETPAAIBHOe YpAaBHeHHe CBOAHTCH IPH HOMOINH IPHÖJHSKEHHBIX 
$opMmyA K HHTErPAAIbBHOMy YPABHeHUM IA pacıpenereunsg CKOpocTeä HA -IIOBEPXHOCTH 
npoßnaa. Pemenne HAXoNuUTcH IyTeMm AONyINeHnA TAKUX TAPaMeTPOB, KOTOPpble TOYHO 
YAOBIETBOPAWT YPABHeHHWO B HEKOTOPoM uucae Toyer. llonyyaerca USBeCTHOe M3 OIIEITOB 
pacıpeperlenne CKOPocTeä C THAPABIHYeCKHM YAapoM. 


1. Grundgleichungen 


Es ist bekannt, daß die im folgenden behandelte stationäre, drehungsfreie ebene Um- 
strömung flacher Körper bei kleinen Machzahlen M (Verhältnis der Strömungsgeschwindigkeit 
w zur Schallgeschwindigkeit c) durch die Laplacesche Gleichung, bei Überschallgeschwindigkeit 
durch die Wellengleichung beschrieben werden kann. Die mathematische Aufgabe wechselt 
danach bei Steigerung der Machzahl im Anströmgebiet M„ von einem elliptischen zu einem 
hyperbolischen Problem. Das in Schallnähe auftretende gemischt-elliptisch-hyperbolische 
Umströmungsproblem hatte bisher allen Lösungsversuchen getrotzt. Es konnten nur exakt und 
näherungsweise Strömungen angegeben werden, die im Unendlichen keine Parallelströmung, 
ergaben. 
Im folgenden soll eine Lösung der Aufgabe wiedergegeben werden, die weniger den reinen 
Mathematiker wohl aber den angewandten Mathematiker und Aerodynamiker befriedigen dürfte. 
_ — Die ausführliche Arbeit und der Vergleich mit den Versuchen folgt in den Physica Acta Austriaca. 
E, Mit w, v als Geschwindigkeitskomponenten und %, % als kartesischen Koordinaten lautet 
die das Problem beherrschende gasdynamische Gleichung (1) und Gleichung für die 
' Drehungsfreiheit (2) 
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Für flache Profile ist mit Vernachlässigung der quadratischen Glieder: 
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mit die gasdynamische Gleichung die Form (4) erhält: 

\ du dv va 

Br — M?) — + — — _ BR} 22 0 SER LEN 
Br \ | re u? U =? (4) 
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Beim vorliegenden Problem kann bekanntlich die Schallgeschwindigkeit und die Machzahl 
Funktion des Geschwindigkeitsbetrages allein angesetzt werden. Es genügt hier linear zu 
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FUN 
nähern, wobei verlangt werden soll, daß die Machzahl in der Anströmung (M = M.) und bei 
Schallgeschwindigkeit (M = 1) richtig wiedergegeben werden soll. Damit ist: b 


— M%, 2 

en -ı (2a N .. . ee (oe 
6 

mit c* als kritischer Schallgeschwindigkeit (bei M = 1). Das erste Glied dieser Entwicklung 


nach —- — 1 ist dabei das Quadrat des bekannten Prandtlfaktors: 


Be=yı — MEN m... (0): 
Besonders einfach werden die Gleichungen für folgende Abkürzungen: u 
s=%; y-= "ß; 

- u ® R 
. BER, SE Us » 9); 
Be Me u (a u & ® | (N. ’ 
% 0% cr — Us 1 on * Ble*— us) i 

ER A B\— —1 } 
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Es ist zu beachten, daß U und V in Schallnähe von der Größenordnung 1 werden. (Esist | 
für M=1; U =1.) Nur weit ab von der Schallgeschwindigkeit sind U und V klein. 


Mit den Symbolen (7) lauten nun Gl. (4) und (2): d: 
ou 0V ou e ) oV Br“ 
FRRER UA rn ar ee N (8); 
20. av 
Sa U a (9. _ 
Y 0% 


Das erste Glied der rechten Gleichungsseite von (8) entspricht der Veränderlichkeit des Faktors 
(1— M?), das zweite Glied den beiden mittleren Summanden der gasdynamischen Gleichung 
Gl. (1). Dieses wird beliebig klein bei Näherung der Anströmung an die Schallgeschwindigkeit 
(Us c*). Die rechte Gleichungsseite von (8) verschwindet weit ab von der Schallgeschwindig- 
keit. Die Resultate entsprechen dann genau der Prandtlschen Regel. & "u 
Das Profil sei durch die Vorgabe von v auf der x-Achse (y = 0; v = v,(x)) gegeben. Dann 
lauten die Randbedingungen im Unendlichen und auf y=0: N 


Ye +y2>o: U>0; v‚>0 
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OT ied 
u 
at: ö 
Bei asymmetrischen Profilen oder Anstellung des Profils wäre V,(x) auf beiden Seiten der z-Achse 
vorzugeben). i SE 
h * * | a re 
br - Der Ausdruck — — 1 ist proportional zu ß? (p= (+1) [2 IR 1) ++. ) Die Rand- 5 
bedingung (10) entspricht also nicht der Vorgabe eines Profils sondern jener einer S char IS 
\ affin verzerrter Profile, wobei die größeren v,(x)-Werte, also die größeren Dickenverhältnisse, 
bei gleichen V,(x) den größeren ß-Werten, also den kleineren Machzahlen der Anströmung ent- 
sprechen. Es sollen in dieser Arbeit nur so flache Profile behandelt werden, bei denen u» ge- 
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y=0: V=V,y(a) = 
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mit den Randbedingungen (10. ' ih 
. . Piese Einschränkung ist für die Durchführbarkeit der folg. 
lich, doch soll auf diese Weise das Wesentliche der Aufgabe 
eifekt getrennt werden. Profile, welche den Gleichungen (11) un 
(10) bei gleichem V,(x) genügen, gehorchen also einem Äh 
welchem dem dickeren Profil die niederere Machzahl der Anströ 
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kung der Geschwindigkeilen am Profil entspricht. Dieses Gesetz wurde von Guderley, 
v. Kärmän und von mir etwa gleichzeitig entdeckt, ohne bisher ausführlich im Druck er- 
schienen zu sein. 

Wird die Stromdichte ou so entwickelt, daß sie im Ansirömgebiet (po u = 0% us) den 
richtigen Wert und — entsprechend zur Prandtlschen Regel — auch die richtige Tangente an 
die exakte Stromdichtekurve hat und und ferner wie diese bei M = 1 ein Maximum aufweist, 
so lautet’die Reihe: 


nee 


Es kann also mit Gl. (7) geschrieben werden: 


BONS, | 
en U = {133 
A) 

U 


Darnach stellt der Ausdruck U — „m wesentlichen die Stromdichte dar, was sich auch 


Zu 


darin zeigt, daß die die Kontinuilätsbedingung darstellende gasdynamische Gleichung (11) auch 


geschrieben werden kann: 
0) Br oV 
En a 
( 2 °y N (14) 
Ein senkrechter Verdichtungsstoß ist ein Sprung von Überschallgeschwindigkeit (U > 1) auf» 
einen Wert Ü < 1 bei gleicher Stromdichte. Es folgt aus 
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Das ist die bekannte Gleichung für schallnahe senkrechte Stöße — nur solche kommen innerhalb 
dieser Theorie flacher Profile in F sagt, daß das arithmetische Mittel der Ge- 
schwindigkeit vor und hinter dem Stoß gleich der Schallgeschwindigkeit ist. 

Man überzeugt sich leicht, daß im Gleichungssystem (11), (12) die charakteristischen 
Eigenschaften eines gemischt elliptisch-hyperbolischen Problems voll erhalten sind. Das System 
ist im übrigen keinen unkontrollierbaren Vernachlässigungen entsprungen. Es gibt einschließ- 
lich der Beziehungen für den Verdichtungsstoß Gl. (15) und (16) die genaue Beschreibung 
der Vorgänge für die flachsten in der Schar der — dem Ähnlichkeitsgesetz entsprechend 
— affin verzerrten Profile. 


2. Eine Integralgleichung 


Bi Es kann nun im System (11), (12) in üblicher Weise ein Geschwindigkeitspotential U = ®,, 

v=®, Seen werden. In der Gleichung für ®: 

“aa ou 

er RN BP a ee N 
RR}: > 


‚stellt dierechte Seite die Quellverteilung in einer Poissonschen Gleichung dar. Die PrandtlIsche 
| > 0 3 & x 
_ Regel, welche der Laplaceschen Gleichung (v 2 = 0) und einer linearen Näherung der 


Stromdichte entsprechen würde, setzt zu große Stromdichten ou voraus. Die wirkliche Strömung 
mit ihrer kleineren Stromdichte entspricht einer Prandtlschen Strömung mit einer Quellver- 
teilung im ganzen Feld. Dabei entspricht die Zunahme der Geschwindigkeit einer Zunahme 
_ der Verdrängung und mithin einer Quelle. Die Quellwirkung verschwindet mit der Entfernung 
vom Körper quadratisch, und das Abklingen weit draußen gehorcht denselben Gesetzen wie 
jenes bei der Prandtlschen Regel oder auch bei der inkompressiblen Strömung. 

Die Lösung der Poissonschen Gleichung bei gegebener Quellverteilung ist bekannt (Cou- 

rant- Hilbert, Methoden der Mathematischen Physik II, S. 230). Zu den Quellen in Gl. (17) 
kommt noch eine Quellbedingung auf der x-Achse, welche die gewünschte v-Verteilung vorgibt. 


gr ou 
einem symmetrischen nicht angestellten Profil ist die Quellbelegung des Feldes (v ) eben- 


eymetrisch zur &-Achse, ändert also die v-Verteilung des Profiles nicht, so daß also die 
{ Bene auf y=0in derselben wese erzwungen wird wie beider Prandtlschen Regel 
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enthalten ist. Die Aufgabe ist auf die Lösung einer eigentümlichen Integralgleichung zurück- Pr 


KR x) = U, gilt: 


das Doppelintegral zeigt dort einen stetigen Verlauf. 
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oder der inkompressiblen Stroee. Man erhält auf diese Weise für das DR ®, eines s sym NER 


metrischen nicht angestellten (nur dieser Fall soll hier behandelt” werden) Pro- 
AR 


files: B 4 “ 
Ua, y) = U,( - 4 [Taen: BEN. 2 ası as) Ei 
%,Y SS ee N . A 
Lee een 2. 
Hierin ist U,(z, y) die durch die Prandtlsche Regel errechnete Geschwindigkeitsverteilung ‘ 
U,(&, 9) = :@ ee |; “ . 
(%, Y) are zer; ve: E-a® ee e B ß ( Ya 
Auf y=0Oist in Gl. (19) der Cauchysche Haupt.wert zu nehmen. Dasselbe gilt auf 
der Geraden 7 = y für das Doppelintegral i in Gl. (18). Es ist definiert durch: | 
1 ‘ + +% au & \ , e u; . 
2: Be TE iR . SS 
= a X Een —y) FT ee 
Bes .: De 
—= — lim : N . U — dnd& + 20),. Kar 
Insnol, See 08 ee e en 
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+ [ rn dnk 
a, ee 
j 
was bei den folgenden Rechnungen ganz besonders zu beachten ist. 3 
In der Nähe des Dickenmaximums des Körpers erweist sich das Doppelintegral als negativ, Ak 
so daß die gesuchte Geschwindigkeitsverteilung U (x, y) dort dargestellt wird als die Prandtl- 
| 


Bi ı ce, 
verteilung U, vermehrt um einen Zusatz, in welchem allerdings wieder die gesuchte Verteilung RR 


geführt. Am Doppelintegralläßt sich nun noch eine partielle Integration ausführen, womit unter 
Beachtung von Gl. (20) SI ee gewonnen wird: i 


2 EN { 
Vay= Ua + Ley Ze ae m den. 


wieder mit einer Gleichung (20) entsprechenden Hauptwertsdefinilion für das Doppelintegral. iR 
Die Integralgleichung (21) stellt die fundamentale Gleichung dieser Arbeit dar. Das Doppel- 
integral hat hierin keinen Sprung an Sprungstellen von U und ist im übrigen in der Nähe des de n 


er Y) 


u 

2 
keitsverteilung kann nach Gleichung (21) dargestellt werden &E U, vermehrt um einen Ausdruck, N j. 

ß 

Kar 


der‘ 2 aus U?/2 und dem Doppelintegral zusammensetzt. Ebenso kann aber die Stromdichte R £ 
AZ aus U, vermindert um das Doppelintegral berechnet werden. ‚Die erste Darstellung | 


hat den Vorteil, daß die Geschwindigkeit sich aus einem verhältnismäßig kleinen Zusatz z zu 17 B N. 


ergibt, während sich der Unterschied von U, und U a an den wichtigen Stellen als ru dd % 
[2 


doppelt so groß erweist. Dieser Unterschied muß in Schallnähe auch sehr genau berec 1 
werden, weil sich sonst die en. aus der San nur el Breenan berec 


Dickenmaximums des Körpers positiv, aber kleiner als . Die gesuchte Geschwindig- % 


den die Berechnung der Stromdichte vor der Berechnung der a an Glkei 
Es werden auf diese Weise die Sprünge in den Stößen in richtiger Gr j 


Es interessiert vor allem die Genre P rofi 
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- Doch ist auch zur Berechnung von U, die Kenntnis der Be indielieverteie: imganzen 
Feld erforderlich. 

Es ist kein Fehler, wenn U mit Gl. (21) und (22) auch im lokalen Überschallgebiet in einer 

dem elliptischen Differentialgleichungstyp entsprechenden Weise dargestellt wird. Jedoch ver- 


bietet sich bei lokalen Überschallgebieten eine Iteration in der Weise, daß ausgehend von der. 


Prandtl-Lösung U, die Werte in der rechten Seite von Gl.(21) eingesetzt werden, um eine 
nächste Näherung für u zu erhalten. Dieser Vorgang entspräche der Methode vonRayleigh- 
Jantzen, welche bekanntlich nicht zu den gesuchten Überschallfeldern führt. 

Eine exakte Lösung von Gl. (21), welche der exakten Berechnung der Strömungen an 
flachen Profilen entspräche, dürfte sehr schwer sein. Eine Iteration mit Gl]. (21) ausgehend 
von einer guten Näherung des Resultates hingegen wäre denkbar. 


3. Näherungslösung der Integralgleichung 


Die Geschwindigkeitsverteilung U auf einer Geraden x = konst. hat stets einen charakte- 
ristischen Verlauf. Sie fällt vom Werte U = U, mit einer durch die Wandkrümmung gegebenen 


| 2 7 e : 
Anfangstangente ((&) BR ) ‚ wird: bei größeren y-Werten flacher und verschwindet weit 


o da 
draußen mit 1/y?. Während in Wandnähe die Wandkrümmung maßgebend ist, ergibt sich für 
den Verlauf weiter draußen eine Bedingung aus der Integration der Kontinuitätsbedingung über 
y. Hier soll (für y > 0) mit dem Ansatz: 


U(z, y) ga U ,(&) 


Yy 2 
gearbeitet werden, in dem U, die gesuchte Geschwindigkeit auf y= 0 und b eine Funktion b(z) 
ist, die so gewählt werden kann, daß die Drehungsfreiheit auf y = 0 oder die Kontinuitätsbe- 
dingung in integraler Form oder beide Bedingungen nur teilweise erfüllt werden. Unter Erhöhung 
des Arbeitsaufwandes könnte die Näherung der Geschwindigkeitsverteilung durch Hinzunahme 
eines zweiten Parameters natürlich noch verbessert werden. 

Der Ansatz Gl. (23) in das Doppelintegral von Gl. (22) eingesetzt, gestattet die Ausführung 
der Integration über 7 und gibt eine Näherungsintegralgleichung für U,(2): 


. (23) 


+ 
22 U, U,(9) En ee 9% 
D,= | be)? b(ö EEE TEA Y 2.) 
mit der etwas komplizierten Funktionsgestalt von E (für 2> 0): 
| 2 4 IE A „2 24 
E(z) aa La de + 24] — [1 — 102? + 52] In z ee 
R — = (1 +22) [25 — 7122 — 22 — | 


_  mittlung von E(z). Die Funktion ist in z 
_ symmetrisch (Bild 1), wächst bei 2=0 
er  logarithmisch über alle Grenzen und ver- 
- schwindet im Unendlichen mit 1/22 ganz 
wie die durch einen Körper in einer Paral- 
— lelströmung verursachten zachrndz 
keitsstörung. 
Bat; Für die Stelle x = E ist wieder dr 
_  Gauchysche Hauptwert zunehmen. Die 
 Singularität von E an dieser Stelle ist 
 logarithmisch, verursacht also keine Schwie- 
 rigkeiten. 
Die Aufgabe ist mithin auf die Lösung 
2 es Näherungsintegralgleichung für U oe &) 


\ für z< 0 bedarf es einer gesonderten Er- 
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Bild 1 


3 ireh die Prandtl- Verteilung U zo( x) in die Rechnung nerhke Bezüglich der Iterationsmög- 

lichkeiten gilt für Gl. (24) dasselbe wie für Gl. (21). 

Be + Zur Berechnung der Geschwindigkeitsverteilung im Gebiete der Übergeschwindigkeiten 

ts 20) werden wegen des starken Abklingens von Z die Untergeschwindigkeitsgebiete (U<0) 
Bene des vorderen und hinteren Staupunktes nur eine untergeordnete Rolle spielen. 

er des Integrals in Gl. (24) daher wichtig sein, vor allem die Geschwindig- 
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keitsverteilung in der Nähe des Geschwindigkeitsmaximums richtig nachzubilden, vor allem 
auch, weil U, quadratisch in den Integranden eingeht. Treten keine Verdichtungsstöße auf, so 
hat U in der Umgebung des Maximums parabolischen Verlauf (Bild 4a). Bei Verdichtungs- 
stößen mit größeren Sprüngen (Bild 4e) genügt eine Näherung der U-Verteilung durch eine 
-Halbparabel, da die niedrigen Geschwindigkeiten hinter dem Stoß im Integral der Gleichung 
(24) nicht ins Gewicht fallen. Nur bei schwachen Stößen ist eine Näherung durch je einen Halb- 
parabelberg auf beiden Stoßseiten erforderlich. Für den Wert des Integrals in Gl. (24) ist ferner 
die Funktion b (x) von Bedeutung, welche ein Maß für die Höhe des Störgebietes darstellt. Nun 
hängt das Abklingen der Störung weiter von den mittleren Profileigenschaften ab. Es ist daher 
naheliegend, mit einem mittleren konstanten b zu rechnen. Die Berechtigung dazu wird sich 
später nur bestätigen. 

Ein Parabelberg oder Halbparabelberg hat dann drei wesentliche Parameter. Die größte 
Höhe gegeben durch U, „, die Breite x, und seine Lage, etwa gegeben durch die Koordinate 2, 
seines Maximums. Damit ist er gegeben durch: 


U = U N SE BE 2 
0 


welcher Ansatz nur für U, >0, also |e — x,| < & gelten soll. 

Nach Einsetzen von Gl. (26) in Gl. (24) und Durchführung der Integralion läßt sich die 
Integralgleichung bei gegebenen b entsprechend den drei verfügbaren Parametern gerade an 
drei Stellen exakt erfüllen, woraus sich U,„,o, %m und x, bestimmen läßt. Die Näherung (26) 
wird im Folgenden nun nicht als eigentliche Geschwindigkeitsverteilung verwendet sondern 


nur zur Berechnung des Unterschiedes von U, und U,,, oder auch zur Berechnung des Unter- 
U; i 
schiedes von U,— ;’ und U,o. Die erste Methode kann wegen der erforderlichen Kenntnis 


von U5/2 nur dort angewendet werden, wo U, durch Gl. (26) gegeben ist, also vor allem im 
Überschallgebiet. Die zweite Methode liefert überall Werte, wird jedoch für Werte nahe an 
U=1 zweckmäßig nicht angewendet. 


Bei konstantem b treten beim Ausrechnen des Integrals in Gl. (24) nur = und damit auch 


x 
nur der Wert ze auf, Es kann damit gesetzt werden: 


LS, = 5) ae Im 5 ( >) 
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Die Funktionen f, und f, sind für den Parabelberg und den Halbparabelberg in Bild 2 und 3 


eingetragen. Sie können für die Rechnungen an ganz beliebigen Profilen verwendet 


werden, solange man sich mit den Näherungen (23) im Feld bei konstantem mittleren 5 und mil 
Gl. (26) zur Berechnung der oberen Funktionen begnügt. Es fällt dabei besonders beim Halb- 


parabelberg auf, daß der Wert am Maximum kaum vom Werte 23 ‚also auch kaum vom Wert 
b selbst abhängt, was eine gewisse Sorglosigkeit in der Wahl dieses Parameters rechtfertigt. 


4. Anwendung auf das Kreisbogenzweieck 


„Bild 4a bis e gibt eine Anwendung auf ein Profil, welches auch bezüglich des Dicken- 
maximums symmetrisch ist und bei Unterschallgeschwindigkeit daher auch symmetrische U- Ä 
Verteilungen liefert !). Es wurde b = 4 gesetzt, was dem Wert von 2) am Maximum beider 
Prandtlverteilung entspricht. Bei der symmetrischen Lösung (Parabelberg) wurde Überein- 
stimmung im Maximum und in jenem Punkte gefordert, in welchem f,[Gl. (26)] verschwindet, 
wo also U,= U5o ist. Außerdem weiß man hier, daß der Parabelberg symmetrisch liegt. Bei 
der Lösung mit Verdichtungsstoß wurde der Nullpunkt des Halbparabelberges mit dem Null- 


23 


punkt der Prandtlverteilung gleichgesetzt, was kaum einen Fehler ergibt. Übereinstimmun Ri 


: ‘) Es handelt sich um die Randbedingung V, = 2r(1— 22) für 0o<x<I1T, also im wesentlichen 
ein Kreisbogenzweieck. Dabei ist 7 ein „‚reduziertes‘‘ Dickenverhältnis, das entsprechend dem Ä 
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gesetz mit dem richtigen Diekenverhältnis des Profils 7’ in folgender Beziehung steht: DE N 
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wurde wieder im Maximum und im Punkte f, = 0 gefordert. Die Rechnungen. können rasch 
durchgeführt werden, nachdem die Arbeit der Berechnung von f, und F, ein für allemal erledigt 


ist. Somit läßt sich auch der Effekt einer Änderung von b leicht studieren. 
Es gibt bei den Verhältnissen, welche Bild 4d zugrunde liegen, eine symmetrische une Eine: oea 
asymmetrische Druckverteilung. Ob die symmetrische richtig wiedergegeben ist, wird noch 
i 2 


genauer studiert werden müssen, denn bei der Berechnung von U, über U, = ergäbe sich k 


wohl auch eine symmetrische Lösung, welche allerdings vor dem Dickenmaximum einen Ver- 

dünnungsstoß und symmetrisch dazu hinter dem Dickenmaximum einen Verdichtungsstoß auf- 

weist, welche also thermodynamisch unzulässig ist (punktierte Kurve). Praktisch ist diese Frage 

allerdings weniger interessant, da bei stärkerer Überschreitung der Schallgeschwindigkeit die 
asymmetrischen Lösungen interessieren. Bei noch hö- 
heren Anströmungsgeschwindigkeiten liefert die Rech- 
nung überhaupt keine symmetrischen Lösungen mehr 
(Bild de). Interessant und praktisch bedeutungsvoll 

- ist die Beschleunigung der Strömung hinter dem Stoß. 
Bei reiner Unterschallströmung ai 4a und b)ergibt 
sich nichts neues. 
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5. Rückblick und Aushlick IR, 


Die hier entwickelte Theorie wird nach der Aufstellung einer Integralg RR Gl. (21) 
gewonnen, indem qualitativ richtige Lösungssätze, welche verfügbare. eier enthalten, 
die Integralgleichung eingesetzt werden und die exakte Erfüllung der Integrs Igleichung in in 
der Parameterzahl entsprechenden Anzahl von Punkten gefordert 
Theorie mit den Versuchen steht allerdings noch aus. Weil sich eı 
Geschwindigkeitsverteilungstypus in einem sehr kleinen Bereich £ 


Sicher wird diese Theorie noch verbesserungsfähigsein. Vor all 
entsprechende Lösungen bei Anstellung des Profils und bei Übers 
Überschallgebieten zu erhalten. © 

Eingegangen am 28. 4. 49, 
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Über die äußere und innere Reibungsschicht bei 
Störungen laminarer Strömungen 
Von H. Holstein in Göttingen 


Das Formelsystem, das zur Durchführung von Stabilitätsrechnungen an laminaren Reibungsschichlen 
nach der Tollmien schen Methode erforderlich ist, wird — z. T. auf neuartige Weise — hergeleitet und 
übersichtlich zusammengestellt. Für die dabei auftretenden universellen Funktionen werden genauere Ta- 
bellen angegeben. Insbesondere zum Fall der angefachten und gedämpften Schwingungen werden einige 
neue Resultate mitgeteilt. Zu Ergebnissen aus dem neuesten Schrifttum zur Stabilitätstheorie, insbesondere 
von C. C. Lin und D. Meksyn, wird kritisch Stellung genommen. 


The system of formulae necessary for stability calculations on laminar boundary layers given by Toll- 
mien, is deduced in a partially new way and ordered systematically. More accurate tables are given of 
Ihe universal functions occurring there. Some new results are published, especially for the case of amplified 
and damped oscillations. Some conclusions given in recent publications on boundary layer stability, such 
as 0. C. Lin’sand D. Meksyn’s papers, are critizised. 

Pour l’examen de stabilitE d’une couche limite laminaire d’apres la methode de Tollmien, un systäme 
de formules est deduit, en parlie d’une nouvelle maniere, et resume bien disposement. Des tableaux plus 
precis soul calcules pour les fonctions universelles qui s’y presenient. Quelgues resultats nowveaus sont 
communiques comprenant le cas des oscillations, renforcees et affaiblies. Les resultats de la theorie de 
stabilitE oblenues recemment par d’autres auteurs, surlout par C.C. Lin et D. Meksyn, sont discules. 

PassnBaerca — OTYACTH HOBEIM CIHOCOOOM — H HATIANHO COCTABIACTCA CHCTEMA POopMyıL, 
HEOÖXOAUMBIX MIA MSCHeNOBAHHNÜ YCTOÜYNBOCTH HAMUHHAPHEX IOTPAHHYHEIX CIIOEB METONOM 
Toasmuna. Jlıa Tam BCTpeyamınuxca BCeoömmmx PyHEUMÜ Nalımca TAÖNHINO YHCTEeHHEIX 
SHAyeHNf, OTIIHYaMINMecH IO CPABHeHNW C TIePBOHAYANBHEMN Öonpmei TOYHOCTBP. Ilpn- 
BONATCH HOBBIe PEe3yAIBTATEI UCYNHCHeHNH TIIABHBEIM OÖPa30M MIA CAYYaeB PA3BHBAWINHXCH 
HIH 3ATYyXamIımXx Konebannü. Kpnurtuyecku OCBETIAMWTCH PesyIbTATbI HOBeHIIHX HSCITENO- 
BAHHÜ B O0CODeHHocTHu paöorz II. I. JIuma u I. Mercnma. 


1. Rückblick und Zielsetzung !) 


Die wichtige Frage nach den Ursachen und den Bedingungen des Umschlags der lami- 
naren Reibungsschicht in die turbulente Form beschäftigt die Strömungsforscher seit geraumer 
Zeit. So hatte Tollmien [13] 1929 das Verhalten kleiner wellenförmiger Störungen, die der 
Hauptströmung überlagert wurden, theoretisch untersucht und u.a. die Indifferenzkurve (oder 
Kurve neutraler Stabilität) für den Fall der parallel angeströmten ebenen Platte errechnet, die 
das Gebiet ‚‚angefachter‘‘ Störungen von den ‚„‚gedämpften‘‘ Störungen abtrennt. Die angefach- 
ten Störungen, wenn sie auf ihrem Wege stromabwärts eine genügend große endliche Amplitude 
erreicht hatten, wurden dabei als Ursache des Umschlags angesehen. Andere Forscher der 
Göttinger Schule, insbesondere Schlichting und Pretsch, haben mittels der Toll- 
mienschen Ansätze eine Reihe von Beispielen durchgerechnet, so daß man jetzt über eine 
genügend lückenlose Serie von Geschwindigkeitsprofilen der laminaren Reibungsschicht verfügt, 
deren Stabilitätsgrenzen bekannt sind. Unter ‚Stabilitätsgrenze‘“ wird dabei die höchstmögliche 
Re-Zahl verstanden, unterhalb deren keine angefachten Störungen vorhanden sind. 

Die Tollmiensche Stabilitätstheorie zur Erklärung des Umschlags der Reibungs- 
schicht erfuhr allerdings ziemlich allgemeine Ablehnung, zumal sowohl amerikanische wie 
englische Windkanalversuche angeliefert wurden, die sie nicht bestätigten. 

G.1.Taylor[11]entwickelte daraufhin 1936 eine anders geartete Theorie der Turbulenz- 
entstehung: Der Turbulenz der Außenströmung wird hier die allein entscheidende Rolle für die 
Entstehung des Reibungsschichtumschlags zuerkannt. Sie bewirkt im Druckanstiegsgebiet 
der Außenströmung an geeigneter Stelle in dem Augenblick, bei dem ihre Schwankungsgeschwin- 
digkeit der Hauptströmungsrichtung gerade entgegengesetzt ist, in Wandnähe eine lokale 
momentane Ablösung der laminaren Reibungsschicht, und diese hat den Umschlag dann un- 
mittelbar zur Folge. Messungen des NACA konnten diese Theorie gut bestätigen. 

Trotz dieser Situation unternahmen es zwei Mitarbeiter der Drydenschen Schule, 
Schubauer undSkramstad [10], im Jahre 1940 ‚experimentell genau so vorzugehen, 
wie Tollmien und Schlichting mathematisch vorgingen“ (Zitataus Dryden[l)) 
Mittels eines magnetisch erregten feinen Metallbandes erzeugten sie Störschwingungen vorgeb- 
barer Fregenz und Amplitude, und ermittelten die Weiterentwicklung derselben stromabwärts 
durch Hitzdrahtmessungen. Die Messungen wurden 1941 abgeschlossen. Das Ergebnis war 

_ eine tatsächlich quantitative Bestätigung der Tollmien-Schlichtingschen Rech- 
nungen. | 
Be; Der scheinbar widerspruchsvolle Sachverhalt, daß somit zwei verschiedenartige Theorien 
“des Reibungsschichtumschlages experimentell bestätigt sind, klärt sich nach Dryden [1] 
1) Für einen ausführlicheren historischen Überblick vgl. die zusammenfassenden Berichte von Dry- 


” en[l] und Tollmien [15]. Indenhistorischen Teilen der Li n schen Arbeit [3] werden dagegen entschei- 
dende Phasen der Entwicklung dieses Forschungsgebietes verkannt. 
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auf folgende Weise auf: Ist die Turbulenz des Außenstromes größer als etwa 0,2% der Außen- 
geschwindigkeit — was insbesondere bei Windkanälen älterer Bauart öfters der Fall ist —, so 
gilt die Taylorsche Theorie; neben dem Druckgradienten der Außenströmung ist der Um- 
schlagspunkt dann durch Intensität und Größe der Turbulenzballen gegeben. Bei Turbulenzen 
von weniger als etwa 0,2% des Außenstromes, d.h. also in turbulenzarmen Kanälen oder bei 
der Bewegung durch ein ruhendes Medium, kommen statt dessen die Tollmienschen Er- 
gebnisse zur Geltung. Oberflächenrauhigkeit, Vibrationen, akustische Luftschwingungen oder 
Außenturbulenz geringer Größe sind dann die kleinen Störungen, die entsprechend den Toll- 
mie'nschen Ergebnissen je nach den Frequenzen, die sie enthalten, und je nach der Form und 
Re-Zahl des Reibungsschichtprofils, das sie durchwandern, gedämpft oder angefacht werden. 
Die angefachten Störungen führen dann früher oder später zum Umschlag, die gedämpften 
dagegen nicht. 

Neben dieser experimentellen Rechtfertigung seiner Theorie durch die Messungen der 
Drydenschen Schule hat Tollmien selbst in einer neueren Arbeit die theoretische 
Rechtfertigung der für seine Rechnungen erforderlichen asymptotischen Näherungen veröffent- 
licht [14]. 

Die Stabilitätstheorie erscheint somit in Zukunft als ein sicherer Bestandteil der Strömungs- 

physik. Es ist daher nunmehr auch angezeigt, einige noch bestehende Unvollkommenheiten 
der älteren Rechnungen, vor allem bei den dabei auftretenden universellen Funktionen und ihrer 
Tabellierung, zu beseitigen. Insbesondere soll auch das neuerdings öfters diskutierte Verhalten 
der sogenannten inneren Reibungsschicht bei angefachten oder gedämpften Störungen weiter . 
geklärt werden. 
b Durch die geschilderte neue Situation bedingt, ist das Interesse an der Stabilitätstheorie 
\  verständlicherweise neuerdings auch in Forschungszentren außerhalb der Drydenschen 
Schule erwacht und hat zu einigen Veröffentlichungen zu diesem Thema geführt [3], [4]. Auch 
zu diesem neuesten Schrifttum werden im folgenden einige Ergänzungen bzw. kritische Be- 
merkungen mitgeteilt. 


2. Die Grundgleichungen der Stabilitätstheorie 
Obwohl dies in der Literatur schon sehr oft geschehen ist, müssen auch hier nochmals 
die Grundlagen der Stabilitätstheorie wenigstens in skizzierter Form angegeben werden, damit 
das Weitere verständlich wird. 
Ausgangspunkt der Rechnungen ist die Störungsdifferentialgleichung 
d? d2U  /da d2 
(U—o) (er N 28 +29)  . 
Sie entsteht aus den Navier-Stokesschen Differentialgleichungen in ihrer zweidimensio- 
nalen Form, wenn man darin die (in komplexer Form geschriebene) Stromfunktion einer stören- 


den Partialschwingung i 
(2, yy = glg) er Zen rg (2) 


einführt und noch die folgenden Voraussetzungen macht: y ist so klein, daß linearisiert werden 
kann; die Hauptgeschwindigkeit U verläuft stets in x-Richtung und hängt nur von der dazu 
senkrechten y-Koordinate ab; alle Größen sind dimensionslos gemacht und die Bezugsgrößen 
zu der Reynoldsschen Zahl R zusammengefaßt. « ist reell vorausgesetzt und definitions- 
gemäß das 2r-fache der reziproken Wellen- 
länge der Störung. e dagegen wird im all- 
gemeinsten Fall komplex angenommen: 


Pe ae U 1 FR): 


Nach (2) ist dann c, die Geschwindigkeit, 
mit der die Phase der Störwelle fortschreitet, 
während c; je nach Vorzeichen Anfachungs- 
bzw. Dämpfungsparameter der Störungs- 
amplitude ist. Der Nullpunkt des recht- 
winkligen x, y-Koordinatensystems wird 
üblicherweise nicht an die Wand gelegt, 
sondern an die Stelle, an der die Geschwin- 
u der Hauptströmung U gleich der 
Br: h ?hasengeschwindigkeit c, der Störungswelle 
En . N 1). U(y) wird dann zweckmäßigerweise von diesem Nullpunkt aus als a: h en 

U)=o+Usy+gp+... 


Bild 1. Zur Festlegung des Koordinatensystems 
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Bei der Lösung von (1) zu vorgegebenem U(,) handelt es sich nun um ein Eigenwert- 
problem, indem nämlich nur für bestimmte Scharen der Parameter &, R, c,, c; Lösungen vor- 
handen sind. Zu jeder der möglichen Parameterkombinationen kann dann die Amplituden- 
funktion p(y) der Störungswelle errechnet werden. Die Besonderheit des Lösungsweges zu (1) 
beruht nun darin, daß die 4 erforderlichen Partikulärlösungen nicht aus (1) direkt gefunden 
werden, sondern aus Näherungsbeziehungen von (1): Man macht z. B. zunächst heuristisch die 
Annahme — die sich dann aber durch die Rechnung bestätigt —, daß es 2 langsam und 2 schnell 
mit y veränderliche Partikulärlösungen von (1) gibt. Setzt man dann noch R genügend groß 
voraus, so kann für die langsam veränderlichen Lösungen die rechte Seite von (1) annähernd 
gleich Null gesetzt werden. So entsteht die „reibungslose‘“ Differentialgleichung 2. Ordnung 

d? d?U h 


(U— c) (er) art es I DEN N ER de (5) 


und aus ihr zwei Partikulärlösungen 9}, 9%. Diese haben die Form 
x 
P3: Yy Dan y" 
n=Uu 


oo Lu 
Pr Bo I v’ "Play 


wobei die hier auftretenden Summen reguläre Potenzreihen sind. Insbesondere ist 
PORN E bi =N AT . “nt are EL . (7) 


gesetzt, während allgemein die a,, d„ Funktionen der vorgegebenen Geschwindigkeitsvertleilung 
U(y) sowie der in (5) noch verbliebenen Parameter sind. Es sei hier noch angemerkt, daß für 
ra Anfachung oder Dämpfung der Ansatz für p, abgeändert werden muß. Darüber mehr 

5. Abschnitt. Der Ansatz in (6) gilt also in der Hauptsache für den neutralen Fall 4 =. 
Bei ER beiden schnell veränderlichen Partikulärlösungen, 9, und 9, genannt, kann die rechte 
Seite von (1) nicht vernachlässigt werden, dafür aber die niedrigeren Differentialquotienten von 
p sowie selbst gegenüber den höheren. Man erhält dadurch die Näherungsbeziehung 


a Dar: 
ao U—) T=0. N TED ON: 
Führt man mit Tollmien [14] statt % eine Koordinate Y durch die Ansätze 
av_1/0=e 
dy Die Y 


N) 


El Il | | 


ein, so wird (für c; =0) dadurch offenbar ein monotoner Zusammenhang zwischen Y und y 
gesichert. Mittels der Reihenentwicklung (4) für die Geschwindigkeitsverleilung wird dann 


1 1 17% ( 1 U . P da ) i 

; = — 4%? 377 1% 02 N RE LP 

Be: Setzt man nun noch zur Abkürzung 
ee. 4. (11) 


und führt eine Funktion F(Y) ein durch 

Rp dy dy 
dy® & Yay 
so wird, wie in [14] ausführlich begründet ist, die Differentialgleichung (8) in asymptotischer 
rt. 5 Näherung .ersetzbar durch = 


dy? 
- Der wesentliche Fortschritt gegenüber dem früheren Bestand ist der, daß diese Dilferential- 
gleichung für den ganzen y-Bereich Gültigkeit hat, nicht nur für kleine y-Werte. 
Wir führen für y und Y noch die Koordinatendehnungen ein 


BR Ex 
u dr 


Sa ae. (12), 


—ieYy-F(Y) — 


(13). 
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Durch diese geht die obige Differentialgleichung über in 
d?F 
dH? 


deren Partikulärlösungen mit F,, F, bezeichnet werden sollen. Aus F,, F, ‚erhält man dann 
mittels (12) zwei Partikulärlösungen 93 9, durch 


nn - VRR. ra. 


Der Vorteil, auf diese Weise für Be y-Werte gültige Funktionen 9;, } finden zu können, 


BE .HR— Be a 


muß also damit erkauft werden, daß diese (durch den Faktor |/ 04 TH -/) IY ihren universellen 


Charakter verlieren, d.h. auch, nicht mehrallgemein tabellierbar sind. Indessen wird man in 
vielen Fällen mit der Näherung 
ey SE Rei 20r, N 


auskommen, die ja nach (10) gerade in dem Gebiet um y =0 gilt, in dem die Lösungen 93, 94 
von besonderer Wichtigkeit sind. Nennen wir diese Näherungen von 93, 94 DZW. 931, Pa1, So gilt 
beispielsweise für p;, nach (15) 


Fra _F,0) Rn ren 


dn? 
und damit nach (14) für 9;ı auch die Differentialgleichung | 
d* Ar N URE Fa _ 09 (18). 


dn* d? 


Die Lösungen von (18) werden uns in Abschnitt 3 ausgiebig beschäftigen. Es sei aber nnchmals 
hervorgehoben, daß die genauere Rechnung, d.h. die Berechnung von 9, bei vorgegebener Ge- # 
schwindigkeitsverteilung U(y), ausgehend von der Näherung d?o;,/dn%, über (17) und (15) 
nunmehr auch möglich ist. 

Jetzt müssen wir nochmals auf die reibungslose Lösung 9, für den neutralen Fall zurück- 5; 
kommen. Wie aus (6) ersichtlich ist, wird die 1. Ableitung von 9, für y= 0 unendlich. Dies 
5 ist offenbar der mathematische Ausdruck dafür, daß bei Abwesenheit von Reibungskräften 
B: ein bei y = 0 strömendes Teilchen eine zeitlich beliebig anwachsende Geschwindigkeitsänderung 
erfährt, da es ja einer Störung mit zeitlich konstanter Phase unterliegt. Um den Sachverhalt 
annähernd richtig wiederzugeben, muß daher p, in der Nullpunktsumgebung eine „Reibungs- 
korrektur“ erfahren. Da nach (6) 


a EN a N zur. (UO 
ist, macht Tollmien [13] für die korrigierte Funktion den Ansatz A 


- n Vie ht Eraser A N (20) 
(ohne daß dadurch eine Entwicklung nach e angedeutet werden soll), der ale in einem schmalen 
Streifen der x, y-Ebene beiderseits y=0 Gültigkeit haben soll. Hiermit BL man in (1) ein, 
ersetzt 9, sowie d2U/dy? durch seinen Wert im Nullpunkt (nach (19) und (4)) und U(y)durch 
‚seine Potenzreihe (4) bis zum linearen Glied. Macht man ferner die a nach Ws 
(13), so erhält man mit. Tollmien als Differentialgleichung für die Korrekturfunktion Der 


Wo ; ® Furt m »D., 
u = IT a se 
Die Größenordnung von e muß offenbar so sein, daß für kleine %, d. h. solange in (6) das logarith- & 


mische Glied noch die beherrschende Rolle spielt, 7 bereits sehr groß ist. Unter diesen Um- & 
ständen erreicht man mittels der Vorschrift Er 


a: n. 


„ [22 


U, , 
Br Trage Bu nIn(e'n) m»). 


als Randbedingung für die Lösung von (21) bei großem nıin der Tat , (20) a 
an die reibungslose Lösung (6). Diese Lösung von en ae im 4. Abschn itt näher 
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Pr ist es Ti 0 11. mien [14, dort @ E sogar gelungen, eine für de ganze Integrationsgebiet 
Pe ültige Partikulärlösung an Stelle von , herzustellen. 

# Wir vermuten, daß es sogar Strömungen gibt, bei denen die Einführung der im ganzen 
trömungsbereich gültigen Funktion ®, zur Durchführung einer Stabilitätsrechnung notwendig 
ist. Die Einführung von 9, durch (23) hätte deinnach nicht nur die Bedeutung einer theoretisch- 
ne, Vereinfachung des Problems, sondern wäre ein wirkliches Erfordernis für manche 
_ Anwendung. y 


3. Die Reibungsiösung Ps 
a) Die 2. Ableitung von 9a 


Nach [14] sind zwei Partikulärlösungen der Differentialgleichung (18) 


4 - m n EX d? 
Bi x a N 5 (en) je) = Fin 


. RN Tee - a2 PIE (N) je) = F,(n) 


Dabei sind die N) (2) die Hankelschen Funktionen 1. bzw. 2. Art von der Ordnung 1/3. Da die 
‚allgemeine Beziehung besteht 


en ER RE, SEE ey; 
die, rechten Striche bedeuten hier die konjugiert komplexen Werte), genügt es, die 
mathematische Diskussion auf 95, (n) zu beschränken. 


Br Nach Watson[18, p. 74] gilt für die Hankelsche Funktion 1. Art von der Ordnung 1/3 
Er die im Endlichen überall reguläre Potenzreihenentwicklung 


:: i j | h DEN 1 \or 1 \9# EIER RAN 
= ie (2) = . - en (— 1)" 57 ; , a ; (28). 


1 1 
a) N ee a 
in | Peatı 3) TIn+1+3) 


(24). 


Hiermit ergibt sich aus (24) nach einigen Umformungen 


de 
I =F, (n): 2 - Ir 
wobei für die en fürn=0 


1 
Sa (3) 2. 316 


at RK 
Sure zn nlaı +1) g Aın En e 3% Asn nn "| SPA (27), 


un Ta N As ri (28) 
ae (3) 
D. ist, während für nZ1 die Rekursionsformeln gelten 
4 #-; n. k 7 1 Dan | : 1 . h ’ 
An = a. A 1,n-1> Ads, = 3n-+1 N A ae Re (29) 
Tı nt. man Roc in (27). 'Real- und Imaginärteil, so wird daraus 
Be “ nv he \ N 5 x 
RTL ZI ERE SER N SERIE, ) 
; Man 22 2 1)" a n n (Aun — sin ’D Ann | 
Ne: gerad | . . (30a) 
x RN. (30b). 


Fe kn 


der komplexen Zahl 


1 
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Mit Hilfe der Formeln (30) sind in Tabelle 3 die Werte der 2. und 3. Spalte für den Bereich 
—5<n<s +5 gerechnet. Mit wachsendem |n| wird die Rechnung nach diesen Formeln zu- 
nehmend mühsam. 

Es gibt indessen eine asymptotische Entwicklung für die Hankelschen Funktionde die 
gerade für große || günstig ist: Nach Watson [18, pp. 196 sq.] ist nämlich 


as ar © Di es a... a1), 
13 ( T*2 : 2 AN = 


Nn=U 
wobei für die Konstanten , gilt 
B=1, | 


En. 1 - 
(31) ist nun allerdings nur eine semikonvergente Reihe: Ihre Glieder nehmen zunächst ab, dann 
aber unbeschränkt zu. Sie nähert die Hankelsche Funktion am besten dann an, wenn man be’m 
jeweils kleinsten Glied abbricht. Der Fehler ist dabei von der Größenordnung des Betrages 
des kleinsten Gliedes. Er wird ersichtlich um so geringer, je größer der. Betrag von z ist. Bei 
festem Betrag von z wird er um so größer, je mehr das Argument von z sich den Grenzen der 
Ungleichung 


191 On ern, tür nzı) | 


EN 


nähert. An den Argumentgrenzen selbst ist die Abschätzung eines endlich großen Fehlers nicht 2 
mehr möglich. m 
Schon hier sei ausdrücklich betont, daß die Heranziehung der asymptotischen Entwicklung 1 
(31) keineswegs eine physikalisch bedingte Notwendigkeit darstellt, sondern nur der bequemeren 
Berechnung der Hankelschen Funktion halber empfehlenswert ist. Grundsätzlich kommt man , 
auch mit der regulären Potenzreihe (26) aus. Es ist daher nicht angängig, aus den Besonder- ’ 
heiten der asymptotischen Entwicklung, wie etwa der Argumentbeschränkung (33), irgend- 
welche Rückschlüsse physikalischer Art zu ziehen. 
Die asymptotische Entwicklung (31) liefert für 2) die Eezjchune: 
- = 0 n | —(-+zn er +) PR , H 
Te -Am=V2e 24 De Bm AR: 3 | a+on NE 
TU Nn=U Y2 R 


während die Argumentbeschränkung (33), wenn wir uns die an sich reelle Veränderliche RAN 


en ee EEE . (35) 


e) 


erweitert vorstellen, sich zu 


ga<d<la ER a ae +38 5:21» (86) 
transformiert. Diese Beziehung macht es notwendig, den Anschluß zu negaliv reellen ze Kt, 
stets über die untere Halbebene zu vollziehen; also wird z.B. DER nn 


} = 
ner (für 7 reell ne BUT >20 2 A: n. 


Bei Trennung von Real- und Imaginärteil in (34) ist es erforderlicil ade dem Vorzeichen von. 
n zu unterscheiden. Mit den Abkürzungen Be 


Po(n) =DeR nl; Pi Bıln) - 5m RE 


: 
a 
x 
x 


gilt dann, für 720 | E 
PrYsır E A x cos ß, — sin ß, » er”, 
ee ef .e (— 1er» Br 1 Se 17. ae 
dn ” 2 IR Y2 na 
2 ad sin Pı° ae Bun an+2 * a enpg Ee Pı + Bunzs : cos E 


d2Qpz1i 3 mr 2 


N=U 


+ Buarı ° cos ßy°n yR — Bun n+2 * 


Erz 
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nd für n<O ist 


d? P31r j 3 
en. 27 Pa er. n| In, fi 
1? Imj1/4 - eP. | 2 (— 1)" 6 Pan cos ßı | 


+ Bunzı a ml ’?-HBin+zsin Pi In®—Bun+s en we | 


d’dzni 13, “ R OL In. re i 
en lie Ina. e# X 1)” In]? B;n sin Pı 
TAT . (394): 


cos sin cosß,-+sin tan 
er BD; n--1 u jn | 3/2 4 Bun+ -9 COS ßı' Inl®+ Bin+3 al 2) 

Mit den Formeln (30) und (39) sind in Tabelle 3 die Werte der 2. und 3. Spalte gerechnet. 
Eine direkte Auftragung der 2. Ableitung von 9,, über den ganzen tabulierten Bereich ist wegen 
der großen auftretenden Ordinatenunterschiede nicht zweckmäßig. Daher wurden die Funktions- 


werte zunächst mit dem Verzerrungslaktor 
a ee ee len in. (40) 
multipliziert, der die Wertunterschiede genügend reduziert. Die Zahl 1,3 im Exponenten von 


(40) erwies sich nach einigem Probieren als günstig für eine übersichtliche Auftragung. Das 
Ergebnis der Auftragung zeigt Bild 2. Einige Werte des Verzerrungsfaktors sind unten im 


V 
el 
2 3 4 
10 10 10 
Bild 2. Darstellung der mit einem veränderlichen Faktor f (n) verzerrten 2. Ableitungen von 9, (n). 


r Zahlenwerte von f (n) sind am unteren Bildrand angegeben 


Bild mit angegeben, so daß auch die Größenordnung der wirklichen Zahlenwerte der zweiten 
Ableitung von p,, daraus entnommen werden kann. 


b) Dieerste Ableitung von 935, 


Aus physikalischen Gründen muß verlangt werden, daß für große n-Werte die Störfunktion 
$3ı sowie ihre Ableitungen nach Null Fr Daher ist zu setzen 


Is u d2p51(£) 
31 En a N AR (4). 


2 Integrationsweg soll hier die reelle Achse gewählt werden. |n,| sei dabei ein 
ckmäßig gewählter Abstand vom Punkt 7 = 0, der auf der positiven wie auf der negativen 
den Bereich der Potenzreihenentwicklung (30) von dem der asymptotischen Ent- 
(39) für den Integranden scheidet. 
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Beginnen wir mit dem Bereich 
©2nz nl 
des Integrationsweges, so stellt sich also hier die Aufgabe, die Entwicklung (34) zu integrieren. 
Hierzu benutzen wir die Identität 


yl-+a—c a : 
iE .eede— ec. ERER ich 3 E21 De Fe 20) aita-8e 


\b’-ce (b-.o% (b - c)® 
1). REN (1+a—o)(1 art (1 -+a—ne) arena! Wi % 
! ip dta—=g +a—2%)-  (1+a—(n-+1) )e), Ba 
HD) KM en. 


. ee “dx const 


Bei ce <0O erhält man für das Integral auf der linken Seite dieser Beziehung offenbar eine 
Näherung, wenn man den Integralausdruck rechts unterdrückt. Dabei läßt sich der erforder- 
liche Gütegrad der Näherung durch die Wahl eines genügend großen n erreichen. Wendet 
man diese angenäherie-Integration auf die asymptotische Entwicklung (34) an, und ordnet 
nach Potenzen von 7], so wird 


Lpzıln) Sn > & et In a | 
—_— O0 — 1 _- ea. Mt >' ERTL | een E z 3 
d n TL e e 7 E= ( 1) 5) D, 7 R (43). 


Nn=V ” 
(gültig für o2nZ |nı|) ‚J 
Dabei sind P, (67); Bı (m) die schon in (38) nl Abkürzungen, während 
>" Ba ER 1—u te Near ned a Kanal genaue wor (44), 
0 für »sO0 Er | 
; | N 
Wr ao (45) 
Een +1) 2a +3 2a +5) -(2u +2» —3)) für ‚>2 
ist. Trennt man in (43) noch Real- und Imaginärteil, so wird 
doa1r 3 \ 2 i 
A: “ TE arg > (ae ee Dun ‚sin By Dinzı een are 
u 7 | Y 2 
I ... (46a) 
+ Dyn+e * cos Bus — Dyn+3 ud N Pı pe] | 
V2 J 
dpsıi 3 & cos ß, — sin 
m” VE ar end nen [Dun >05 9, — Dim Sea ve | 
n=U | Y2 . (46b) 


Dasein pur + Dan FE zer (garig fr »an> 


Für den anschließenden Teil des Integrationsweges, |7ı!>2n>=— |yıl führen wir die 
Hilfsfunktion ein 


d? 4 i 
In) = De: (gültig für mla2n2— m) »-- :»- (47). 
: (N) 
Die Integration von (27) liefert dann 
EI tan EHEN A | 
I)=— Intl, 12 ©, any 
was nach Trennung von Real- und ah, ergibt | 
1 4A A 
a) DV anni (Lan sn ©. Ass ) 
a den PIEIH In+1 ana 
gerade 
© Br 1 | 
— Zn? Aura An 2. Au anl 
2 ni ar Snt1 9° 72 3n+2 7 


ungerade 


Yus 
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A EBEN. 
Ne Intl. ın Be ER 
An) 27 DrTEE e- per | 
| Fr ; N (40B). 
# n® sur. („San — Bar Am: ) 
R +25 ee an+2 
& i nühgerade 
In Tabelle 1 sind in Spalte 2 und 3 Wertepaare dieser Funktion für den Bereich |] <5 wieder- 
R gegeben. Für die erste Ableitung von 9,, ergibt sich dann in diesem Bereich nach (41) und (47) 
Ilm) ; 1 EN 
ln CRM — Ilm) +Im) (gültig für plan 2m) 60. 


wobei die Konstante (dpsı/dn)un,n mach der asymptotischen Formel (43) gerechnet wird. 


-Tabellel. Hilfsfunktionen fürg,(m) im Bereich 
‘der Potenzreihe (27) 


ER 2. | 3. 4. 5. 

- an I 5m I nm | Mm | Min 
B: Bei: —56,6495 | — 9,2087 | + 3,5089 | -+21,5517 
3 —,5 | —ı5,2712 | —23,8585 | —13,6281 | +11,5117 
Br‘ 1 + 4,2542 | —15,2208 | —15,5449 | + 1,4109 
Ei 3,5 + 8,8977 ' — 5,3783 | —11,8510 | — 3,5991 
Dr Ex + 7,5882 | + 0,2870 | — 7,6128 | — 4,6927 
Rt: a5 + 5,1170 | + 2,4221 | — 4,4390 | — 3,9035 
2 2 + 3,1319 | + 2,6328 | — 2,4074 | — 2,5880 
Ri 278 + 1,8448 | + 2,0491 | — 1,1888 | — 1,4000 
BR 1 + 1,0356 | + 1,2617 | — 0,4831 | — 0,5712 
% —0,5 + 0,4665 | + 0,5407 | — 0,1138 | — 0,1265 

) (0) 0 0) 0 
+0,5 | — 0,4040 | — 0,3287 | — 0,1037 | — 0,0909 
F x 221 — 0,1312 | — 0,4685 | — 0,3912 | — 0,2970 
| SER — 0.9595 | — 0,4776 | — 0,8183 | — 0,5376 
j 9 — 1,0873 | — 0,4241 | — 1,3338 | — 0,7643 
| +2,5 — 1,1370 | — 0,3619.| — 1,8925 | — 0,9604 
air — 1,1421 | — 0,3187 | — 2,4634 | — 1,1295 
Bi % +3,5 .| — 11314 | — 0,2987 | — 3,0320 | — 1,2830 
® +4 — 1.1210 | — 0,2941 | — 3,5949 | — 1,4308 
Ei +45 | — 11158 | — 0,2956 | — 4,1539 | — 1,5781° 
+5. —_ 1.1145 | — 0,2977 | — #7113 | — 1,7265 


Für den letzten Teil des Integrationsweges, also für —||2n=— ©, muß wieder die 
asymptotische Entwicklung (34) herangezogen werden. Führen wir hier überall mittels (37) 
|n| statt 7 ein, so können wir die Hilfsfunktion definieren | 


In] 
L(nı}) = IKACKE PREV MR EINEN 


R Die Integration mittels der Identität (42) ergibt dann mit den Abkürzungen wie bei (43) 


(4 
i a TE 2/1 — a\3+3R RE U 
Lin!) za ei. CR y>i Ba en Ben BE er 102). 


” 


Bis v: N a) \ "| cos B,+ sin Pı 

 EtmeV 2: mer er > 18 - Im Du 

“ ne a Bi 2,‘ Fe Be." y2 (53a). 
— Dans sin Bı* ml ®— Dante =... Iml-® + Dinss* cos Bı Rn 


“ & 


2. u l.ugsäb). 


Be ARE sin A 
*.c08 Bı ml + Din+2* ze Br Im| —D;n.+3 ' sin ßı° Im] u 
{ nr nn; re / gt er k PR . 
Dee we . a. Br k Pr 3 
ERE Pu, En i 
& , e 


R 
Gr 
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In Tabelle 2, Spalte 2 und 3 sind Wertepaare dieser ah ktion für den Berch Ms: 2n2— 8 
errechnet. Für die erste Ableitung von 93, ergibt sich dann nach (41) und (52) r. 


{ dv (Fe) HLmd Elm) (gültig für mlanz—o) MN, 
dn \ dn/ın) 
wobei die Konstante (dpz1/dn)— m) aus (50) berechnet wird. E 


Tabelle2. Hilfsfunktionen für g„(n) im Bereich 
der asymptotischen Entwicklung (34 


% 
N 2. | 0: 4. MN) 6, 

n | zn) | Zn) | Irland | Nithn) 

F —d + 58,86 + 9,81 + 13,65 + 25,12 
—5,5 + 98,49 — 75,74 | + 55,00 + 12,72 

, —6 + 16,01 — 270,0 + 93,43 — 70,05 
—6,5 — 453,8 — 414,1 + 6,44 — 253,2 
Se —1424 + 265,5 — 452,0 — 353,1 
A — 1485 -+3233 —1293 + 400,7 
—$8 +4784 +7679 — 906,4 +3195 


Damit sind die Formeln für die Berechnung der ersten Ableitung von @;, ermittelt. Mittels 
(46), (50), (54) wurden Zahlenwerte für do;,/dn im Bereich —8<sn<s +8 errechnet und in 


N Tabelle 3, Spalte 4 und 5 eingetragen. 
A Für die graphische Darstellung von dQ;,/dn wurde wie bei Bild 2 wieder der Verzerrungs- 


faktor f aus (40) mitverwendet. Bild 3 zeigt das Ergebnis. 


Bild 3. Darstellung der verzerrten 1. Ableitungen von 9, (n). Verzerrungsfaktor wie bei’ Bild 2 
ec) Die Funktion 9 { 
Da 9;, für große n nach null en ee muß ähnlich an) gez a 


Pzıln) = 


[art & 951 (8) de. 


Auch für diese Integration werde die reelle . als Ian 


Bereich ., 
oZnZ2hm| u 


ist dann die asymptotische Entwicklung (43) zu int er n. S 
u; Beantitet (42). Die Rechnung ergibt 5 2 z = 


P51(N) 17 Ba .e-BiM. und. > = 


2 


und 


MR u 1 


N—U 


cos 9, — sin 
Bae cos ßı 7 ea Sn 


n 


Pam ei 5/4. ee’ » > (— 1)” i Da 


Nn=U 


—Guntısinßrn 


a den anschließenden Bereich |nı| zn — |nı| führen wir die Hilfsfunktion ein 


Die Integration von (48) liefert dann 


oo 


NM (n) I ea HF BnEnEnF)) 


was nach Trennung von Real- und Imaginärteil Setbe, 


N no u 


m 


ungerade 


Bau: - => em ei 


a ur) — 


' Trennt man in (56) Real- und Imaginärteil, so wird 


Ps1ır (N) DE m: es» E> (—1jr mr: Ga 


Bl) —1) 


‚eos ßı+sin P, 


cos ß, — sin ßı 

1 
cos sin 5 

2 Gunte A n3—Gan+s cos ßı 7 i | 


(gültig für e>2n2 |n!). 


BEL u Eee ec 6 Par, Taf wie 


a) - ns: 


1 
LE Eee 


a ala > | 
Ent) (En+3) "| 


Be 
en DE “ Re: en+d ent)" | 


ER ne = ne eu lie Dr ET 2n . 
Gn+1) "RER 12 8n+2)(an+3) 


en) (8n+2) re Kerze] @nz3y" 


(58). 


(59a), 


. (59b) 


(62a). 


(62b). 


In Tabee 1sind i in Spalte 4 und 5 Wertepaare dieser onetion für den Bereich In] <5 Matten 


Fi u JE - 


eben. ‚Für: Baugelhat Set sich dann in diesem Bereich a (55) und (60) 
NM (\n)+ Mn) | 


— m] -( (ml m) — 


(ig für lem) 
Formel (56) berechnet wird. 
e 


s, also für —|n,|2n2 —®, muß wieder die 
ng 6 herangezoge werden. Wir definieren dazu die Hilfsfunktion 


ir P- Pi \ 2) =l 


N a inüschenn 
Teil des Integrationsweg 


. (63), 


wa 7 N. “ 
ee FE 
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Die Integration mittels der Identität (42) ergibt dann mit den Abkürzungen wie bei (56) 


j PR] \2+n a 
N (|n]) -Y2 eihim em Imj-5%% Sn vs) Gy’ Ulla N os) 


N=U 


(gültig für — |n1|272 —) 
Trennt man hier noch Real- und Imaginärteil, so wird 


[3 a 5 cos ß, — sin | | ES, 

N,(|n))= Ye 2 Narr Dr nl-sn Gun sin Bı — Gun+ı Pr\,j-or ; 
| m ” ar 

— (Uun+2 - cos Bı nl? — Gunts cos ßı+ sin Pr | por \ er 


3. BR rl 3 
5 S cos sin ßı , | 
GER Dr (17 aeg Gr 008 Bı+ @untı neinfı Ulak | | a 
| > (Gb). 
cos Pı—sin Pr yjsel | %. 
y2 : a 
In Tabelle 2, Spalte 4 und 5, sind Wertepaare dieser Funktion für den Bereich -52 7,2 —8 
errechnet. Für @,, ergibt sich dann in diesem Bereich nach (55) und (64) 


ir 


3 
8b = Ein 


N=U 


+ Gun » sin Br In? — Gun+3 


t 


d 
Ps1(n) = P51(— mil) -(e) " + Zn.) (nl — ml) —N (ml) HN (ml) . . (67), 
wobei die Konstante 95, (—|n1|) nach (63) berechnet wird. 

Damit sind die Formeln für die Berechnung von ,, ermittelt. Mittels (59), (63), (67) 
wurden Zahlenwerte von 9,3, im Bereich —8<n<+8 errechnet und in Tabelle 3, Spalte 6 
und 7, wiedergegeben. % 

Für die graphische Darstellung von @;, wurde wie bei Bald 2 und 3 wieder der Verzerrungs- 
faktor f aus (40) mitverwendet. Bild 4 zeigt das Ergebnis. 


BER 


Tabelle 3. Die Reibungslösung 9, nebst 1. und 2. Ableitung ig 
1. 2. | A 4 | 5. | 6. \ 1. 
de? 2 r . ! 
n Psır d? pzıi doysır dopsıi Dar | Psıi 


r N. 
x rk, 


0,8395 | - 
0,2989 | 
0,0299 


Me a a BI u ann 5 U U 


0,000034 + 4 0,000012 
0,000011 007 | — 0,000001 
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Bild 4. Darstellung der verzerrten Funktion 9, (n). Verzerrungsfaktor wie bei Bild 2 


d) Die Tietjenssche Funktion D(n) 
Die vorstehend wiedergegebenen Formeln und Tabellen für die Funktion o,, sind erforder- 
lich, um den Verlauf der Amplituden der Störschwingung zu ermitteln. Für die Bestimmung 
der Indifferenzkurve dagegen benötigt man nur den Verlauf des Quotienten 


Pam _p D,+iD 68 

— deos/dn (n )= u [ \ ( ) 
für negative n-Werte. D(n) wurde zuerst von Tietjens [12] tabuliert. Nach ihm haben 
noch andere Autoren [S, 6, 3, 4] genauere Tabulierungen, teilweise auch in anderen sich z. T. 
gegenseitig ergänzenden n-Bereichen, vorgenommen. In Tabelle 4 sind diese Tabulierungen 
nebeneinandergestellt. Dabei ist der Vollständigkeit halber D (n) auch mit unseren neuen, aus 
Tabelle 3 zu entnehmenden Werten, nochmals berechnet worden. Bild 5 zeigt die für die Sta- 
bilitätsrechnung benötigte Auftragung von D(n) in Form eines Polardiagramms. Die Kurve 
ist mit Hilfe unserer neu berechneten Werte gezeichnet. Daneben sind die von Tietjens 
berechneten Werte eingetragen. Für große negative n-Werte gilt, wie man leicht aus (56) und 
(43) entnimmt, 


Da rDin) 2. na BER EN. (69), 
| Man sieht, daß dieser — dem ee Abschlußteil in Bild 5 entsprechende — Bereich von 
nn bei unserem letzten tabulierten Wert n = —8 nahezu erreicht ist. 


Hier seien noch einige kritische Bemerkungen angefügt: Gelegentlich einer Kritik von 

- Ulrichs [16] Stabilitätsrechnungen an Absaugeprofilen erwähnen drei amerikanische For- 
scher, Hahnemann, Freeman und Finston [2], eine ‚von Tollmien und Schlichting 
Be: angegebene unrichtige Kurve‘. Obwohl diese Bezeichnung an sich nichts davon ahnen läßt, 
möchten wir doch der Sachlage nach vermuten, daß dabei die Tietjenssche Funktion mit 
den von Tietjens selbst stammenden Werten gemeint ist, die Tollmienund Schlich- 
5 ing in ihren ältesten Arbeiten zur Stabilitätstherorie auch benutzt haben. Die Verfasser [2] 
_ sprechen dann die Vermutung aus, daß alle mittels dieser Tietjensschen Zahlenwerte ge- 
h fundenen Rechenergebnisse — eben weil diese Zahlenwerte „unrichtig‘‘ seien — von vornherein 
schlechter sind als die Stabilitätsgrenzen, die mittels einer neuerdings von Lin [3] abgeleiteten 
tformel unmittelbar ohne Stabilitätsrechnung gefunden werden können. Offenbar ist den 
'erfassern [2] dabei entgangen, daß Tollmien [13] schon in seiner ältesten Arbeit mit den 
lie Y ensschen Werten für die, ebene Platte die gleiche Stat äts ae R* = 420 erhielt 
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Tabelle 4 Die Tietjenssche Funktion D(n) nach verschiedenen Autoren 


L. angew. Ba ie Dr 
Nr.1/2 Jan. RR, 1950 


Dr (n) Di(n) 

>| 12,821 SA: 5. 101723. 3. 4. 5. 6. 

0 | 0,702 | 0,672 670 | 0,6724 | -0,425 | -0,387 0,387 | -0,3882 
-0,5| 0,785 | 0,770 0,769 | 0,7698 | -0,411 | -0,380 0,381 | -0,3808 
-1,0 | 0,920 | 0,892 0,89161 | 0,892 | 0,8916 | -0,389 | -0,350 -0,35025 | -0,350 | -0,3503 
1,2 0,94763 : -0,32772 
-1,4 1,0076 -0,29698 
-1,5| 1,043 | 1,023 1,038 | 1,0389 | -0,297 | -0,261 0,278 | -0,2780 
-1,6 1,0709 ur: -0,25614 a 
28 1,1366 -0,20293 
-2,0 | 1,206 1,202 | 1,2029 | 1,202 | 1,2029 | -0,147 0,185 | -0,13482 | -0,135 | -0,1348 
2,2 1,2672 | -0,04897 | \ 

2,4 1,3255 -0,05748 
-2,5| 1,357 | 1,358 1,350 | 1,3507 | 0,108 | 0,124 0,119 | 0,1193 
2,6 1,3721 -0,18728 
28 1,3987 0,34216 
-3,0 | 1,400 | 1,397 1,3937 | 1,393 | 1,3937 | 0,515 | 0,519 | 0,52173 0,522 | 0,5217 
3,2 | 1,3423 | 0,72064 
234 1,2277 0,92456 
-3,5 | 1,180 | 1,139 1,142 | 1,1424 | 1,130 | 1,023 1,020 | 1,0202 
-3,6 1,0369 1,1054 
-3,8 0,77338 1,2209 
-4,0 | 0,460 | 0,493 0,47208 | 0,480 | 0,4722 | 1,250 | 1,214 | 1,2288 | 1,224 | 1,2289 
4,2 0,19732 | 1,11558 
4,4 0,01056 £ 0,91568 
—4,5 |-0,0405 -0,0413 | 0,8080 0,8196 0,8034 
: 0,09936 0,66585 | 
4,8 0,07090 0,47400 
-5,0 | 0,0057 0,0056 0,3645 | 0,3643 | 0,3644 
-5,5 | a | | 0,196 | 0,2393 0,2286 | | 0,236 
28, | 0,292 0,2552 0,284 
| | | | | 0,301 0,2927 IR | 0,2946 
7,0 0,2767 0,3021 0,3032 
1,5 0,2554 0,2901 
8,0 0,2447 0,2738 
1. Tietjens 1922 [12], S.19. - 4. Lin 1945 [3], p. 140, Tab. 1. 
2.Schlichting 1935 [8], S. 73, Tab. 2. 5.Me ksyn1946[4], p. 400, Tab. 3. 
3. Pretsch 1941 [6], 8. 62, Tab..3. 6.Holstein 194. 
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wie Lin [3] mit Benutzung seiner sehr genauen Tabulierung von D(n). Daß Schlichting 
[7] mit den gleichen Werten von Tietjens für das gleiche Problem statt dessen R* = 575 
errechnet hat, kann also sicher nicht an den Ungenauigkeiten der Tietjensschen Zahlen- 
werte liegen. Schlichting [8], Lin [3] und Meksyn [4] bezeichnen übrigens überein- 
stimmend die Abweichnug ihrer Tabellierungen für D von der Tietjensschen als „un- 
bedeutend‘. Auch die Betrachtung von Bild 5 vermittelt diesen Eindruck. Wenn Schlich- 
ting, der ja wohl die bisher größte Zahl von Anwendungsbeispielen zur Stabilitätstheorie 
durchgerechnet hat, in einer späteren Arbeit [9] doch wieder die alten Tietjens schen Zahlen- 
werte zugrundelegt, obwohl er inzwischen selbst bessere (vgl. Tab. 4) angegeben hatte, dann 
darf man doch wohl annehmen, daß er die Unterschiede zwischen beiden Tabulierungen als 
praktisch belangslos im Rahmen der Rechengenauigkeit erkannt hat. 

Ohne also die Möglichkeit der Zuverlässigkeit der Lin schen Faustformel in Frage stellen 
zu wollen, müssen wir daher doch nach obigem wohl den Schluß ziehen, daß man für diese kein 
Kapital herausschlagen kann, indem man die Zuverlässigkeit der mit den Tietjensschen 
Zahlenwerten vorgenommenen Stabilisätsrechnungen von vornherein anzweifelt. 


4. Die Reibungskorrektur 9sı im neutralen Fall 
a) Die 2. Ableitung von 95 
Zur Beseitigung der Konstanten U”/U; in der Differentialgleichung (21) führen wir ein 


Mi U‘ 3 
7 Dias ge Yaı (70). 
0 
Dann wird aus (21) 
d? 951 U P5ı 7 
rer = 3 1 
dn* dn? er 
und aus der Randbedingung (22) 
Pırn'in(e:n) En Er TEE ER SEN 
Für die 2. Ableitung von 9;, ergibt sich dann nach (72) die Vorschrift 
0 Pal 1 er 
Rn n>1) (73) 


Die zu (71) homogene Differentialgleichung ist identisch mit der Differentialgleichung (18). 
Daher erhalten wir mit Tollmien[14, S. 71, (8)]als Lösung von (71) für die 2. Ableitung von 


P51(n) 


n n 
Al 
a BE [ro de—F,(n) IEAG a. Ara), 
2 6 | [ee a | 
durch die auch die Randbedingung (73) erfüllt wird. F, und F, sind dabei wieder die durch (24) 
und (25) definierten Funktionen. Man könnte demnach grundsätzlich aus den Formeln des 
vorigen Abschnitts für 9;, nun 9;, mittels (74) und (25) formelmäßig ausdrücken. Einfacher 
gelangt man dazu jedoch auf einem direkteren Wege: Man transformiert die Differentialgleichung 
(71) auf die Veränderliche 1/7 und geht dann in diese mit einem Potenzreihenansatz ein. So er- 
hält man die asymptotische ee für is n 
Ip ee In — 1)! ie 7R 
DR: RTTasI) en+1), ER Een is ei 70), 

wobei ersichtlich die a (73) befriedigt ist. Die rechte Seite von (75) ist wieder, 
wie wegen des engen, aus (74) ersichtlichen Zusammenhanges mit den Formeln für o,, auch 
unmittelbar einleuchtend, eine semikonvergente Reihe, die für die Darstellung der Funktion 
d° 9,,/dn? nur dann brauchbar ist, wenn man beim kleinsten Glied abbricht und einen Fehler 
von der Größenordnung dieses Gliedes als zulässig ansieht. 

Zur ‚Erläuterung der natürlich auch hier bestehenden Argumentbeschränkung für n er- 
weitern wir wieder nach (35) n ins Komplexe. Nach Watson [18, p. 198] gilt nun für F, die 
ei eehränkung 1 

a gr< v = 


so daß sich für die asymptotische I; (75) nunmehr wegen (36) die einschneidendere 
Argumentbeschränkung 


Ze 1 
re Se > > eo 25176) 
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Geht man andererseits in (71) direkt mit dem Potenzreihenansatz 


d?o; < 
Pr u De 7) 
i N=U 
ein, so bleiben die Koeffizienten @,, a, olfen, während 
D 
(4, = Te 
Fe 0 re (78) 
n—83 .. 
= 1 ———— f 
An nn (für n>5) 


wird. a, und a, können aus (74) bzw. deren Differentiation nach 7 gewonnen werden: Die 
Rechnung mittels der Formeln (27) und (25) ergibt 


Se NEN RAN a) 2) | 
(Te ee Au 


an a ER Ki; gun En er I ) | 
Fe N — 7 er As ke: 12 din ARE sın iD dn CN, 


Somit sind auch die Koeffizienten der regulären Potenzreihe (77) bestimmt, die der Potenzreihe 
(27) für @,, entspricht. Wegen der Mitbenutzung von (74) ist auch die Erfüllung der Rand- 
bedingung (73) gesichert. Mittels der Koeffizientengleichungen (78) und (79) läßt sich (77) in 
der Form schreiben 


co 


do; a: ak 
a BE a RL (80 
P=0,1,2 
mit 
| [5 
lasntol = — (PH) (PH P+N:-(p+3n— 2! (n21; p=0,1,2) (81). 


Ontp)! 
Die Schreibweise (80) hat den Vorzug, daß aus ihr die Aufteilung in den reellen und imaginären 
Anteil unmittelbar ersichtlich ist. 

Wie man (75) und (80) entnimmt, ist der Realteil der 2. Ableitung von 95, eine ungerade, 
der Imaginärteil eine gerade Funktion in n. Es genügte daher, diese Funktionen für „>20 zu 
berechnen. In Tabelle 6, Spalte 2 und 3 sind die errechneten Zahlenwerte für den Bereich 


7 5 } 
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0<n<ß8 eingetragen. Dabei wurde bis 7 = 6 die Berechnung nach der Potenzreihe (80), der 
restliche Bereich nach der asymptotischen Reihe (75) vorgenommen. Bild 6 zeigt den Verlauf 
der beiden Funktionen. 

Anhangsweise sei hier eine mathematische Besonderheit mitgeteilt, die bei der Aufstellung 
der Koeffizientenformeln für die Potenzreihe (80) beobachtet wurde: 

Nach obigem ergaben sich alle Koeffizienten a„ mit m = 2 (mod 3) der Potenzreihe (80) 
schon allein durch Einsetzen des Potenzreihenansatzes in die Differentialgleichung (71), ohne 
alle Randbedingungen. Nach (78) erhält man leicht für diese Teilfolge der Koeffizienten 


n!» 3% 
Ganta SIE a SERRSAT 23.) 
Führen wir hier m=3n +2 statt n ein, so entsteht aus (82) 
- 
Fr EI 
= mr RL NR N... (BBY: 


Man ist versucht zu fragen, ob diese Formel nicht auch für Werte m == 2 (mod 3) anwendbar ist. 
Für m =0 und m =1 ergibt sie nach einfachen Rechenformeln der /-Funktion 
| a, =i- 3723 T'(1/3) = i- 1,28790 | 
D-7 3316 
= — — — = 0,938 
4 = Faß 0,938893 | 


Rechnet man mittels dieser Werte und der Formel (79) rückwärts versuchsweise die ersten 
Ableitungen von 93;, im Nullpunkt aus, und setzt für A,o As0 die in (28) angegebenen Zahlen- 


werte ein, so wird 
) . Bst — 1,115355 | 


ARE 


a $ . (85). 
Kan au Nom hl 0,298859 
dn /n=0 6 i 


Diese Werte sind in der Tat in vollkommener Übereinstimmung mit den entsprechenden in 
Tabelle 3, Spalte 4 und 5 numerisch ermittelten Zahlen, wodurch auch die Richtigkeit der 
Formeln (84) sichergestellt ist. 

Auf die Frage der Verallgemeinerungsfähigkeit dieser eigenartigen ‚Koeffizienten- 
suchmethode“ sei hier nicht näher eingegangen. Wir haben jedenfalls zweifachen Vorteil aus 
ihr gezogen: Einmal ergeben die Formeln (85) genaue Zahlenwerte für 2 Konstanten, die bei der 
Integration der 2, und der 1. Ableitung von 9;, auftreten. Ferner kann die Tabulierung von 
%,ı nun unabhängig von der von 9,, mittels der nach (84) beliebig genau berechenbaren Koeffi- 
zienten 4, 4, vorgenommen werden. 

1 


by Die 1.Ableıtungvon 9,, 
Die Differentiation der Randbedingung (72) nach n liefert für die 1. Ableitung von 931 
die Bedingung 


ea einn+ine+i EEE DR SR Ep 
Diese wird offenbar dann erfüllt, wenn man die asymptotische Entwicklung (75) zu 
| AR: ds =] 1 1 u “\n (In —1)! _8n 8% 
E hr u ren BEN... (87) 


integriert. 

Die Bags, wie das hier wieder auftretende logarithmische Glied bei Vorzeichenwechsel 
von n zu behandeln ist, wurde bisher hier nicht behandelt. Da indessen wegen der Argument- 
€ en (76) die Beziehung (37) auch hier wieder gilt, folgt auch eindeutig 


Inn =In m —i-r u Dam... Seal. (88% 


Damit ist die „‚Übergangssubstitution“ in der reibungslosen Lösung 9, bei y =0 eindeutig 


Hier sei erwähnt, daß die Beziehung (88) auch bei solchen Rechenbeispielen nicht über- 
sehen ı werden darf, bei denen sonst die Reibungskorrektur Pa vernachlässigbar ist; denn die 


F 
e 
» 

u 
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reibungslose Lösung 9; enthält ja nach 6) ein mit dem Faktor In y behaftetes Glied, ad für 
y<0 gemäß (88) in Iny— iz übergeht. Bis auf den hier auftretenden Imaginärteil "vom 
Betrage z auf der negativen Seite von y sind die reibungslosen Lösungen @, und 9, im neutralen 
Fallrein reell. Man bezeichnet diesen Sachverhalt daher auch als „Phasensprung‘‘ der reibungs- 
losen Lösung im Nullpunkt. 


Zur Berechnung der 1. Ableitung von 95, im Anwendungsbereich der Potenzreihe führen 


wir die Hilfsfunktion ein j 


“ 


8,4766 terem ergibt sich die für die Tabulierung wichtige 
Beziehung 


ä | u 
(da) ke -Pum)=—5 Be NEE 
hie “2 R* | 


"® 
Pu... . ee: 
1 | 
Die Integration von (80) liefert dann 
Ö . RR Ayntpl j ; 
Dorn ante amt... 2.22. 00). 
=D aaa (90) 
9=0,1,2 : - - 
Tabelle 5. Hilfsfunktionen für 9,(n) In Tabelle 5, Spalte 2 und 3 sind Wertepaare die- 
Jam Bereich der'Pobenzr oil o (TOjEEzEEEr Funktion für den Bereich 0O<n<6 wieder- 
; a = 5 ER: gegeben. Nennen wir allgemein n, die Grenze 
- = ' : " zwischen dem Anwendungsbereich von (87) und 
ıa|PArm| Fon | am | Prim) | Pin | dr (m) im) _ (90), so wird im Anwendungsbereich der Potenz- 
Ä n N STE PEESERSETT TERRA ER Hreihe,ses 
0,5 | 0,1141 | 0,6236 | 0, 0182 0,1584 IPzı > (2) — Pn)+P(n) . (9), 
10 0,4197 | 1,1359 | 0,1463 | 0,6048 dn dn /n=n, 
s 0, 3266 1 ‚4725 0,4559 1,2646 0 
2,0 IR ‚2352 1 ‚6339 0,9730 2,0478 5 ee 7) v2 
2,5 | 1,5775 | 1,6716 | 1,6797 | 2,8782 wobei der erste Ausdruck rechts nach (87) be- 
an } ‚3323 1, 6484 |- 2, ‚9357 3,7095 ß reehnet wird. . z 
23 | 20 | East gu | Sr Wie man leicht aus (87), (8) und (00) er- 
45 | 22617 | 1.5771 | 5.6455 | 6.1153 kennt, ist der Realteil von (91) eine gerade Funk- 
5.0 | 23614 | 1,5739 | 6.8019 | 6.9029 ; tionin n, während der umx/2 vermehrte Imaginär- 
2 on en Su 1,6898 teil eine ungerade Funktion in n ist. Aus letz- 
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a u gene 
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Es genügt also, die Tabulierung der 1. Ableitung von 93, auf positive 7-Werte zu beschränken 
Tabelle 6 zeigt in Spalte 5 und 6 Zahlenwerte im Bereich 0O<n<8. Bild 7 zeigt den Verlauf 
der beiden Funktionen. 


‘ 


Tabelle 6. Die Reibungskorrektur 9g(n) nebst 1. und 2. Ableitung. 


1... 2, 


d? Qzır 
dn? 
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integriert. 
Für den Potenzreihenbereich führen wir die Hilfsfunktion ein 


on) u; P(&)d 


| 3. 
| = 
| d? YPaıi 
dn? 


1,2879 
1,1677 
0,8594 
0,4878 
0,1777 
—0,0061 
—0,0704 
—0,0635 
—0,0348 
—0,0122 
—0,0021 
0,0000 
—0,0005 
—0,0011 
—0,0008 
—0,0006 
—0,0005 


Falm)n net: Inn+ > (-)" 
n=1 


4. | 5. | 6. ” 
no een dyıai rn ine Sei 
dn dn 
—0,7506 —1,5708 0 —0,9389 
—0,6366 0,9472 0,1437 0,3120 
—0,3309 —0,4349 0,3953 —0,0271 
0,0760 —0,0983 0,8293 —0,1526 
0,4845 0,0631 ‚4709 —0,1548 
0,8268 - 0,1008 2,3021 0,1099 
1,0817 0,0776 3,2826 0,0640 
1,2633 0,0423 4,3711 —0,0342 
1,3990 0,0176 5,5379 —0.0198 
1,5110 0,0063 6,7659 —0,0144 
1,6107 0,0031 8,0467 —0,0122 
1,7052 0,0029 9,3760 —0,0107 
1,7919 0,0028 10,7504 —0,0092 
1,8719 0,0024 12,16 —0,0079 
1,9460 0,0019 13,62 —0,0068 
2,0149 0,0016 15,11 —0,0059 
2,0795 0,0013 16,64 0,0052 


c) Die Funktion 9, 


Die Randbedingung (72) für große n wird offenbar dann erfüllt, wenn man die asyMmpto- 
tische Entwicklung (87) zu 


(In — 2)! 


nl 


1—3n 


(93) 


(94). 


ER RN 8 Darstellung der — um ein lineares Glied 


Er bzw. rin — Funktion 73, (n). Die 


0 Fr Kurve ist eine unsciderd die gesipbelte, eine gerade Funktion in 7 


N, H ; 
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j EN | Gyn+p| q 
Bi .. n+1—D —_ 13 nt Fame. 7 Se). 
j On) EN Bn+p+1) @n+p+ 2" | h | 
p=0,1,2 ; : 


In Tabelle 5, Spalte 4 und 5 sind Wertepaare dieser Funktion für den Bereich 0<n< 6 wieder- 
gegeben. Für den Potenzreihenbereich von @3, (m) ergibt sich dann nach (91) und (94) 


Nu UNE | 
Pan ine = Palm) —nuine+||( ie e me Pa) nm) rem) (6) 


Br eine Darstellung, die frei von dem Parameter e ist. Der Term 9,, (11) — nı : Ine wird dabei nach 
w | (93) berechnet. 


Man erkennt leicht, daß der Realteil von 9,,(n) eine ungerade Funktion ist, also 
insbesondere 
Para). Vals RN 220 
ist, eine für die Tabulierung wichtige Beziehung. Dagegen wird f 2 
e— IT 2 
Pair Su j 
eine gerade Funktion in n. Tabelle 6 zeigt in Spalte 6 und 7 Wertepaare von 9,, für den Bereich 
0<n<sB8. Bild zeigt den Verlauf der beiden Funktionen. 
9.1 wurde bereits von Schlichting [8] tabuliert. Wir hielten jedoch eine genauere 
Tabulierung für erforderlich, um dem heutigen Stand der Stabilitätstheorie gerecht zu werden. 
Zudem haben die von uns benützten Formeln den Vorzug größerer Übersichtlichkeit. ; 
d)\ Bemerkungen zu MeksynsKritik De N, Bi 


Meksyn [4] löst die Störungsgleichung (1) für den Fall des ebenen Kanals bei aus- 
gebildetem Parabelprofil. Er kommt dabei seiner Ansicht nach zu einer anderen Übergangs- | 
substitution als Tollmien und folgert merkwürdigerweise daraus sofort, daß Tollmiens 
Übergangssubstitution falsch sei. Wir kamen vorstehend in Gleichung (88) auch zu der Toll- F 
mienschen Übergangssubstitution. Wir folgern hieraus jedoch nicht, daß Meksyns 
Übergangssubstitution falsch sei, sondern haben uns statt dessen der Mühe unterzogen, Mek- 
syns Rechnungen mit den unseren zu vergleichen. Dabei ergab sich der folgende Sachverhalt: 

Nach einem bekannten Satz darf man der Lösung einer inhomogenen linearen Differential- 
gleichung stets jede Lösung der entsprechenden homogenen Differentialgleichung super- 


ponieren. Daher ist in unserem Falle auch Ser | | a 

Bor = Paıt k*Ygı WR EL 274 Ver Tee | 
stets eine Lösung von (71). Ersichtlich wird auch die Randbedingung (72) für „>0 durch die 
Superposition nicht gestört, weil ja 9,5, hier sehr stark abklingt, wie (56) zeigt. ‘Wählt man 
speziell ’ x In r 


} Be ee En 
N [3 — 1) +7 V3 KO] er. 
so zeigt die Rechnung, daß der Ansatz (98) genauestens das Konjugiert-Komplexe dr von 
Meksyn[4]gefundenen und p. 397, Tabelle 1 tabulierten Funktion 7 M+n ergibt. Mek- 
syn war von dem zu (2) konjugiert komplexen Ansatz ausgegangen) a 

"Meksyn hat also lediglich eine andere: Aufteilung des Fundament 
tikulärlösungen zu (1) erhalten. Uns erscheint allerdings die Tollmien 
wir vorstehend gefolgt sind, insofern praktischer, als bei ihr ,, nicht von 
im Negativen sehr große Funktion 95, verdeckt ist, sondern gerennt vo 


k 


5. Die innere Reibungsschicht bei angefachten und gedäm 
Während bisher im wesentlichen der Spezialfall der ‚neutrale 
wurde, der zur Errechnung der Indifferenzkurve führt, sollen in diesem ] 
Einzelheiten über den allgemeinen Fall c; =0 ausgeführt werden, d« 

Vorzeichen von c; zeitlich angefachte oder gedämpfte Störschwin; 

dabei verschieden vor, je nachdem c; endlich groß oder so klein 

berücksichtigt zu werden brauchen. FRE 
2) Nach Fertigstellung dieser Arbeit erhielte N wi 

syns [5], in der er seine Kritik an, Toll mie 


Ba. ” Nr.1 an En los Holst e im, Be die BerbungescHioEr bei Störungen laminarer Strömungen 45 | 


a) Der Fall endlich großer s,-Werte 


In diesem Fall hat die reibungslose Differentialgleichung (5) bei U = c, keine Singularität. 
Statt (6) erhält man zwei Partikulärlösungen 


Y 
Da: yr, wert (mas Dit 0... 2 AERLUCN: 
v U 1 n=uU 


Unter Verwendung der Reihenentwicklung (4) für die Geschwindigkeitsverteilung erhält man 
dann für 9 stetigen Anschluß an die neutrale Lösung @, in (6), wenn man 


€ 


TE und EA N SE ER NL 
setzt. Als weitere Koeffizienten seien noch angegeben 
U" — 1026; @U’+U” 
en zu 2 
a, 3U R a”, 67 (102). 
Die Koeffizienten von 9% sind dann nach (100) ebenfalls bestimmbar: 
.2U, 1 U, 
bi ==; Ge: =. — 5(e+:%) a 7371033 
usw. 


Wie man sieht, ist hier im Gegensatz zu 9% kein stetiger Übergang c6;— 0 zu der neutralen 
Lösung 9, möglich. Hieran ändert sich auch nichts, wenn man die Integrationsgrenzen des 
Integrals in (100) anders wählt. Wir werden daher, wie schon erwähnt, anschließend den Fall 
kleiner e;-Werte gesondert entwickeln. Hier seien nur noch einige Unterschiede zum neutralen 
Fall angemerkt: 

Bei der Lösung (100) wird natürlich eine Reibungskorrektur 93, für die Umgebung von 
'y = 0 überflüssig; denn die Einführung einer solchen (vgl. Abschnitt 4) wurde ja erst notwendig 
durch die im neutralen Falle hier vorhandene Singularität. 

Die einzige imaginäre Größe in der reibungslosen -Differentialgleichung (5) ist die Kon- 
stante © c,. Daher wird jede Lösung von (5) mit 6; > 0 konjugiert komplex zu der entsprechen- 
den Lösung von c;<0. Dieses einfache Verhalten der angefachten und gedämpften Teillösungen 
zueinander überträgt sich nicht auf die Reibungslösungen 9, und Y,, deren Differentialgleichung 
(18) ja schon im neutralen Fall komplexe Koeffizienten enthält. 

Was an Stelle des ‚„Phasensprungs‘‘ der Teillösung ®, des neutralen Falles für den Fall 
6; #0 tritt, läßt sich nicht allgemein formulieren; denn jetzt ist 9X eine für positive wie ne-. 
gative Werte von y komplexe Funktion, wie aus (103) entnommen werden kann. 


b) Der Fallkleiner ,-Werte 


c; werde nunmehr so klein vorausgesetzt, daß bei den reibungslosen Lösungen nur in c; 
lineare Glieder berücksichtigt zu werden brauchen. Führt man formal mit Schlichting 
statt y die komplexe Veränderliche 


Ve ee Rt... (104) 


ein, so kann man nun die Entwicklung (4) für die Geschwindigkeitsverteilung in eine solche 
nach Potenzen von y* umbilden. Setzt man also , 


Ge ns... (105) 


und wählt g so klein, daß nur in g lineare Glieder berücksichtigt zu werden brauchen, so ergibt 
‚sich sofort die Koeffizientenbeziehung 

BU’); ET: 

=; Ei Dre = VER UND .g (gültig für n20) .. . . (106). 
Legt man nun g fest durch AU 


a ee Cr, 


E; n; nach (3) und (106) fürn-0 5 
Di Dalaran = +, =6 Bus, /. Seine > 1. . (108), 


= nibungsoe a nizleich "wieder formal die Gestalt (5) erhält, wenn man 
} er 


Ex 
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IE er N 
dort die Differentiationen nach y* statt nach y angeschrieben denkt. Die partikulären Lö- 


ge: sungen dieser Differentialgleichung lauten also analog (6) ei Ha! 
F- 
oo oo ’ TA \ . i > 2. 
= aD y*, BD by ter -In y* Te (109), > 
HU Been=V r 2 
4 


wobei die Koeffizienten mit denen in (6) identisch sind. 


y% benötigt analog dem neutralen Fall eine Reibungskorrektur 9}, für die Umgebung von 
y* =0. Mit den Einführungen CR 


a ie 20 9 


ı 


Br * U 29 | i 
f n PB: { Y ae nn >. (110) En 
; erhält man dann nach (104) # 
3 | nr ee 1 X 
Setzt man noch analog (70) ; 

x TH 

Ze a ee ‘a 


Tabelle . Der Imaginärteil der 2. Ableitung von 9# in Polarkoordinaten 


ER 90° |—67,5°| —A5° |—22,5°| 0° | -+22,5° |  +45° 467,52. |. 4.808. 908 

ba | Ä 

0 1,2879 | 1,2879 | 1,2879 | 1,2879| 1,2879 | 1,2879 1,2879 1,2879 1,2879 | 1,2879 

0,5 | 0,921 | 0,932 | 0,970 | 1,044 | 1,1677 | 1,355 1,597 1,820 1,894 | 1,915 

1,0 | 0,693 | 0,698 | 0,715 | 0,755 | 0,8594 | 1,123 1,705 2,575 2,962 | 3,061 2 0 
1,5 | 0,543 | 0,541 | 0,529 | 0,510 | 0,4878| 0,572 1,277 3,485 . 4,340 | 5.261 N 
1,75 | 0,486 | 0,483 | 0,460 | 0,406 | 0,3182 | 0,224 0,726 3,882 6,307. 12:7.1922 Ser 
2,0 | 0,443 | 0,430 | 0,403 | 0,331 | 0,1777 | —0,135 | — 0,107 4,067 8,297 | 9,830 „4 
2,25 | 0,405 | 0,392 | 0,353 | 0,270 | 0,0689 | —0,451 | — 1,228 3,769 10,995 | 13,854 F 
2,5 | 0,369 | 0,359 | 0,313 | 0,221 | —0,0061 | —0,684 |— 2,577 2,541 14,620, }.19,926 7209 
3,0 | 0,317 | 0,300 | 0,256 | 0,167 | —0,0704 | —0,805 |— 5,251 |—. 6,122 25,620 | 43,741 ER 
3,5 | 0,276 | 0,260 | 0,214 | 0,120 | —0,0635 | —0,515 |— 5,572 |— 33,839 41,800 | - 

4,0 0,233 | 0,183 | 0,096 | —0,0348 | 0,000 |— 0,111 |—-102,766 51,857. 

4,5 0,201 | 0,165 | 0,082 | —0,0122 | 0,239 11,271 |—234,089 | —17,069 

5,0 0,184 | 0,145 | 0,078 | —0,0021 | 0,134 18,860 |—367,326 | * CR 
5,5 0,168 | 0,134 | 0,081 | —0,0000 | —0,131 3,406 |— 56,783 |... Ba 
6,0 0,157 | 0,139 | 0,105 | —0,0005 | —0,234 |—40,584 | 1024,514 ae 2 


‚so erhält man die zu (71) m 2 
ganz analoge Differentia- 
:gleichung er 7 RE 


iR Min —in) für >00 


z Mar } ı% 
Bild 9. Darstellung des Imaginärteils der 2. Ableitung von 


x ei ’ 
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bei festgehaltenem Betrag |7*| und veränderlichem Argument # tabelliert, weil dann die für 
den neutralen Fall durchgeführten Rechnungen mitverwertet werden konnten. Auf graphischem 
Wege wurden dann mit Hilfe dieser Zahlenwerte Kurven über n bei y = const ermittelt. Diese 
sind in Bild 9 und 10 wie- 
dergegeben. 

Die Integrale der in 
Bild 9 und 10 darge- 
stellten Funktionen über 
den ganzen Bereich von 
n= tr ® . bis er 
bedeuten bekanntlich im 
neutralen Fall den ‚Pha- 
sensprung“ —rz. Da es 
sich hier um gerade Funk- 
tionen handelt, genügt es, 
die Integrationen von +% 
bis O0 durchzuführen. Der 
Integrationsweg ist dann 
nach Bild 11 je nach dem 
Vorzeichen des jeweils 
vorgegebenen y=y,durch 
a bzw. b gegeben. Statt 
dessen kann auch z.B. —n 
längs der Wege q,, Q, 43 
bzw. b,, b,, b, integriert 
werden. «a, und 5, sollen voraussetzungsgemäß im Unendlichen liegen, so daß diese Teile des 
Integrationsweges wegen des Abklingens der Funktion nach Null keinen Beitrag mehr liefern. 
Der B:itrag von a, bzw. b, wurde in (92) zu —/2, eben der Hälfte des Phasensprunges, errech- 
net. Der Beitrag längs der imaginären Achsenteile a, bzw. b, ergibt sich aus (91) und (90), 
wenn man dort +%-y, bzw. —iy, statt n einführt. Man sieht dann sofort, daß (90) nur 


Bild 10. Darstellung des Imaginärteils der 2. Ableitung von gi mn —iy) für y<0 


Bild 11. Zur Bestimmung des Integrationsweges für a bei Anfachung (a) oder Dämpfung (b) 


reeller Werte fähig wird, so daß auf a, bzw. b, kein Beitrag zum Imaginärteil von dg,,/dn ge- 
liefert wird. Wir erhalten also das wichtige Ergebnis, daß bei schwacher Anfachung oder 
Dämpfung das logarithmische Glied in 9} wie im neutralen Fall durch 
In y* = In (y—i:g) =In|y*| —i- Vaoeoee..... (115) 
gegeben ist. ? 
Auch hier läßt sich wie in (23) eine für das ganze Integrationsgebiet gültige Funktion ®} 


zusammenstellen: 
n oo * = 
Di= I bu: yike or, (2) ED er RER, 30) 
‚ vn] 


n=U 


Wie die für endliche c;-Werte vorher abgeleitete Entwicklung p} in (100) bei abnehmendem 
6; an (116) anschließt, kann allgemein nicht eingesehen werden, sondern muß an Hand von 
„Beispielen bei gegebener Geschwindigkeitsverteilung untersucht werden. 
Die asymptotische Entwicklung in der Form (75) von p}, galt, wie erinnerlich, nur mit der 
Einschränkung, daß |n*| genügend groß ist und das Argument von n* die Sektorengrenzen (76) 


- 


nicht erreicht. In Bild 12 sind diese Einschränkungen durch Schraffierung des — als nicht 
abgeschlossen anzusehenden — Gültigkeitsbereichs von 9,, angedeutet. Die Vorschrift ds 
Abklingens der 2. Ableitung von p#, wie 1/n*(vgl. Gleichung (73)) gilt also nach Bild 12 z.B. 

- nicht längs der negativen y-Achse. Dies kann aber auch nicht verlangt werden; denn wenn 
(73) für jedes Argument von n7* gelten würde, wäre d?9%,/dy** im Unendlichen regulär. Daes 
nun im Endlichen aber auch überall regulär ist, wie die Potenzreihenentwicklung (80) zeigt, 


Z 


Bild 12. Veranschaulichung des gemeinsamen Gültigkeitsbereichs (schraffiert) der 
Hankelschen Funktionen 1. und 2. Art von der Ordnung 1/3 bei asymptotischer 
Darstellung in der komplexen n*-Ebene , 


so wäre es dann nach einem bekannten Satz der Funktionentheorie eine Konstante, was offenbar 
nicht möglich ist. Daher kann d?p%,/dn*? im unendlich fernen Punkt der komplexen n*-Ebene 
nicht regulär sein. In der Tat wächst z.B. d?9,,,/dn?, wie sich auch in Bild 10 andeutet, längs 
der negativen y-Achse über alle Grenzen, während es längs Linien y = const der Vorschrift (73) 2 
u folgend nach Null abklingt. 
c) Stellungnahme zu Lins Aussagen über die innere 
| Reibungsschicht = 
Wie wir früher näher ausgeführt haben, hat die Partikulärlösung 9, (y) der reibungslosen De 
Differentialgleichung (5) im neutralen Fall eine Singularität im Nullpunkt U=c. Dadiesse Bi 
Singularität in der wirklichen — reibenden — Strömung nicht vorhanden ist, ist man zu dem 
Schluß gezwungen, daß an dieser Stelle nicht reibungslos gerechnet werden darf. Es wurde 
daher mittels der Reibungskorrektur p,,} eine verbesserte Lösung ®, (Gleichung (23)) geschaffen, 
die im ganzen Strömungsgebiet, insbesondere also auch bei y = 0, regulär ist. Die Reibung ist 
also in der Tat im neutralen Fall im Punkte y = 0 eine wesentliche Größe, und man ist dher 
berechtigt, diesen Sachverhalt dadurch auszudrücken, daß man von der Existenz einer „inneren 
Reibungsschicht“ um den Nullpunkt herum spricht, die in dem gegebenen Reibungsschicht- 
profil U(y) eingebettet ist. In den ersten Teilen von Abschnitt 5 wurde dargelegt, daß die 
Singularität der reibungslosen Differentialgleichung — und damit die innere Reibungsschicht — 
bei endlich großen Anfachungen oder Dämpfungen verschwindet, daß sie dagegen bei kleinen 
Anfachungen und Dämpfungen noch vorhanden ist. AR } f: 


hier verwendeten Hankelschen Funktionen 1. und 2. Art betonten, sind diese Entwicklungen _ 
kein notwendiges Erfordernis für die Durchführung der Stabilitätsrechnungen. Sie haben nur. 
. . em 


Potenzreihenentwicklungen von der Art (26). 
Lin [3] hält irrtümlicherweise offenbar die Verwen 


esnunL i n, die Gültigkeitsgrenzen der asymptotischen Entwicklungen der Jar 
nen physikalisch zu interpretieren. Wir halten es für nötig, hierauf ets 
Wir halten uns an Bild 12. Denken wir uns 7 und y bzw. durch yund 


hi 


(110)), so schrumpft wegen der Kleinheit von e der 
4 I 


Radius des Kreisboge 
EI. We ei y RR 
- 


Pi A 
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grenze der asymptotischen Entwicklungen der Hankelschen Funktionen andeutet, etwa auf 
Null zusammen, während die Sektorengrenzen für den Gültigkeitsbereich des Arguments er- 
halten bleiben. In dem so geänderten Bilde durchstößt dann eine Parallele zur reellen Achse 
in der oberen Halbebene (= Integrationsweg bei Dämp“ung) zweimal die Sektorengrenze, die 
reelle Achse selbst (= Integrationsweg im neutralen Fall) berührt die Sektorengrenze im Null- 
punkt und eine Parallele zur reellen Achse unterhalb derselben (= Integrationsweg bei An- 
fachung) verläuft ganz im Gültigkeitsbereich der asymptotischen Entwicklung. Jedes Er- 
reichen der Gültigkeitsgrenze beim Durchlaufen des jeweiligen Integrationsweges wird nun von 
Lin mit dem Vorhandensein einer inneren Reibungsschicht identifiziert. Nach Lin gibt es 
also im Falle der Dämpfung zwei innere Reibungsschichten, die symmetrisch zum Nullpunkt 

- liegen, im neutralen Falle eine, im Falle der Anfachung gar keine. Diese Neuheit kann man 
p- 133 der Linschen Arbeit in besonders hervorgehobenen Lettern nachlesen. Diese vermeint- 
lichen Tatsachen werden dann noch im Linschen Text durch Gleichgewichtsbetrachtungen 
zwischen Wirbeltransportgrößen der Haupt- und der Störströmung, der Diffusion von Rotation 
durch Reibung u. dgl. in den drei Fällen ‚‚erhärtet“. Tatsächlich zeigen die Bilder 9 und 10 
ganz deutlich, daß sowohl bei neutralen wie auch bei schwach angefachten oder gedämpften 
Störungen eine einzige innere Reibungsschicht existiert. 

Wenn auch für den Leser, der nicht nur Lin [3]sondern auch die Originalveröffentlichungen 
der Schöpfer der Stabilitätstheorie etwas kennt, keine Gefahr besteht, etwa auf den vorstehend 
geschilderten Fehlschlü sen weiterzubauen, so besteht diese doch um so mehr für alle anderen, 
und im besonderen Maße auch für mehr mathematisch als physikalisch vorgebildete Forscher. 
So hat neuerdings eine im übrigen sehr verdienstvolle Arbeit von W. Wasow [17], die im 
wesentlichen mathematisch orientiert ist, in der Tat die Linschen Fehlüberlegungen nochmals 
abgedruckt. Wir hielten es daher für unsere Pflicht, die Kritik an dieser Stelle etwas aus- 
führlicher anzusetzen. 
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Spannungsfunktionen für die von Kugel- und Kegelflächen 
begrenzten Körper und Kuppelproblem 
Von A. Timpe in Berlin-Charlottenburg i we 


In vorliegendem Aufsatz wird der auf eine erzeugende biharmonische Funktion gestützte Lösungs- 
$ schlüssel, den Verf. in Bd. 28 (1948) dieser Zeitschrift für die achsensymmetrische Deformation von Um- 
4 drehungskörpern in räumlichen Polarkoordinaten entwickelt hat zur Gewinnung von typischen Lösungen 2 

für die von Kugel- und Kegelflächen begrenzten Körper nutzbar gemacht, insbesondere einer Gruppe von 
„„Brückenlösungen‘‘, die für das Kuppelproblem (verteilte Belastung oder konzentrierte Last in der Kuppe) 
: in Frage kommen. Anschließend Kegelprobleme. 
\ In v. 28 (1948) of this journal, the author has developed a method of computing the awially-symmetrical 
af deformation of solids of rotation in polar coordinates, based on a ‚‚generating‘‘ biharmonic function. The 
method is now used to get Iypical solutions for bodies bounded by spherical and conic surfaces, especially a 
group of „‚bridge-solutions‘‘ that are impgglant for the cupola problem (distributed load or load at the top). 
2 Subsequently, cone problems are treated. - 
5 Dans le vol. 28 (1948) de ce journal l’auteur a developpe une methode pour determiner en coordonnees 
polaires la deformation symötrique par rapport & l’awe de solides de rotation, se basant sur une fonction 
h , «generatricey biharmonique. Ici il empioie cette methode & atteindre ü des solutions Iypiques pour solides, 
dont les surfaces sont spheriques ou coniques, et sp£cialement & un groupe de «solutions de pont), qui sont 
importantes pour le probl&me de la coupole (charge distribuee ou au cime). Alors des problömes de cone sont 
traites. 


Br. 28 (1948) 3Toro »xypHana aBTOPoM ÖBIT PaspaboTaH MeTon ONpeneleHnnd OCecHM- 
METPu4YHEIX mehopmannä Te BpAmmeHHA B HOJAPHBIX KOOPAHHATAX, OCHOBBIBAIMINMÜCH H& 
IIPuMmeHeoun „oÖpasylolmux‘‘ ÖHrapMOHHYeCKHX PYHKIHÜ, TOT METON HCHOABSYETCH BECK 
AAA NONYYeHHA TUNMYHBIX PeIIeHHÄ MIA TEIL, OTPAHHYEHHEIX CDePmyeckuMH HH KOHHYEC- 
KUMH HOBePXHOCTAMH, HU B OCOÖCHHOCTH AA TMOAyYeHHA HEKOTOPOH TPyImkI „‚UePeKRPEI- 
BAIOIIMX‘‘ PeimeHHÜ, HMEIWIIHK BHAYEHHE MIA PacyeTa CBONOB (C PacıpenereHHuoH Harpy3RKoK 
HH C COCPeNOTOYeHHOH HATPY3KOH B BEPIIHHe). 3ATEM PACCMATPHBAMTCH CIIYYaH KOHH- 
yJecKUX Teil. . { 


In dem Aufsatz ‚‚Torsionsfreie achsensymmetrische Deformation von Umdrehungskörpern 
und ihre Inversion‘‘t) hat Verfasser einen Lösungsschlüssel für das achsensymmetrische Elasti- 
kr: zitätsproblem entwickelt, der demjenigen analog ist, mit dem er früher?) im Anschluß an I. H. 
R Mitchell?) einegrößere Zahlvon Beispielenfür das ebene Problem behandelthat und der durch 
die späteren Arbeiten von A. und L. Föppl *) für manche weitere Probleme nutzbar gemacht 
worden ist. Kennzeichnend für diesen neuen, auf räumliche Polarkoordinaten gestützten Lö- 
sungsschlüssel für die achsensymmetrische Deformation ist die besondere Eignung für Um- 
drehungskörper, die von Kugel- und Kegelflächen begrenzt sind. In der Praxis spielen derartige 
achsensymmetrisch beanspruchte Körper eine außerordentliche Rolle. Auch die Materialprüf- 
versuche und die der Bestimmung von Elastizitätskonstanten dienenden Belastungsproben 
werden vorzugsweise an solchen Versuchskörpern angestellt. Die unserem Ansatz als Grund- 
lage dienende erzeugende Funktion F, die als Funktion des Radiusvektors r und der Poldistanz x 
erscheint und die der Air yschen Funktion des ebenen Spannungsproblems analog ist, genügt 


h 3 5. 2 
N (in 2) 
der Bipotentialgleichung F7?V7?F= 0, wo jetzt V?F = ur ae) Be: 0 Ne 2. 
; r or: ‚2sin x 0% 


Aus ihr leiten sich die Verschiebungskomponenten eines Punktes P in Richtung wachsender 
r- bzw. x-Werte ab durch die Formeln: 


. 


ER NR AA N Ne h — (sing 29); 5 
or ı m = 2 Ber 0% 2 
oF m 1 ( = 
U NEN TEE N... EE r-—|* 
% X m—2 r 9 °x x 


der durch P gehenden Achsenebene und dem Nullmeridian, verschwindet. 
Die zugehörigen Formänderungskomponenten sind: 


Die Verschiebungskomponente us im Sinne wachsenden Winkels ©, wo © der Winkel zwischen 


ou, 1 ou u il Ur 
Dehnungen: ge, m —. Ze sur. Sn 
5 u Ta u 


2) BB ur (5. Abna 


Winkelverzerrungen: Yey =Yye, =0; yy=r- & ( Er 
* we 


or\r 
v 1 1) } 
Kubische Dehnung: A1= Psıny : E (r? Uy sin 4) B= . rg sin | 


*) Z. angew. Math. Mech. 28 (1948), S. 161-165. 3) Lond Mat > Soc. Pr 
2) Z. Math. Phys. 52 (1905). = Dre Ya LI 
A u ee 
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‚ Zugehörige Spannungskomponenten ; 
mE A mE A mE A 
em an 
0 


"mtl m + 1\m — j 
mE (3). 


> Tex Bu Y(m - (m 1) ) "Iry 
Wiei in dem Aufsatz von 1948 betont, ist dann für die Behandlung von Körpern, die von Kugel- 


und Kegelflächen begrenzt sind, zweckmäßiger Ausgangspunkt der auf zwei harmonischen Funk- 
tionen 9, @* gegründete Ansatz 


Tyo = Tgr 


re De. 0 u re DRAN 

der allgemein die biharmonischen, d.h. die der Bipotentialgleichung genügenden Funktionen F 

erfaßt. Gelegentlich empfiehlt sich in diesem Sinne auch der ebenfalls auf 2 harmonische Funk- 

tionen gegründete Ansatz 
F=z2:9* +9 =o*.rcosx-+p. 

Dabei kommen für Reihenentwicklungen als geeignete Elementarlösungen der Potentialgleichung 


in Frage: 
1 
ı=r”- P, (cos z); ve -Pn(cosy); | 
’ (9) 
V,=r.Qn(cosY%); Ys= EN -Qn (cos x) | 
4 


wo die P,„ (cos y) die für 0 Sy Sr durchweg regulären Kugelflächenfunktionen erster Art, die 

Q„ (cos x) Kugelflächenfunktionen 2. Art, die für y = 0 und y = divergieren, im übrigen regu- 
lärsind. Dasich aber, wies. Zt. betont, bei den Gliedern zweiten und niederen Ranges Überschnei- 
dungen ergeben, erfordern diese eine besondere Untersuchung. Dieser Aufgabe dient die vor- 
‚liegende Arbeit). Es zeigt sich, daß die Ergebnisse, die ich als ‚„‚Brückenlösungen‘‘ bezeichnen 
möchte (von denen einige bekannt sind), für manche praktischen Belastungsfälle in Frage j 


kommen. 
I. Brückenlösungen des Bipotentialproblems 
a 1 2 
(Aı) FR Ban): U = 7; uy =d. 
4 2 2 | 
EORT Er a a Pre BT Br SAN Yy=: 
AmE 1 2mE 1 
y=— { = -=0 Tym0. 


mir Kar 
Auf die Kugelfläche vom Radius r wirkt zentrosymmetrisch die Gesamtspannkraft 


EmE. 1 
— 16 .—, 
ri r 
NT: mer Bar = 6, Her; uUy = 0. 
(A,) rd. en, BNE:9: A=6: —) 
ee N v AT ’ 0 — 5, Br 3 Yrx = A 
2 mE 
| ‚ 9, —0y a - BE. ry=0 


Die or ahion von A, und A,liefert die bekannte Lösung für das Problem der Hohlkugel, die 
aufr =r, durch den Druck p,, auf r = r, durch den Druck p, beansprucht wird. 
En 2 pn mn +! nn(dı— 29) 1. 0. 


= 4 Ar 
mE n—rn 2 mE nn a 


Ft) — e er * 
(Bi) 2. — PF=0; 200%, u = —. 
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EEE 


—n 2cosyx 4m —1) cosy S3m—4 singy 
m Pe Zyaa  \ Be. 2 12, ee 
(B,) ha a m— 2 2’ 0% m—2 r 
4A'm —1)  cos% MEENCO" X 2 ! 2<inX 
= BEL; . 2a en Pr = &9;5 A= za 0a u 72 
2m(2m—1)-E cosyx mE cosy mE siny 
= ht RER, OT 3 = 08; 1 eg Dr: 
(m +1) m — 2) r m+1l r m+1l r 


Dieser Spannungszustand wird nach Boussinesq ‚Applications des potentiels‘ als der 
„erste Typ einfacher Lösungen‘ bezeichnet. 


er 1 m — 4 3m — R 
(B) F=— rcosyg; V?F=d5rcosx; Un zrrcosy; u = I r®sin x. 
_ - Be Br, ee 
m—4 ER 2 (2m — 3) Nr 2(m+1) 5 
& =2: —— 1°C JEy= ——— = $& rsi 
ö m—2 %°% m— 2 le AT TE % 
2m dm mE Mm 2 
= — 7C0SY; = —-— rC008SX4; 098 = ——1C005Y; Ty—= - rsiny. 
er RR, I RE a N A a) k 


Man bestätigt leicht, daß die Potentialgleichung 
une (e.P) 1 R) | : > 
2 et N ee 
BETH Br up öx un °x N 
1 
von den Funktionen 9, =Intg z und 9, = 2 In tg = befriedigt wird. Demgemäß sind 
außer ihnen gemäß (4) auch ro, = r?In 5 und ro, =r:Intg 3 als Spannungsfunk- 


tionen heranzuziehen. Dadurch erhalten wir die Spannungsverteilungen der Gruppe D,, D,, 
DD... D;: 


D ER A. pen _ , Br , RANGE ARRE, 
(D,) F=r?2-Intg 9: v?F 6intgs; u =ArIntgn; SE 
Earl 2cos SX 2cosx X 
—2lntg<; — . CASE k; 
&r ntg ER A an eg tg 5° &8 a +21ntg 2 
] 2 
A= 6Intg 2° Yry sin? y 
2m 2 2mE 2mE cosyx, 
Kr -] t a =. ee RES I! 
HET N NN In tg En, 177 
2 mE X 2 mE cos, mE 1 
(ee = > t GESRLCHE b %y>= "RR . 
m—2 2. m+1 su2y’ m+1 siny 
D,) F=-. ft, Pr —0: or RR En © 
(D,) 5 nig 5: FARO: w— Intg 5; Uy = —a en 
22 X 2 cocx 2 If I 4 
&r Int Zee: = —« . eu == Has, N 
85; &% 73 Sina 7a Intg 5; 8e I he Inte; 
8 1 
A ==); u a N 
IT 5 sing 
N mE. 4 Sy mE 2[rosx 
ee 
mE | cos x mE 4 i 
OA = — I — Sie. une 0 Finden , 8) menu 
I m+tir may un ig 2); Irx m+1 r sing 
D I —a 2m _.N. RR Pd 1 
(D3) 3m —4 Inte; BE ; W=0; u, m 
l cosy Fr 
Er Em 7 Sr, 
& X sing’ @ 2; d=0; yry Fat.) 
’ mE 1 cosx mE 1 
ol; Te z——_ . 4 En { 


c —— NO 
% m+l rn sin? x’ ® 


‚ke na pi 2 
‘8 (H}) Pr #2. 4. (oosginigk +1); v2F =0; RE 


‘ 1 \ 
a (eis —sin Intg 2: 
A — 
rt sin2y 


6 | 3 
a losz migitt); 4-4 + lit eosznig 2); 


3 1 
8 = 7:008%' In tg, ur (2 —cigy); 7 A=0; 


6 ; X 
Bee ca etgxy —sin x ntgn 2 


6mE 
0%, = — 


L Bee 
4 (eosz In1g% +1); 


m +1 

$ rn u tm ale net +1); 
3 0 4: (30057- IntgZ +2 — ctg? ): a 
F nee A g 
m+1 rt? 2 < 
(B,) F = cosy:Intg =; ver = —z(eosg-mig&+1); £ & 
% = — un, 2 (cosy: -Intg— N). un 2 (eigy-sing-intg2). E 
„Er. 4 (cos z Inig +1); | ! 

7, ee: a na a lesr Inte, +1); 

ae ara wanez.): 


3 (cos2- Inte, +1); = 3 (eig2—sinz- -Intg “a 


2m(2m—1)E 1 i 74 ; 
= Ka m Bor a6: Inte, u: 


mE 1 z( ) nA 1 |: 
E23 le AS ng ni iR m—2 rsın?y.’ 

mE 1 mE 1 ß y\ 
er m+1 0,0 (nern ae 2), en aletez-sinzintg?) 


Dre cos x: nt; ni; u N (oszinte£+1); 


Beifügung des Faktors r? und anschließende Weglassung 
alfunktion D, durch Beifügung des Faktors 2= r -cos x. 

E en des Faktors 2 = r - cos He 
0 zusammen. Doch we (wie hier) A=0 schon 


“ ” „.J 
RN — (sin ny: 7) - 
7 Rs 
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Ama, 1 . 
ir - I. — en; A=0N): =0. 
a 2 m—2 rsin:y z % Be 


0: Am—l) mE, 1 An 
zT) De mA rsin?y ar 2°. A 


ee 1 X 
(H,) F=trlosynig2 +1). 
Diese Funktion genügt derDifferentialgleichung F?2V?F = 0; denn V?F=5r (cos; In tg 3 2% ı) $ 


2 Dr: 
Nach dem bei H, erhaltenen Ergebnis ist aber: 7? (r -cosy-Intg e Sep andererseits ist 
2 ui 
bekanntlich F’r= = 


Für diese Spannungsfunktion ergibt sich: 


ww a »[eosz: in tg 3 +1); uUy= ae : (etigy —sinzin tg = b 
N. "—ör(cosz: -Intg 4 2 +41); 

4. ni lu ts Ar Enz 

en Yry= 2 (aeg) 

GH, = Sue ‚(eosz: In.tg > 241); 

rl mg lH) 

0 "5 rleoszintg &41) A m 

pe u; rleigz—sinz-Intg a 


= In[r(1 +cosy)] =Inr+In (1 +cosy), eine der Potentialgleichung genügende 
Funktion: 


0) 
or? nsiny dy Fa ahUr 


EN 
/p 
2, 
Be 0 
sing dy nr 2 «sin 2 
Als Spannungsfunktion angesetzt, liefert F: 


a ee er 
or m—2 rsing 9y\l+cosx 
2 1 g Bi 
nn m 7 (2sin? 2\- 2 2,1 Zn, 
r m—2 7-siny.dy 2 r,.m—2ır m —2 7 


Br, m ‚1.0 > 
4 ray m—2 r dr 1+cosy) r 1-+cosx 
„m, 1, Zn Pe RE 
m—2 r L+esy mar 1 cos 
m — 1 1 sinx 
Für F= 22 BE EEE . 1 ZA 
ür BT "Tin r+In (1-4 cos y)Tergibt sich also u 73 72° Aue TFcosx und 
1 1 1 1 RARBSInEN 
gg g= on ze 
Y 272 cos? % 5 2,72. cos? £ 2 EI 
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2: De su Dr 
mer: x mi 1+cosx m’ e em+1 1+coy nr’ 
X min 1 Be, 
Dieser PanungszustZie wird Bach Boussinesgq als „zweiter Typ einfacher Lösungen“ 
bezeichnet. 


eG) Z=72.In[r( a] =? In Bar r2- In ( x); gemäß (4) aus T, Beleen 
VP?F=6-In[r(l-+cosy)]+ 4. 
W=6r: mer6r: TR Dr er 


warn. 1 B) uk) 
Bu Ir. 
m—2 ar er 1-+co0sy Be _ az N 
); 
< sinyx m 1 N sing ) 1 sinx 
ee el EEE 85 I EEE De Kae, 
= n 1+cosy m—2 r 9r ‘ 1 +cosyx K 1-+cosy 
3m sinx a>m—4 sinx 
’ + —:rf: — m... 
m— 2 1+cosy m—2 Be 
3m—8 : Im—4 M— 
= 2:1 1 —;',= 2. 1 ; 
| 2, n[r( (bei 008 2], Be a EZ —+2-In[r( +c0s%)] + — =; 
Bi: 5m—4 
za las Se he +2-In[r(1+cos n)+— Be 


m— 2 1-+cosy 
2 m+l). sinX 


A=6-In[r(1-+cosy)] +10; ioce BE DE 
mE 


8 
0, — (m+1) (m 2) 2. -In[r(l cos oe: 


Bi 
= en, In[r (1-+ cos Allan 10} + mi -In rd +cosy)]+ 


er BR Da 
m—2 1I+cosy m—2)’ 


5, (6 In[r(1+.cosx)] + 10 Di 


= mE a 
u; —— 1 
Re 06= mr pres 5516 -In[r( (1 Heosyl+ 10} vo 2ln[r(l+cosg)]+ 
b \ , ImMm— ae I, cos X en ER 
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.& F= =— ee In sing, eine der Potentialgleichung genügende Funktion. 
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0 4 (m —1) 1 3 Mm AN 
armen Tumor ma 
4m—1) ctgy 3m 1 2. 4(m —1) ctgx 
= — rn ER 
m—2 r m—2 19 r m— 2 r 
: 1 4(m —1 
Ba 2mE SE Ba mE 23m + udn |: 
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mE 2(m—1)mE ctgx 


u : [2 —3m—Am—1)dtgx; ya 


00 7 (m+1)(m— 2)  m+i)m—2) 


— 2 s r e x & 
(U;) = _ er Infrsiny], eine der Potentialgleichung genügende Funktion. 
1 1 
Ur Be) Un EX 
1 ctg?y ee 1 2ctgy 
Du ro 
mE 1 mE ctg?y mE ET 
Bu = a 10 Fe sine 
m+1 r m+1l r m+1 r?sın“y 
mE ctgy 
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(U,) F=r?-In[rsiny], gemäß (4) aus U, abgeleitet. 


ee 


u = 2r-Infesing)+ 5: r; u RT oigr. 

6 = 2-In[rsing)+ Rn; 24 = 2 Inlrsin ne 
co = 2-Inlrsing]— ZI og +5; 

A=6-In[r-siny]+10; Yyry=2:cigy. 

ug a . E -In[r - sin A =, 

2 a Intrsin a EHE 

08 12 Info sin ZH ne gt a l8r. 


Spannungsfunktionen für die Kugelzonenschale mit schubspannungsfreier 

sphärischer Begrenzung 

Wenn die sphärischen Grenzflächen r = r, und r = r, eines hohlen Kugelkegels mit den 
konischen Begrenzungsflächen y=%, und =, einen Abstand Ar =r,— r, haben, der im 
Vergleich mit r,, r, selbst klein ist, haben wir eine Kugelzonenschale vor uns. Wir wollen den 
Fallins Auge fassen, daß die beiden Kugelflächen schubspannungsfreisind und daß auf die Fläche 
r = r, auch keine Normalspannungen wirken, während wir das Gesetz der Verteilung der Nor- 
malspannungen über die äußere Kugelfläche r =r, uns noch vorbehalten. 


Ansatz ; 
Wir ziehen als Spannungsfunktion zunächst folgende Kombination heran: 
U=%k-U, kV, aD Do ea (6), 


wobei wir zur Vereinfachung der Schreibweise von den homogenen Parametern k,, ka, kz, k, von 
vornherein k, = 1 setzen wollen. Wir stellen aus hernach ersichtlichem Grunde die Bedingung 
voran, daß für die innere sphärische Begrenzung r =r, die Verteilung der Normalspannungen 


c, von der Winkelkoordinate y unabhängig sein soll. Das drückt sich aus durch eine Beziehung, 


durch die auch zwei mit derartigen Gliedern gekoppelte, an sich nur von r abhängende Glieder 


betroffen werden, nämlich durch 
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4 NER 2(m —2) 1 
Rt er... 0, . Rfard 
rag 7 8 =, nn woraus ka a; 
Soll U weiterhin für r =r, und r = r, Schubspannungen %,y = 0 liefern, so muß gelten: 
2 2(m—1) k %, 2 2(m—1) k k: 
ter oa he 0; —_ a > ur 
| at m—2. r, mt x ei m—2ı2 rn rn her 
oder 
2(m—N).k Kk 2 :2m—2) 1, 2(m—]1) u, k, 2 2m-—2) 1 
m—2 nn n m+1l rn’ 7 Pe ee ae mı+1‘ 7 
Hieraus ergeben sich die Parameterwerte \ 
IE a Im— N r 
> m— 2 > 2 (m — 2) de -2)|; n=—2-(} s) 2(m WER 
2ume—- 1) 7, tr, set Ti \Tk ME En Va N = 


Bei der durch diese Multiplikatoren gekennzeichneten Spannungsfunktion U fassen wir zunächst 
nur diejenigen Glieder in den Normalspannungen o, ins Auge, die von der Winkelkoordinate x 
nicht abhängen und auch nicht mit x-Gliedern gekoppelt sind, die also auf der Begrenzungsfläche 
r=r,alseinzige verblieben sind. Die entsprechenden Werte sind: 


, E (m — 97: N. 2 

aufop.=y.: en. ı ky ; 1 Ba, en | 
DE: ER Be 2 ent, mt 1 

4 (m — J% j B 9 

aufr=r;: hasst = im Br k; RR RS U | 
m+1 m— 2 a m-2.n. 5 m+1 


Wir können diesen Anteil der Oberflächenspannungen kompensieren durch Überlagerung von 
Lösungen A,, 4,, die ja nicht mit Schubspannungen verknüpft sind (s. oben). Zu diesem Zweck 
haben wir eine Spannungsfunktion K, = U -+a,: A, +a, A, anzusetzen, in der die Koeffi- 
zienten a], a, den Bedingungen genügen: 

mE nn 1 2%; GE 1 “ AmE a, mE 


ee BET Zus‘ wen ® 25 
m+1 m — 2 ER are m+1 m+1l rn m_ 20% 
ai mE nn li ı 12005 eg ee 4mE a, mE 2 
m+llm-2 2 m-2n nt m+l m+1 r!' m—2 ” 
Daraus ergibt sich: 
m+1l r:r m—2 L(r,)-r;, — L(r-)* r} 
N ER A BE DALE TU NET, 
EN Te a mE le 
wenn ) 
a ke Re, 3m +2 
DE 2m — 2. nur ml 


Die Kugelschale, auf die sich die so fixierte Spannungsfunktion K, bezieht, ist für die Grenz- 
fläche r =r, normal- und schubspannungsfrei. Die Grenzfläche r =r, ist schubspannungs- 
frei, erfährt aber Normalspannungen entsprechend den in U,, U, enthaltenen durch unseren 
Ansatz nicht kompensierten Gliedern, nämlich 


De ar a In[rı sin 24 De 2 In[r, sin y] = la a n) In[r, sin z]. 
Die Basstung der Kugeloberfläche zerfällt also in einen konstanten Bestandteil 

4nE [1 1 

= mr) a, 


und einen Bestandteil, der Insin x proportional, der also mit Annäherung an den Polxy=0 
logarithmisch divergiert. 
Das vollständige Ergebnis für das vorgelegte een piohlen K stellen wir nochmals 


zusammen: E 
Ä 2 m—4 4 3m —4 k m —4 
—_I— 4 2 5 ir rn a 7 v3 N 
= = kır)in abe sin x] m —2 E m— 2 har r au ER 
4=|7- a 7 Ir F- rieeigy 
ze 1 3m-+2 
SEE + rung! 


mE (2  2(m+1) 


: m 
- k.) In[» - sin y] 2 


B RT m+ıllr m — 2 2 r3 r3sin? x "RR 

E re rn ee 
6 = = ie + un. t) In [r sin y] — .: ; a: 2 .ctg?y + a \ Er 
’ te on ctg? a 
Et | 


Alle Größen divergieren für —0, so daß die Lösung a an sich nur für das Problem der Kuel: 
B. zonenschale y, Sy <S%,»rı Sr Sr, herangezogen werden kann. Aber man kann im Fall 

der Kugelhaube 0 <y < x, auf sie zurückgreifen, wenn es sinnvoll ist, aus dieserein koni- 
Be sches Kernstück 0 <y<y, auszusondern und durch Material von 
höherem Elastizitätsmodul zu ersetzen, das bei großer Deformation hohe 
h; Spannungen, insbesondere Scherspannungen 7,, aufzunehmen vermag. ‚Wir kommen hierauf 
u bei der Besprechung des Kuppelproblems zurück. 


Ansafzıll 
Wir ziehen als Spannungsfunktion folgende Kombination heran: 


Kr =hD, +1,D, +1D; +4uD,. | 
Von den homogenen Parametern I,» I,, 15, 1, können wir I, = 1 setzen. Wir wollen wieder die Be- a 
i dingungen einführen, daß die Kugeloberflächen r =r, und r =r, beide schubspannungsfrei A 
Br sind und daß die innere Kugelfläche r = r, auch normalspannungsfrei wird. Für den letzteren | 
Zweck haben wir, da den Funktionen D, und D, keine Radialspannungen 0; zugeordnet sind, 
einfach in der Spannungsfunktion D, + % D, den Parameter 1, zu bestimmen dureh, 


» AT m+1 RER 


een; also = 2m) 


Aus den Bedingungen dafür, daß auf die Flächenr =r,undr = nz keine Tangentialspannungen | 
_  „. wirken, folgen dann die en ar, et? 


I 
1 4» BEN =— 4 = — A u‘ > br: R 
S + 2rı 2 N Bi Eur 2 N ER RR 
Hieraus bestimmen sich die Parse 1; u. ie ER a Ma 2 Re FOR 
h; 81 MEERE NE 
N . Ik EA ra, SR ET a SPEER Tree A ee: 
sul De L,; lı 2, (nt Be Pur je N ar: 
Die damit festgelegte Kombinations- Spannungsfunktion Kır liefert: 
ai ei 4 Are 
ern) ne: 


2 mE ee Rn: mE 1 
Emltaraatai Sy mt 
E +! 3 cos x 
08 = m+1 rag ae sin? x 
mE 41 1 
127 = al - a4]. siny' iA 


Das ER Ergebnis ist, daß bei dem vorliegenden 
5 le Kugelfläche de nunsegei ist, während die 
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sind, die also für y=% verschwinden und für x—0 logarithmisch un- 
endlich groß werden. Durchweg divergieren die Verschiebungs- und die Spannungs- 
komponenten für „—0. Der Ansatz Il ist demnach an sich, ebenso wie Ansatz I, nur auf das 
Problem der Kugelzonenschale x; Sx S% rı Sr <r, anwendbar, läßt sich. aber auch durch 
Einführung eines Kernstückes 0 <y < x, aus Material von höherem Elastizitätsmodul (Kugel- 
kegelstumpf 0 <y <y,) für die Kugelhaube nutzbar machen. 


Kuppelproblem 
Bei den in der Technik verwendeten Kuppeln handelt es sich meistens um Kugelschalen» 
die auf der Außenseite achsensymmetrisch verteilte Lasten bzw. im Scheitel stärker konzentrierte 
Lasten aufzunehmen haben. Beispiele hierfür sind die in den Ansätzen I, II behandelten Be- 
lastungsfälle X}, Kr, die wir jetzt zu K=c, K, + Cı Kır kombinieren. Die auf die schub- 
spannungsfreie Außenfläche wirkenden Normalspannungen sind dabei gegeben durch Ausdrücke, 


die einerseits auf ihr konstant, andererseits durch Glieder von der Form c; :Insin x, ey Intg & 
gekennzeichnet sind, wo = 
mE 4 m+1 ) 2 mE | A| 

ne AL A I. EEE RN ER 
e N ER m—2 ka]; Mm Tm—2 i 

Von den ersteren sehen wir bei der jetzt zu erörternden Problemstellung ab, da sie durch bekannte 

einfache Lösungen des achsensymmetrischen Deformationsproblems, nämlich A,, As, erfaßt 

und bestimmten auf die Innen- und Außenfläche wirkenden konstanten Normaldrucken 7, Pa 

angepaßt werden können. Dann kommt in Frage auf r = r, die Normalspannung 


cin (2 cos SE Cr Incos + (67 + Cr) Insin z = cı In 2-+ (cr — Gyr) In cos x + 


+(& +6) Insin 2 ; 


Wählen wir 5, = — c;, so wird die Belastung durch das Gesetz o,—= cı In2+2ciIn cos £ 
gegeben. 


Wählen wir c, = cy, so gilt das Belastungsgesetz 


D|N 


%, = uin2+2- cr -Insinz. 


Im ersteren Falle ändert sich die normale Oberflächenspannung von c; In 2 bis auf O0, wenn x 
vom Wert O0 bis zum Wert wächst. Im zweiten Fall ändert sich, wenn x das gleiche Intervall 
durchläuft, die normale Oberflächenspannung von — © bis auf c;, (1—1In 2). In beiden Fällen 
bleibt aber die Divergenz der Ausdrücke für die Verschiebungen, und — im allgemeinen — auch 
die Spannungen im Falle x = 0 besteher, und das gilt natürlich auch für andere Kombinationen 
von K7, Kır; d.h. für ihre Verwirklichung ist notwendig, daß in einem kleinen Kuppenbereich 
0 <y<gı ein konisches Kernstück aus Material von höherem Elastizitätsmodul eingeführt 
wird. Die Erfassung einer solchen in der Kuppenlinie y = 0 zur Auswirkung kommenden kon- 
zentrierten Kraft Zist das Typische des Kuppelproblems. Wir wollen weiterhin, 
um praktisch vorkommenden Belastungen nahe zu kommen, den ersten der oben besprochenen 
Sonderfälle zugrunde legen, wo die Normalspannung o, auf der äußeren Kugelfläche von 0 gegen 
6; In 2 konvergiert, wenn x vom Wert $ gegen den Wert 0 strebt. 

Das hier supponierte Gesetz für die auf die äußere Kuppeloberfläche r = r, wirkenden 


2 

kurve heraus, die in y = 0 mit horizontaler Tangente beginnt und bis zu ihrem in y = erreich- 
‘ten Endpunkt in dem Polarkoordinatenbild durchweg konkav einwärts gekrümmt ist. Es dürfte 
- für diejenigen in der Praxis vorkommenden Belastungsfälle eine brauchbare Unterlage geben, wo 
- es sich nicht um gleichmäßigen (hydrostatischen) Druck auf r = r, handelt. Das Wesentliche, das 
in der Divergenz der Spannungen für y=0 zum Ausdruck kommt, ist die Erfassung der 
inder Kuppe zur Auswirkungkommenden konzentriertenn KraftZ. 
Es steht nichts im Wege, die verteilten Oberflächenspannungen 9, zu beseitigen durch Über- 
lagerung einer Spannungsfunktion, die durch Ansatz einer Reihe der regulären mit cosnx 
bzw. sin nx proportionalen Spannungsfunktionen entsteht, die also die Form der allgemeinen 
für den Polarkoordinatenansatz geeigneten biharmonischen Funktionen haben: # = r?-9*+9, 


Normalspannungen 0, = c;In2+2cyIn cos = c; In (2 cos? kommt auf eine Belastungs- 


wo @,9* harmonische Funktionen des Typs V, =": P„ (cos xy); V;, = - P„ (cos x), unter 
P„(cos x) Kugelfunktionen erster Art verstanden ®). 


} 8) Vgl. A. Timpe, a.a.O., 8. 165. 


yrtl 


funktion ein: 


g Proceedings of the London Math. Soc. 32 (1901), 
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Kugelkegel- und Hohlkugelkegel-Probleme N: 

Bei den die Kugelschale betreffenden Spannungsproblemen sind wir besonders auf den 
Fall einer im Scheitel sich auswirkenden konzentrierten Kraft Z eingegangen und haben dabei 
den Gedanken herangezogen, den Bereich 0 <xy <x, durch einen entsprechenden Kugelkegel- 
stumpf aus anderem Material zu ersetzen, der die der Kuppe zuzumutenden außerordentlichen 
Spannungen aufzunehmen hat. Wir wollen jetzt das Spannungsproblem für diesen Kugelkegel- 
stumpf mit der Maßgabe ins Auge fassen, daß es angesichts der Kleinheit der Dicke r, —r, der 
Kugelschale nur darauf ankommt, denResultierenden der genannten Spannungen ge- 
recht zu werden. Wir bilden zu diesem Zwecke zunächst eine Kombination 8 der Spannungs- | 
funktionen A,, As, Bj, Ba, B;, T,, die darin übereinstimmen, daß sie im Bereich 0 Sy s%ı i 
reguläre Spannungskomponenten ergeben. Die sechs Multiplikatoren, mit denen diese Spannungs- 
funktionen eingeführt werden, unterwerfen wir den Bedingungen: Die aus B,, B,, B; resul- 
tierenden Normalspannungen o, sollen für r=r, und r =r, sich annullieren, so daß die der 
Gesamtkombination S entsprechenden o, auf den begrenzenden Kugelflächen jedenfalls konstant 
ausfallen. Die dritte und die vierte Bedingung gehen dahin, daß die o, auf r = r, verschwinden 4 
und aufr = r,als Resultierende die Kraft Zliefern. Damit wird aus statischen Gründen zugleich 
der Bedingung genügt, daß die Resultierende aus den tangentialen Randspannungen 7,, dieser 
Kraft entgegengesetzt ist. Die 5. und die 6. Bedingung gehen dahin, daß die auf dem Kege- 
stumpfmantel ringförmig verteilten, aus den Spannungen o, sich ergebenden Resultierenden und $ 
resultierende Momente übereinstimmen mit den entsprechenden Größen für die die Mutterrolle 
spielende Kugelschale. 

Weiterhin seien die Abmessungen auf den Kegplmantellinien vergleichbar mit den sphäri- 
schen. Man kann dann die bei der soeben behandelten Aufgabe herangezogenen Funktionen 
B,u. T, so verknüpfen, daß sich de SpannungsverteilungZ, in dem Conus 
ergibt,in dessen Spitze 0 eine konzentrierte Kraft Zin Achsenrichtung angreiftund dessen 
Mantely=x, spannungsfreiist. Setzt man L, = B,+%kT,, so wird 


a ME} O8 N ee er 
Be Te m+1 r?2 1+cosx j 
und A 
a ER mE. sing 7 m en h 
KT ML Nr m+1 r?2 1+cosy’ 
R k 


und beide Größen verschwinden für x = x, wenn 1+ 


BE k=—(l+cosg. 
| I ERS 0, also | (1+cosxo). 
Mit der Spannungsfunktion Z, = B,— (1 +cosy,): T, ergibt sich 9: # 


er 17%Bm—4. 1+cosxol .... 
m rt lt; a nn 
__mE 1 Be mE 1 ar ne SE 
r m-t1 72 ZT 00sy—tltesgol; 0% ee 
mE A en mE 1 1i+cosy.]2 2% Rn 
0 = —— : — |— cos y— 01. vr a NT Re er 
EEE RN IFeosy.’ Fey aan "| el ld: 


Die konzentrierte Kraft in der Kegelspitze ergibt sich durch Integration über die Spannungen, iu, 
die auf die von der Kegelfläche begrenzte Kugelhaube wirken, zu RT RE ® 


x 


m (m —1) 
& N (m +1) m —2) SR A 
Wir vermerken, daß die auf die Kugelfläche r — konst. wirkenden Spannungen dem Quadrat 
des Radius umgekehrt proportional sind. | Nein 2 a N ee 

‚.Es handele sich weiterhin um den durch 4,<x<S%s gekennzeichneten Hohlkuge 
bei dem die Mantelflächen y=yx, und x = x frei von Normal- und Schubspan 
Wir ziehen dann außer den oben benutzten Spannungsfunktionen B,, Tı noch 
Fall der Kugelzonenschale verwendeten) Funktionen D,, U, heran und führen : 

RN 


ZN N 


I;=aB, +bT,+0D,+dD,. 
mE Ku Bi. SE 
m+1 1a = wır 


°) Von anderen Ansätzen aus wurde die gleiche Sp: r 
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Bei Weglassung des Faktors 


s . 
« PD > ERS; y 
a 4 JEIR h 
r A Eu} .r CH y I 
ah AURR x 
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1 ( | cos x cos xy 
0, = — |—a:cosy—b-: ——— . Y4. to2 
AR e 1+cosy " siny a .) i 
Ft. inx sin x ) 
up= — a-siny—b- ——— wet d.cte 
RN | x 0.2 sin? x cIEX 


b c d-cosx 
—a— ——+— I #=0 
l+(05xy sin?y sın? x 


Formuliert man sie für y=x, und x = x, so erhält man zwei lineare Gleichungen für die vier 
Parameter a, b, c, d bzw. für-die maßgeblichen drei Verhältniszahlen, so daß über eine der Ver- 
hältniszahlen noch willkürlich verfügt werden kann. Das kommt darauf hinaus, daß man sich 
noch verschiedenen Situationen bezüglich der Oberflächenbedingungen für die Spannungen oder 
Verschiebungen anpassen kann. Man hat: 


4(m+-1)(m —1) cosyx m+t1 1 m+1 3m —4 sin 
re. ee es hy ar ER a ee 
I m(m — 2) E us ) mE r’ “x mE m—2 Pi 
Br, | 2(2 m—1) | 1 cos x cosx R 
2 er ne u er m er; 
Be = |—a:cos N oy., .d | 
ep x 1+cosyx sın?y sin2yl’ 
1 4 sinx PAR, 
| ey en R 


Oft wird die Sachlage dadurch gekennzeichnet sein, daß die auf die beiden begrenzenden Kugel- 
flächen wirkenden Radialspannungen von der Winkelkoordinatey un- 
abhängig sind, was der Festsetzung «a = 0 entspricht. In diesem Falle geht die Bedingung 
dafür, daß eine Mantelfläche x = konst. spannungsfrei ist, über in 

b c d-cosyx 


eo LEE ABA EAST A 
2° cos 9 4 sin 9 cos 5 4sin 5 co 9 


oder wenn von den homogenen Parametern b =1 gesetzt wird: 


ed, 


0—d-c0sy—2sin? &. 


Die Anwendung auf die beiden Mantelflächen y=xı und x =xzdes hohlen Kegelstumpfes 
liefert für die verbliebenen Parameter die in die Spannungstunktion 7T,+cD,;,+d U, einzu- 
setzenden Werte ; 


FR 4 } 
cos Yı° sin? 2 — cos X," sin®— sin?2 _ sin? #1 { 
2 2 2 
G= 2: ———; d=2:- 


cos X — C08 X COS Xı —C0sXg 


Damit hat man die Lösung für das Problem des hohlen Kugelkegels, in dessen Spitze eine kon- 


zentrierte Kraft angreift und dessen Mantelflächen x = %,, X = x3 Spannungsfrei sind, während 
die auf die sphärische Begrenzung wirkenden Spannungen dem Quadrat des Radius umgekehrt 
proportional sind, oder, was auf dasselbe herauskommt, die entfernten Kegelteile in Ruhe ge- 
halten werden. Die Spannungsfunktionen D;,, U, erweisen sich als das für den Hohlkegel wesent- 


‚liche Gegenstück zu der für den Kegel charakteristischen Kombination der vnBoussinesq 
- aufgezeigten Spannungsverteilungen, die als erster und zweiter Typ einfacher Lösungen bekannt 
‚sind. 


Zur Behandlung von Aufgaben für Kegel und Hohlkegel, bei denen die Mantelspannungen 
Gy, Try als Funktionen von r gegeben sind, sind Reihenentwicklungen heranzuziehen, die sich 


- auf die durch (4) und (5) erklärten biharmonischen Funktionen stützen. Die erfolgreiche Ver- 
_ wendung von Spannungsfunktionen bei Zylinder- und Polarkoordinaten legt den Gedanken 
_ nahe, in allgemeineren Fällen krummliniger Koordinaten nach dem entsprechenden Gegenstück 


— einer von einer dieser Koordinaten unabhängigen erzeugenden Funktion — zu forschen. 
3 


\ x k 
‘ (kr y fen: 7 i Ar RT 
‘2 x 4 ac Wal > 
LE f BT Rn . 5 N WE a 
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Ein Polarplanimeter zur Bestimmung des po- 
laren Trägheitsmomentes. 

Definiert man die Stellung des Systems Polarm— 
Fahrarm in bezug auf ein rechtwinkliges Koordinaten- 
system durch die Winkel p und y gemäß Bild 1, so 
kann man diese Winkel als neue Koordinaten in dem 
vom Fahrstift überstrichenen Bereich auffassen. 

Es gilt dann 


2=NCo8sp-L Teos(a ty) : : "E (1); 

y=psinpg+Iisin(P+Y) 

typen ia 2 plcosy - (2), 
dx = [—psinp—Isin (p + y)ldp 
r —Isin(p+ y)dy (3) 
{ dy=[pcosp +lcos(p + y)ldp ; 
€ +Tcos(p + y)dy 
5. 
hr 
ie. 
f yt--------#----7- 
We \ 

ns 
N | pzconst. 


0 


Bild 1 


Die Funktionaldeterminante wird 


du da 


ua so. ket'ya 


Die Fläche F wird demnach 
F=f[fdady=plf/snydopdy: - - (5). 
Führt man bei konstantem p die Integration über y 
aus, so wird 
F=—pl/f (cosy, — cos y,) dp. 
%,, Y, sind. die y-Werte auf der Randkurve von F. Die 
weitere Zusammenfassung liefert für die Fläche F das 
Kurvenintegral 
; F=—plßcosy do 
Liegt der Pol innen, betrachtet man zunächst die 


Differenzfläche zwischen der vorgegebenen Kurve und 
einem Kreis y, = const. 


TED BE ru) 


Wählt man y, = 3: so wird nach (2) der Radius des 


Kreises gleich Yp* +1? und die Differenzfläche wegen 
cos y, = 0 wieder gleich dem Integral (6). Die Gesamt- 
fläche ergibt sich durch Hinzufügen der Kreisfläche zu 


F=—pl$cosy dp+r(p? +12) (Polinnen) (6a). 
Für das polare Trägheitsmoment in bezug auf den 
Pol (Nullpunkt) ergibt sich mit (2) 


Jp=/f(+y?)dady=fSf(p + +2plcosy)dady 
=(e+PR)F+2pXlffcosysinydydp 


2 ]2 
= (pt +) F— 5 (cos 29, — cos 2y,)dp 


dh . 212 5 
p=(pP +B P— 5 beos2ydp - : (M. 


Liegt der Pol innen, so findet man in derselben Weise 
wie bei der Flächenbestimmung 


v ae 
„=(pP+B) F—5 Br) — u; cos2y dp (7a). 


Er / A 
u P A 
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Die in (6) und (7) auftretenden Kurvenintegrale lassen 
sich nun mit Hilfe von Meßrollen mechanisch be- 
stimmen, die mittels ‚‚Bollenarmen“‘ ebenfalls mit dem 
Polarm verbunden sind. Die entsprechenden, Größen 


Bild 2 


seien p7, Ip und yz (vgl. Bild 2). Außerdem möge die 
Rollenachse mit dem Rollenarm den festen Winkel & 
einschließen. Bezeichnet man die Anzahl der Rollen- 
umdrehungen mit U, den Rollenhalbmesser mit o, 
das Bogenelement der von der Rolle durchfahrenen 
Kurve mit ds, sowie den Winkel zwischen Rollenachse 
und Kurventangente mit ö, so wird 
2rodU =sin öds 8), 
da die Rolle nur die zur Rollenachse senkrechte Kom- 
ponente des Bogenelementes ds mißt. Nun ist aber 
6=ptyrtema (9), 
wenn man mit & den Neigungswinkel der Kurven- 
tangente bezeichnet. Aus (8) wird somit 
2rodU=sin(p+yr+ e)cosads 
— cos(p+yr+e)sinads 
=sin (pt yYyR+E)de—cos(p+yrR + E)dy 
= -PrCos(yr+ e)dp-IpRcosedp 
 —Ipecosedyn, 


wenn man dx, dy gemäß (3) einsetzt. Integriert man 
jetzt über eine geschlossene Kurve, so ergibt sich 


2re(U,-U)=-Pr$ cos(yr + e) dp h (10) 


— 10088 (95-9,) 
wobei 9, — 9, =0 bzw. —2r ist, je nachdem der Pol 
außen oder innen liegt. 

Setzt man jetzt 
vrte=y: 0: : (11), 
so stimmen die Kurvenintegrale in (10) und (6) überein. 
Die Beziehung (11) bedeutet nun, daß Rollenachse und 
Fahrarm parallel sind, da ja y5 + & der Winkel zwi- 


schen Bollenachse und Polarm ist. Durch weitere 
Umformung von (10) erhält man die üblichen Plani- 
meterbeziehungen. 

Um nun das Kurvenintegral in Formel (7) bzw. (7a) 
zur Bestimmung des polaren Trägheitsmomentes zu 
erhalten, muß für die entsprechende Meßrolle 

Yr+te=2y (12) 
gesetzt werden. Das bedeutet aber im wesentlichen, 
daß die Rollenachse um einen doppelt so großen Winkel 
gegenüber dem Polarm ausgelenkt wird wie der Fahr- 
arm. Mechanisch läßt Sich das in derselben We 
bewerkstelligen wie beim Potenzplanimeter. In der 
Nullstellung (y= 0) sind Polarm, Fahrarm sowie 
Bollenachse parallel. PETE 

Es sei in diesem Zusammenhang noch ds 
gewiesen, daß sich das Flächenmoment : 
Weise durch die Differenz 
momente, bezogen auf zwei 
drücken läßt. . 
Berlin. 


“ie, No Kae ee! 
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BUCHBESPRECHUNGEN 


Arno Ziegler, Einfache Formeln für 
Träger über 2—10 Felder. 99 $S. mit 48 Bei- 
spielen und zahlreichen Abbildungen. Bonn 1949. 
Verlag Ferd. Dümmler. Preis brosch. 7,80 DM. . 

Das Buch enthält eine Zusammenstellung der 
- Stützenmomente durchlaufender Balkenträger mit 
2—10 Feldern und konstantem Trägheitsmoment für 
die üblichen Belastungsfälle. Der Verfasser wieder- 
holt damit Angaben, die sich in zahlreichen Lehr- und 
Taschenbüchern finden. Die Veröffentlichung würde 
also nur eine Berechtigung haben, wenn die Aus- 
wertung der Rechenergebnisse oder ihre Darstellung 
eigenartig ist. Beides trifft nicht zu. Das Buch ge- 
hört zu den zahlreichen Veröffentlichungen der tech- 
nischen Literatur, die dem Leser Rezepte ohne Be- 
gründung geben. Einzelne Unstimmigkeiten dürfen 
bei der Anwendung nicht übersehen werden. 


Dresden. KR. Beyer. 


Dr. Lothar Heffter (Prof. an der Universität Frei- 
burg i. Br... Kurvenintegrale und Be- 
gründungderFunktionentheorie.48S8. 
mit 7 Abb. Berlin—Göttingen—Heidelberg 1948. 
Springer-Verlag. Preis brosch. 5,40 DM. 

Dieses ‚‚Lehrheft‘ stellt eine abschließende Zusam- 
menfassung von Gedankengängen dar, die der Ver- 
fasser in einer ganzen Anzahl von Arbeiten in den 
Jahren 1902—1941 veröffentlicht hat. Es handelt sich 
dabei vor allem um die möglichst einfache Formu- 
lierung ‚einer hinreichenden und notwendigen Bedin- 
gung dafür, daß /(z) = u(2,y) + iv(zy)ß= 2 + iy] 
in einem Bereich @ analytisch, d.h. durch Potenz- 
reihen darstellbar ist. Im Gegensatz zu dem in den 
Lehrbüchern üblichen Vorgehen, den Cauchyschen 
Integralsatz aus dem analytischen Verhalten von /(z) 
herzuleiten, macht Verf. die eindeutige Integrierbar- 
keit zum Ausgangspunkt der Untersuchung und ge- 
langt so zu einem Satz, der in der Form der Aussage 
dem Satz von Morera ähnelt, in seinen Vorausset- 
zungen jedoch weniger fordert. Für eine in @ stetige 
Funktion /(z), für die jedes Integral, genommen über 
ein beliebiges in @ gelegenes Rechteck (mit achsen- 
parallelen Seiten), verschwindet, wird mittels eines 
Hilfssatzes die Cauchysche Integralformel bewiesen und 
daraus die Potenzreihendarstellung abgeleitet. Durch 
Beschränkung auf die achsenparallele Integrierbarkeit 
ist der Begriff des allgemeinen Kurvenintegrals ent- 
behrlich geworden, und.damit hat das Ergebnis einen 
Grad von Einfachheit erlangt, der es geeignet macht, 
neben der Cauchyschen und Weierstraßschen Auffas- 
sung als dritter Ausgangspunkt der Funktionentheorie 


zu dienen. Allerdings kann sich Bef., gerade auch auf - 


Grund der vorliegenden Arbeit, nicht der Ansicht des 
Verf. anschließen, daß die Integrierbarkeit eine schwä- 
chere Bedingung für /(z) bedeutet als die Differenzier- 
barkeit oder Entwickelbarkeit, eine Tatsache, die wohl 
nur zutrifft, wenn man eine einzelne reelle Funk- 
tion betrachtet, nicht aber, wenn zwei reelle Funk- 
tionen u und v in der hier angegebenen Weise ge- 
koppelt sind und damit eine komplexe Funktion 
bilden. 

In einem einleitenden Abschnitt sind, zum Teil ohne 
Beweis, die Grundlagen aus der Theorie der reellen 
Funktionen dargestellt, soweit sie benötigt werden, 
so daß die Arbeit auch schon von Studenten höherer 
Semester gelesen werden kann. Außerdem sind noch 
Abschnitte iiber Beziehungen von Kurvenintegralen 
einmal zu Treprenintegralen und zum anderen zu 
Stieltjes-Integralen mit stetiger Belegungsfunktion hei- 
gefügt, die im Zusammenhang mit dem Hauptthema 
stehen. 


Dresden. G. Opitz. 
Dr. R, Rothe (ehem. Prof. a. d. Techn. Hochschule 
Berlin, Höhere Mathematikfür Mathe- 
matiker, Physiker und Ingenieure. 
Teil III: Flächenim Raume. Linien- 
integrale und mehrfache Integrale. 
Gewöhnliche und partielle Diffe- 
rentialgleichungen nebst Anwendun- 
gen. (Teubners Mathematische Leitfäden , Band 23.) 
4. Aufl. 236 S. mit 167 Abb. Leipzig 1949. Verlags- 
gesellschaft B. G. Teubner. Preis kart. 7,— DM. 


Teil IV, Heft 5/6: (Unter Mitwirkung von Studien- 
rat OÖ. Degosang) Ubungsaufgaben und 
Lösungen zu Teil III. (Teubners Mathe- 
matische Leitfäden Band 37/38.) 3. Aufl. 108 S. mit 
59 Abb. Leipzig 1949. Verlagsgesellschaft B. G. Teub- 
ner. Preis kart. 3,60 DM. 

Mit dem jetzt in vierter Auflage vorliegenden 
Band III und der dritten Auflage des Heftes 5/6 des 
Bandes IV, das den dazu gehörenden Übungsstoff 
bringt, ist die Neuauflage des ausgezeichneten Rothe- 
schen Unterrichtswerkes bis auf dıe Formelsammlung 
abgeschlossen. Auch diese beiden Bände sind ein fast 
unveränderter Abdruck der vorigen noch von Rothe 
selbst durchgesehenen Auflage. Eine besondere 
Empfehlung der überall geschätzten Bände von Rothe, 
die zu den besten deutschen Unterrichtswerken zählen, 
erübrigt sich. Da der Preis im Vergleich zu ähnlich 
inhaltsreichen Werken mäßig ist, darf man hoffen, 
daß auch diese Auflage in viele Hände kommen wird. 

Dresden. : Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


Prof, B, v. Borries, Die Übermikroskopie. 
Einführung, Untersuchung ihrer Grenzen, und Abriß 
ihrer Ergebnisse. 416 S. mit 225 Abb. Berlin 1949. 
Verlag Dr. Werner Saenger. Preis geb. 48,— DM. 


Proceedings of Symposia in Applied Mathematics, 
Vol. I: Non-Linear Problems in Mechanics of 
Continua. VII + 219S. New York 1949. American 
Mathematical Society. Preis geb. 5,25 $. 


Dr. habil. F. Ringleb, Mathematische Formel- 
sammlung. (Sammlung Göschen Bd.51.) 5. ver- 
besserte Auflage. 2748. mit 37 Abb. Berlin 1949. 
Verlag Walter de Gruyter &Co. Preis brosch. 2,40 DM. 


Dr. Fritz Chemnitius, Kleines Handbuch 
der Infinitesimalreehnung. VIII+1848. 


14,— DM. 


r 


Pößneck 1949. Rudolf A. Lang-Verlag. Preis geb. 


F. Mühlig, Der 24 m-Interferenzkomparator 
des Geodätischen Institutes in Potsdam. 
(Veröffentlichungen des Geodätischen Institutes in 
Potsdam Nr.2.) 508. mit 30 Abb. Berlin 1949. 
Akademie-Verlag. Preis brosch. 7,— DM. 


Dr.-Ing. Werner Matz, Die Thermodyna- 
mik des Wärme- und Stoffaustau- 
schesin der Verfahrentechnik. XVI+ 
355 S. mit 114 Abb. Frankfurt/Main 1949. Verlag Dr. 
an Steinkopff. Preis brosch. 26,— DM, geb. 
28,— DM. 


Dr. Carl Ramsauer (o. Prof. a. d. Techn. Univ. Ber- 
lin, Physik—Technik— Pädagogik(Er- 
fahrungen und Erinnerungen) (Wissen- 
schaftliche Bücherei) IV - 130 S. mit 4 Abb. Karls- 
ruhe 1949. Verlag G. Braun. Preis brosch. 4,20 DM. 
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x 
Dr. H. Hasse (o. Prof. a. d. Univ. Berlin, Direktor am 
Forschungsinstitut für Mathematik der Deutschen Aka- 
demie d. Wissenschaften zu Berlin). Zahlentheorie. 
XII + 470 8. Berlin 1949. Akademie-Verlag. Preis 
brosch. 44,— DM., geb. 47,— DM. 
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“N Tagung der deutschen Mathematiker-Ver- 
einıgung in Köln. 


Die Dautsche Mathematiker-Vereinigung hielt vom 
18.—23. S>ptember ihre diesjährige Tagung in Köln ab, 
an der eıne große Anzahl von Mathematikern nicht nur 
der Westzonen, sondern auch der russisch besetzten 
Zone und von Berlın, sowie einige Ausländer teilnahmen. 
Inder E.öffnungsansprache wies der Vorsitzende, Prof. 
D:, Kamke, aut dıe Schwierigkeiten hin, die für die 
mathematische Forschung durch die erheblichen Lücken 
in der ausländ.schen L.teratur an unseren B.blıotheken 
entstehen, sowie auf dıe Gefahr, die in der Abwanderung 
vieler jüngerer Forscher in das Ausland besteht. Er 
bat dıe Kultusministerien, das Ihrige dazu zu tun, daß 
die Lücken in den B.bliotheken schnell aufgefüllt wer- 
den und daß den deutschen Hochschulen ıhr wissen- 
schaftlıcher Nachwuchs erhalten bleibt. 

Nach einem einleitenden Vortrag von ©, Müller- 
Hannover begann dıe w.ssenschaftliche Tagung mit 
Vorträgen über Grundlagenfragen (Kneser-Tübıngen, 
Lorenzen-Bonn, Martin-Köln, Schröter-Berlin). Gıup- 
pentheorie und moderne Algebra stehen nach wie vor 
im Vordergrund des Interesses, wie die Vorträge von 
Keller-Dresden, Krull-Bonn, Pickert-Tübingen, Schö- 
neborn-Bonn, Wever-Göttingen zeigten. In einer Reihe 
von Vorträgen wurden Fragen aus folgenden Gebieten 
behandelt: Theorie der Funktionen einer und mehrerer 
komplexer Veränderlicher (Cremer-Aachen, Grunsky- 
Tüb.ngen, Lammel-Tutzing, Sommer-Münster, Stein- 
Münster, Thimm-Bonn), Geometrie, insbesondere D.ffe- 
rentialgeometrie (Herrmann-Wolfenbüttel, Rembs-Ber- 
lin, R.tter-Berlin-Charlottenburg, Sparner-Freiburg, 
Strubecker-Karlsruhe),Summierungsverfahren(Knopp- 
Tüb.ngen, Meyer-König-Stuttgart), D.fferential- und 
Integralgleichungen einschließiich Potential- und Ope- 
ratorentheorie (Hamburger-Ankara, Maruhn-Jena, 
Mönnig-Aachen, Sauer-Frankfurt, Schäfke-Berlin, H. 
L. Schmid-Berlin, A. Schmidt-Jena, H. Schmidt- 
Braunschweig, Ullrich-Kiel, Wielandt-Mainz), An- 
gewandte Mathematik (Collatz-Hannover, Marre- 
Ehingen, Unger-Darmstadt). Walther-Darmstadt be- 
richtete über die neuere Entwicklung der Rechenauto- 
maten in Amerika und England, dıe dort in vollem 
Gang ist und für die große finanzielle Mittel zur Ver- 
fügung stehen, eine Entwicklung, an der Deutschland 
wegen Mangels an Mitteln zur Zeit nicht aktiv teil- 
nehmen kann. Außerdem wurde über eine Reihe von 
E.nnzelfragen vo getragen (Belck-Frankfurt, Dörge- 
Köln, Jaeger-Göttingen, Kappos-Erlangen, Lorey- 
Frankfurt, Noack-Köln, Richter-Haltingen, Schneider- 
un Göttingen, Svenson-Erlangen, Töpfer-Solingen, Witt- 
Hamburg). Auffallend war, daß kein Vortrag zahlen- 
FE theoretischer Natur im eigentlichen S:.nne war. Im 
- ganzen g’nommen fiel das hohe Niveau der Vorträge 
auf. Auch eine Reihe von jüngeren Mathematikern 
lieferten wertvolle Beiträge zu der Tagung. 

In der Mitgliederversammlung im engeren Sinne 
wurde unter anderem über Referatenorgane, den Jah- 
resbericht der DMV und die Diplomprüfung gesprochen. 


GmbH., Berlin NW 7, Schiffbauerdamm 19, 


c 
a 


Nachrichten 
w h 


Jos. E. Hofmann, Die Entwicklungsgeschichte 
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der Leibnizschen Mathematik während des 

Aufenthaltesin Paris (1672—1676). VIII + 2528. 
mit 27 Abb. München 1949. Leibniz-Verlag. Preis 

geb. 26,— DM. ef | 


NACHRICHTEN 


Dem letzten Thema wurde ein großer Teil der Ver- 
handlungen gewidmet, mit dem Ziel, bei der erforder- 
lichen Umarbeitung der Diplomprüfungsordnung des 
Jahres 1942 möglichst zu einheitlichen Anforderungen 
zu kommen. An diesen Verhandlungen nahm auch Mi- 
nisterialrat Gedigk von dem Kultusministerium Nord- 
rhein-Westfalen teil. Kamke, 


Frühjahrstagung der Gesellschaft für angewandte Ma- 
thematik und Mechanik. 


Aus verschiedenen Gründen erschien es zweckmäßig, 
im Herbst 1949 keine Tagung der Gesellschaft für an- 
gewandte Mathematik und Mechanik zu veranstalten, 
‘sondern sie auf das Frühjahr 1950 zu verschieben. Die 
"Tagung sollin Darmstadt nach Ostern 1950 stattfinden, 
und zwar von Montag, 17. April bis Mittwoch, 19. April. 
Anreise und Begrüßungsabend Sonntag, 16. April. Ge- 
schäftssitzung der beiden Zweige der Gesellschaft 
Dienstag, 18. April vormittags; abends voraussic htlich 
Empfang durch die Stadt. Ein Tag ist für moderne 
Rechenhilfsmittel (Integrieranlagen und Rechenauto- 
‚maten) und Tabellierungen in Aussicht genommen, 
‘wobei Professor D. R. Hartree (Cambridge) und 
Dr. J. C. P. Miller (London) Referate über auslän- 
dische Entwicklungen erstatten werden. Es wird ge- 
beten, schon jetzt Anmeldungen (insbesondere für Vor- 
träge), Quaitierwünsche, Nachfragen u. dgl. an den 
örtlichen Tagungsleiter, Herın Pıcf.Dr.A.Walther 
(16) Darmstadt, Fichtestraße 32, Telefon 936, zu 
richten. R. Grammel. 


Hannover: Herr Prof. Dr.-Ing. Flachs- 
barth, seit 1947 Rektor der Techn. Hochschule 
Hannover, wurde für die Amtszeit 1949/50 zum Rektor 
wiedergewählt. nr 


Frankfurta.M.: DerPrivatdozent der physika- 
lischen Chemie (Anwendungstechnik) Dr. Richard 
Klar an der Univ. Frankfurt a. M. wurde am 31.8. 
1949 zum a. o. Professor ernannt. _ } 3 

Dem o. Prof. an der ungarischen Universität in 
Oluj/Rumänien, Dr. 2 
Ordınariatfür Mathematik an der Techn. Hochschule 
Miskole/Ungarn übertragen. Br ve 


Druckfehlerberiehtigng 
0. Emersleben: Die Schwingungsdauer eines um- 


laufenden Pendels als Analogie zum Potential eins 


Kreises. Z.arg°w. Math. Mech. 29 (1949) Ba 3 
* stat 


S.279 in Fußnote 4 heißt es: „Cl. Schaefer“ : att 
„C. L. Schaefer“. IE: VAR, R ae 
S.280 (Formel 13a) im Argument des Tangen 
Zu 0 


n n N 
87T statt rn heißen. u SE 
S.2s0 (Formel 17c): Der Zähler 4 und d 

strich über &, sind ausgefallen. 
S. 281 (Formel 23a): Der Punkt zwis 
der Klammer ist zu streichen. 
2 . 
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